Abstract
The numerical evaluation of Bessel function integrals may be difficult when the Bessel function is rapidly oscillating in the interval of integration. In the method presented here, the smooth factor of the integrand is replaced by a truncated Chebyshev series approximation and the resulting integral is computed exactly. The numerical aspects of this exact integration are discussed.
Similar content being viewed by others
References
M. Abramowitz and I. Stegun (Eds.),Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables, Dover Publ., New York (1965), p. 364.
W. L. Anderson,Fast Hankel transforms using related and lagged convolutions, TOMS 8 (1982), 344–368.
M. Branders,The asymptotic behaviour of solutions of difference equations, Bull. Soc. Math. Belg. 26 (1974), 255–260.
M. Branders and R. Piessens,Algorithm 001, An extension of Clenshaw-Curtis quadrature, J. CAM 1 (1975), 55–65.
C. W. Clenshaw,Curve fitting with a digital computer, Computer J. 2 (1960), 170–173.
C. W. Clenshaw and A. R. Curtis,A method for numerical integration on an automatic computer, Num. Math. 2 (1960), 197–205.
P. Cornille,Computation of Hankel transforms, SIAM Rev. 14 (1972), 278–285.
B. Davies,Integral Transforms and their Applications, Springer-Verlag, Heidelberg (1978), pp. 237–247.
P. J. Davis and P. Rabinowitz,Methods of Numerical Integration, Academic Press, New York (1975), pp. 118–132 and pp. 178–180.
J. Denef and R. Piessens,The asymptotic behaviour of solutions of difference equations of Poincaré:s type, Bull. Soc. Math. Belg. 26 (1974), 133–146.
A. Erdélyi,Asymptotic representations of Fourier integrals and the method of stationary phase, J. Soc. Industr. Appl. Math. 3 (1955), 17–27.
A. Erdélyi, W. Magnus, F. Oberhettinger and F. G. Tricomi,Tables of Integral Transforms, vol. 2, McGraw-Hill, New York (1954), p. 5.
W. Gautschi,Algorithm 236: Bessel functions of the first kind, Comm. ACM 7 (1964), 479–480.
W. Gautschi,Computational aspects of three-term recurrence relations, SIAM Rev. 9 (1967), 24–82.
W. Gautschi,On the construction of Gaussian rules from modified moments, Math. Comp. 24 (1970), 245–260.
W. M. Gentleman,Algorithm 424, Clenshaw-Curtis quadrature, Comm. ACM 15 (1972), 353–355.
P. Linz,A method for computing Bessel function integrals, Math. Comp. 26 (1972), 504–513.
I. M. Longman,A short table of ∫ ∞x J 0(t)t −n dt and ∫ ∞ Jx 1(t)t−n dt, MTAC 13 (1959), 306–311.
I. M. Longman,A method for the numerical evaluation of finite integrals of oscillatory functions, Math. Comp. 14 (1960), 53–59.
Y. L. Luke,The special functions and their approximations, vol. 2, Academic Press, New York and London (1969), p. 53 and p. 334.
J. Oliver,The numerical solution of linear recurrence relations, Num. Math. 11 (1968), 349–360.
R. Piessens,Gaussian quadrature formulas for integrals involving Bessel functions, Microfiche section of Math. Comp. 26 (1972).
R. Piessens and M. Branders,Approximations for Bessel functions and their application in the computation of Hankel transforms, Comp. & Maths. with Appls. 8 (1982), 305–311.
R. Piessens, E. de Doncker, C. Überhuber and D. Kahaner,QUADPACK: A Subroutine Package for Numerical Integration, to be published by Springer-Verlag, Heidelberg, pp. 66–67.
M. A. Snyder,Chebyshev Methods in Numerical Approximation, Prentice-Hall, Englewood Cliffs, N.J. (1966), p. 103.
G. N. Watson,Treatise on the Theory of Bessel Functions, Cambridge University Press, 2nd ed. (1958), pp. 537–538.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Piessens, R., Branders, M. Modified clenshaw-curtis method for the computation of Bessel function integrals. BIT 23, 370–381 (1983). https://doi.org/10.1007/BF01934465
Received:
Revised:
Issue Date:
DOI: https://doi.org/10.1007/BF01934465