A method for computing Bessel function integrals (Q1098552)

An algorithm is given for calculating integrals containing Bessel functions of integer or integer plus one-half order. The calculation in the Bessel functions involves a one-dimensional Fourier sine or cosine transform followed by evaluation of the coefficients of the Chebyshev series of the Fourier-transformed function. In the case of spherical Bessel functions, the calculation consists of the Fourier sine or cosine transform and evaluation of the Legendre expansion coefficients.