Bessel functions J_n(z) and Y_n(z) of integer order and complex argument - du Toit, C. F. (Comput. Phys. Comm., 1993)
Double-precision Pascal, maximum accuracy 7D.