十進BASIC 第2掲示板 過去ログ

LET N=8

PUBLIC NUMERIC D,ws
LET D=1.432
LET L=5e-3
LET wd=0.3e-3
LET ws=L/(N-1)

!     /2π
!p(h)=| {1/y1(x,h)-1/y2(x,h)}*dx
!     /0

LET ndiv=100000 !計算分割数

PUBLIC NUMERIC h
LET h=wd/2
LET V1=Simpson(0,2*PI,ndiv)
LET h=ws-wd/2
LET V2=Simpson(0,2*PI,ndiv) !理屈上はV1=-V2になる。


LET Capa=1/(V1-V2)

PRINT "V1=";V1,"V2=";V2
PRINT "Capa=";Capa

END


EXTERNAL FUNCTION f(x)
LET y1=SQR(h^2+(D*SIN(x/2))^2)
LET y2=SQR((ws-h)^2+(D*SIN(x/2))^2)
LET f=1/y1-1/y2
END FUNCTION


EXTERNAL FUNCTION Simpson(a,b,N) !シンプソン法 1/3則
LET HH=(b-a)/(2*N)
LET S1=0
FOR i=1 TO N !h/3*{f(a) + 4*Σ[i=1,n]f(a+(2*i-1)*h) + 2*Σ[i=1,n-1]f(a+(2*i)*h) + f(b)}
   LET x=a+HH*(2*i-1)
   LET S1=S1+f(x)
NEXT i
LET S2=0
FOR i=1 TO N-1
   LET x=a+HH*(2*i)
   LET S2=S2+f(x)
NEXT i
LET Simpson=HH*(f(a)+f(b)+4*S1+2*S2)/3
END FUNCTION

　EXTERNAL FUNCTION p(h)  !シンプソン公式で積分
　DEF y1(x)=SQR(h^2+(D*SIN(x/2))^2)
　DEF y2(x)=SQR((ws-h)^2+(D*SIN(x/2))^2)
　LET dx=2*PI/ndiv
　LET sump=0
　FOR n=0 TO ndiv-2 STEP 2
　   LET x0=n*dx
　   LET x1=(n+1)*dx
　   LET x2=(n+2)*dx
　   LET fy1=1/y1(x0)+4*1/y1(x1)+1/y1(x2)
　   LET fy2=1/y2(x0)+4*1/y2(x1)+1/y2(x2)
　   LET sump=sump+fy1-fy2 !※←←←←←
　NEXT N
　LET p=sump*dx/3
　END FUNCTION

　EXTERNAL FUNCTION p(h)  !シンプソン公式で積分
　DEF y1(x)=SQR(h^2+(D*SIN(x/2))^2)
　DEF y2(x)=SQR((ws-h)^2+(D*SIN(x/2))^2)
　LET dx=2*PI/ndiv
　LET s1=0
　LET s2=0
　FOR n=0 TO ndiv-2 STEP 2
　   LET x0=n*dx
　   LET x1=(n+1)*dx
　   LET x2=(n+2)*dx
　   LET s1=s1 + 1/y1(x0)+4*1/y1(x1)+1/y1(x2)
　   LET s2=s2 + 1/y2(x0)+4*1/y2(x1)+1/y2(x2)
　NEXT N
　LET p=(s1-s2)*dx/3
　END FUNCTION