title  "BESSEL FUNCTIONS"


header <napatool.hdr>

fs  1.0e6

directive TIME_OFFSET  1.0e-6   //  avoid TIME = 0!

dvar  freq   1000.0
dvar  omega  _2pi_ * freq

node j0 dalgebra  bessel_j(0, omega * TIME)
node j1 dalgebra  bessel_j(1, omega * TIME)
node j2 dalgebra  bessel_j(2, omega * TIME)

node y0 dalgebra  bessel_y(0, omega * TIME)
node y1 dalgebra  bessel_y(1, omega * TIME)
node y2 dalgebra  bessel_y(2, omega * TIME)

node i0 dalgebra  bessel_i(0, omega * TIME)
node i1 dalgebra  bessel_i(1, omega * TIME)
node i2 dalgebra  bessel_i(2, omega * TIME)

node k0 dalgebra  bessel_k(0, omega * TIME)
node k1 dalgebra  bessel_k(1, omega * TIME)
node k2 dalgebra  bessel_k(2, omega * TIME)

output "stdout"  j0 j1 j2 ...   // bessel 1st kind
                 y0 y1 y2 ...   // bessel 2nd kind
                 i0 i1 i2 ...   // modified bessel 1st kind
                 k0 k1 k2       // modified bessel 2nd kind

terminate TIME >  1.0e-3