Print"Rs(単位 Ω)=(=γ*RtL)=",Rs Print"Cs(単位 F)=(=θ*CtL)=",Cs Input "CtLを入力してください(単位 F)．",CtL Input "RtLを入力してください(単位 Ω)．",RtL Print "Emeas/V" For t = 1 to 1300 n=0 While n <=100 n=n+1 Let a=n*PI() Let c=Cs/CtL DEF f(x)=TAN(x)-c/x DEF g(x)=COS(x)^(-2)+c*x^(-2) LET x=a LET E=0.000001 DO LET x1=x-f(x)/g(x) IF ABS(x1-x)<E THEN EXIT DO LET x=x1 LOOP S=0 S=S+cos(x1)*EXP(-t*X1^2/RtL/CtL)/(x1*((CtL+Cs)*sin(x1)+x1*CtL*cos(x1))) Wend Emeas=Ei+1/1000*(Iapp*(Rs+2*t/(2*CtL+Cs)+RtL*CtL^2*(8*CtL+5*Cs)/(6*(2*CtL+Cs)^2)-2*RtL^2*S)) Print a Next t END
