Rx = 8.880E+02 ! [Ohm] Lx = 2.820E+00 ! [Henry] Cx = 4.700E-05 ! [Farad] *afu,rad PI = ACOS(-1) p1 = -Rx/2/Lx+SQRT(Rx**2/4/(Lx**2)-1/Lx/Cx) $p2 = -Rx/2/Lx-SQRT(Rx**2/4/(Lx**2)-1/Lx/Cx) U = 4.000E+03 ! [Volt] l = 3.000E-02 ! [second] q = 6.000E-03 ! [second] prah=2000 *DIM,AK,,500 *DIM,BK,,500 *DIM,FK,,500 *DO,I,1,500 k = I+prah x1 = PI**2*k**2+l**2*p1**2 $x2 = PI**2*k**2+l**2*p2**2 t = q akq = U*(l*(exp(p1*t)*p1*x2-exp(p2*t)*p2*x1)*COS(k*PI*t/l)-k*PI*(exp(p2*t)*x1-exp(p1*t)*x2)*SIN(k*PI*t/l)) bkq = U*(l*(exp(p1*t)*p1*x2-exp(p2*t)*p2*x1)*SIN(k*PI*t/l)+k*PI*(exp(p2*t)*x1-exp(p1*t)*x2)*COS(k*PI*t/l)) t = 0 ak0 = U*(l*(exp(p1*t)*p1*x2-exp(p2*t)*p2*x1)*COS(k*PI*t/l)-k*PI*(exp(p2*t)*x1-exp(p1*t)*x2)*SIN(k*PI*t/l)) bk0 = U*(l*(exp(p1*t)*p1*x2-exp(p2*t)*p2*x1)*SIN(k*PI*t/l)+k*PI*(exp(p2*t)*x1-exp(p1*t)*x2)*COS(k*PI*t/l)) jm = Lx*(p1-p2)*x1*x2 ak(I) = (akq-ak0)/jm $bk(I) = (bkq-bk0)/jm fk(I) = k/l/2 *ENDDO