Numerical Calculation of Bessel Functions
  A new computational procedure is offered to provide simple, accurate and flexible methods for using modern computers to give numerical evaluations of the various Bessel functions. The Trapezoidal Rule, applied to suitable integral representations, may become the method of choice for evaluation of the many Special Functions of mathematical physics. 
Charles Schwartz
math.NA physics.comp-ph