High-precision evaluation of the Bessel functions via Hadamard series

Richard B. Paris

Research output: Contribution to journalArticle

  • 3 Citations

Abstract

We present a method of high-precision calculation of the Bessel functions using Hadamard series. Such series are absolutely convergent expansions involving the normalised incomplete gamma function P(a, z) = γ (a, z)/Γ(a) and possess early terms that behave like those in an asymptotic expansion. In the case of real variables the function P(a,z) acts as a smoothing factor on the terms of the series. We show how these series representing the Bessel functions of complex argument can be chosen so as to produce rapidly convergent series that possess terms decaying at the geometric rate ϑk , where 0 < ϑ < 1 and k is the ordinal number of the series. We give numerical examples with ϑ = 1/2 , 1/3 and 1/4 . The theory is extended to cover the confluent hypergeometric functions 1 F 1 (a; b; z) and U(a, b, z), thereby dealing with many of the special functions arising in mathematical physics.
Original languageEnglish
Pages (from-to)84-100
Number of pages17
JournalJournal of Computational and Applied Mathematics
Volume224
Issue number1
DOIs
StatePublished - 1 Feb 2009

Fingerprint

Series
Term
Bessel functions
Real variables
Physics
Incomplete gamma function
Confluent hypergeometric function
Special functions
Asymptotic expansion
Smoothing
Cover
Numerical examples
Evaluation

Cite this

Paris, Richard B. / High-precision evaluation of the Bessel functions via Hadamard series.

In: Journal of Computational and Applied Mathematics, Vol. 224, No. 1, 01.02.2009, p. 84-100.

Research output: Contribution to journalArticle

@article{45927fcf7e6f45c08f2f9ea1640a9cbb,
title = "High-precision evaluation of the Bessel functions via Hadamard series",
abstract = "We present a method of high-precision calculation of the Bessel functions using Hadamard series. Such series are absolutely convergent expansions involving the normalised incomplete gamma function P(a, z) = γ (a, z)/Γ(a) and possess early terms that behave like those in an asymptotic expansion. In the case of real variables the function P(a,z) acts as a smoothing factor on the terms of the series. We show how these series representing the Bessel functions of complex argument can be chosen so as to produce rapidly convergent series that possess terms decaying at the geometric rate ϑk , where 0 < ϑ < 1 and k is the ordinal number of the series. We give numerical examples with ϑ = 1/2 , 1/3 and 1/4 . The theory is extended to cover the confluent hypergeometric functions 1 F 1 (a; b; z) and U(a, b, z), thereby dealing with many of the special functions arising in mathematical physics.",
author = "Paris, {Richard B.}",
year = "2009",
month = "2",
doi = "10.1016/j.cam.2008.04.025",
volume = "224",
pages = "84--100",
journal = "Journal of Computational and Applied Mathematics",
issn = "0377-0427",
publisher = "Elsevier",
number = "1",

}

High-precision evaluation of the Bessel functions via Hadamard series. / Paris, Richard B.

In: Journal of Computational and Applied Mathematics, Vol. 224, No. 1, 01.02.2009, p. 84-100.

Research output: Contribution to journalArticle

TY - JOUR

T1 - High-precision evaluation of the Bessel functions via Hadamard series

AU - Paris,Richard B.

PY - 2009/2/1

Y1 - 2009/2/1

N2 - We present a method of high-precision calculation of the Bessel functions using Hadamard series. Such series are absolutely convergent expansions involving the normalised incomplete gamma function P(a, z) = γ (a, z)/Γ(a) and possess early terms that behave like those in an asymptotic expansion. In the case of real variables the function P(a,z) acts as a smoothing factor on the terms of the series. We show how these series representing the Bessel functions of complex argument can be chosen so as to produce rapidly convergent series that possess terms decaying at the geometric rate ϑk , where 0 < ϑ < 1 and k is the ordinal number of the series. We give numerical examples with ϑ = 1/2 , 1/3 and 1/4 . The theory is extended to cover the confluent hypergeometric functions 1 F 1 (a; b; z) and U(a, b, z), thereby dealing with many of the special functions arising in mathematical physics.

AB - We present a method of high-precision calculation of the Bessel functions using Hadamard series. Such series are absolutely convergent expansions involving the normalised incomplete gamma function P(a, z) = γ (a, z)/Γ(a) and possess early terms that behave like those in an asymptotic expansion. In the case of real variables the function P(a,z) acts as a smoothing factor on the terms of the series. We show how these series representing the Bessel functions of complex argument can be chosen so as to produce rapidly convergent series that possess terms decaying at the geometric rate ϑk , where 0 < ϑ < 1 and k is the ordinal number of the series. We give numerical examples with ϑ = 1/2 , 1/3 and 1/4 . The theory is extended to cover the confluent hypergeometric functions 1 F 1 (a; b; z) and U(a, b, z), thereby dealing with many of the special functions arising in mathematical physics.

U2 - 10.1016/j.cam.2008.04.025

DO - 10.1016/j.cam.2008.04.025

M3 - Article

VL - 224

SP - 84

EP - 100

JO - Journal of Computational and Applied Mathematics

T2 - Journal of Computational and Applied Mathematics

JF - Journal of Computational and Applied Mathematics

SN - 0377-0427

IS - 1

ER -