Abstract
I consider an arbitrary-precision computation of the incomplete Gamma function from the Legendre continued fraction. Using the method of generating functions, I compute the convergence rate of the continued fraction and find a direct estimate of the necessary number of terms. This allows to compare the performance of the continued fraction and of the power series methods. As an application, I show that the incomplete Gamma function Γ (a, z) can be computed to P digits in at most O(P) long multiplications uniformly in z for Re z > 0. The error function of the real argument, erf x, requires at most O(P 2/3) long multiplications.
Access this chapter
Tax calculation will be finalised at checkout
Purchases are for personal use only
Preview
Unable to display preview. Download preview PDF.
References
M. Abramowitz and I. Stegun, eds., Handbook of special functions, National Bureau of Standards, 1964.
W. Gautschi, ACM TOMS 5, p. 466 (1979); ibid., “Algorithm 542”, p. 482.
A. R. DiDonato and A. H. Morris, Jr., ACM TOMS 12, p. 377 (1986).
W. Gautschi, ACM TOMS 25, p. 101 (1999).
D. M. Smith, “Algorithm 814”, ACM TOMS 27, p. 377 (2001).
N. M. Temme, SIAM J. Math. Anal. 10, pp. 757 (1979).
H. C. Thacher, Jr., “Algorithm 180”, Comm. ACM 6, p. 314 (1963).
D. M. Smith, Math. Comp. 52, p. 131 (1989).
B. Haible and T. Papanikolaou, LNCS 1423, p. 338 (Springer, 1998).
W. H. Press, S. A. Teukolsky, W. T. Vetterling, B. P. Flannery, Numerical recipes in C, 2nd ed., Cambridge University Press, 1992.
W. B. Jones and W. J. Thron, Math. Comp. 28, p. 795 (1974).
A. H. Karp and P. Markstein, ACM TOMS 23, p. 561 (1997).
Sh. E. Tsimring, Handbook of special functions and definite integrals: algorithms and programs for calculators, Radio and communications (publisher), Moscow, 1988 (in Russian).
J. H. McCabe, Math. Comp. 28, p. 811 (1974).
D. A. Field, Math. Comp. 31, p. 495 (1977).
F. W. J. Olver, Asymptotics and special functions, Academic Press, 1974.
Author information
Authors and Affiliations
Editor information
Editors and Affiliations
Rights and permissions
Copyright information
© 2003 Springer-Verlag Berlin Heidelberg
About this paper
Cite this paper
Winitzki, S. (2003). Computing the Incomplete Gamma Function to Arbitrary Precision. In: Kumar, V., Gavrilova, M.L., Tan, C.J.K., L’Ecuyer, P. (eds) Computational Science and Its Applications — ICCSA 2003. ICCSA 2003. Lecture Notes in Computer Science, vol 2667. Springer, Berlin, Heidelberg. https://doi.org/10.1007/3-540-44839-X_83
Download citation
DOI: https://doi.org/10.1007/3-540-44839-X_83
Published:
Publisher Name: Springer, Berlin, Heidelberg
Print ISBN: 978-3-540-40155-1
Online ISBN: 978-3-540-44839-6
eBook Packages: Springer Book Archive