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