Computation of Jn(x) by Numerical Integration
It is shown to be practical to compute Jn(x) by
numerical integration of its integral representation 
using the trapezoidal rule. The error in this
approximation was studied empirically.
CACM April, 1969
Stroud, A. H.
Kohli, J. P.
