High-precision evaluation of the Bessel functions via Hadamard series

  • Authors:
  • R. B. Paris

  • Affiliations:
  • Division of Complex Systems, University of Abertay Dundee, Dundee DD1 1HG, UK

  • Venue:
  • Journal of Computational and Applied Mathematics
  • Year:
  • 2009

Quantified Score

Hi-index 7.29

Visualization

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)=@c(a,z)/@C(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 @q^k, where 0