cylBessel_jn

Bessel function of integer order

Returns Bessel function of order n, where n is a (possibly negative) integer.

The ratio of jn(x) to j0(x) is computed by backward recurrence. First the ratio jn/jn-1 is found by a continued fraction expansion. Then the recurrence relating successive orders is applied until j0 or j1 is reached.

If n = 0 or 1 the routine for j0 or j1 is called directly.

real
cylBessel_jn
(
int n
,
real x
)

Bugs

Not suitable for large n or x.

Meta