Computation of bessel functions in light scattering studies

Computations of light scattering require finding Bessel functions of a series of orders. These are found most easily by recurrence, but excessive rounding errors may accumulate. Satisfactory procedures for cylinder and sphere functions are described. If argument z is real, find Y(n)(z) by recurrence to high orders. From two high orders of Y(n)(z) estimate J(n)(z). Use backward recurrence to maximum J(n)(z). Correct by forward recurrence to maximum. If z is complex, estimate high orders of J(n)(z) without Y(n)(z) and use backward recurrence.