  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.
