De. Amos, A REMARK ON ALGORITHM-644 - A PORTABLE PACKAGE FOR BESSEL-FUNCTIONS OF A COMPLEX ARGUMENT AND NONNEGATIVE ORDER, ACM transactions on mathematical software, 21(4), 1995, pp. 388-393
Algorithm 644 computes all major Bessel functions of a complex argumen
t and of nonnegative order. Single-precision routine CBESY and double-
precision routine ZBESY are modified to reduce the computation time fo
r the Y Bessel function by approximately 25% over a wide range of argu
ments and orders. Quick check (driver) programs that exercise the pack
age are also modified to make tests more meaningful over larger region
s of the complex plane.