THE NUMERICAL COMPUTATION OF BESSEL-FUNCTIONS OF THE 1ST KIND AND 2ND KIND FOR INTEGER ORDERS AND COMPLEX ARGUMENTS

A rigorous algorithm for the computation of integer order Bessel functions of the first and second kind with complex arguments is discussed. The algorithm makes use of backward recurrence for the computation of Bessel functions of the first kind where applicable, and of Hankel's asymptotic expansion for large arguments. The argument is limited only by the ability of the computer to handle large numbers. Orders up to a few thousand have been computed to 13-digit accuracy. © 1990 IEEE