On the Evaluation of Bessel Functions

In this: publication
By this: publisher
In this Subject: <a title="Mathematics and Statistics" href="/content/subcat?j_subject=195&j_availability=">Mathematics and Statistics
By this author: Matviyenko G.

Author: Matviyenko G.

Source: Applied and Computational Harmonic Analysis, Volume 1, Number 1, December 1993 , pp. 116-135(20)

Buy & download fulltext article:

Price: $52.63 plus tax (Refund Policy)

Abstract:

In the present paper we describe an algorithm for the evaluation of Bessel functions J(x), Y(x) and H^(j)(x) (j = 1, 2) of arbitrary positive orders and arguments at a constant CPU time. The algorithm employs Taylor series, the Debye asymptotic expansions, and numerical evaluation of the Sommerfeld integral, and is based on the following two observations. (1) The Debye asymptotic expansions, contrary to what appears to be a popular belief, are not expansions in inverse powers of (large) parameter but turn out to be uniform expansions in inverse powers of (large) parameter g₁ = (x - )x^1/3 for x > and (large) parameter g₂ ( - x)/^1/3 for x < . (2) For x and such that both Taylor and Debye expansions do not provide a specified accuracy Bessel functions can be computed at a constant CPU time via (numerical) evaluation of the Sommerfeld integral along contours of steepest descents. In addition, in Appendix B we obtain certain new estimates concerning decay of the functions J(x) and -1/Y(x) of fixed x and large , and in Appendix C we show that functions J(x) of integer provide the solution for a certain system of coupled harmonic oscillators.Copyright 1993, 1999 Academic Press