Abstract
We introduce an algorithm for the evaluation of the Incomplete Gamma Function, P(m, x), for all m, x > 0. For small m, a classical recursive scheme is used to evaluate P(m, x), whereas for large m a newly derived asymptotic expansion is used. The number of operations required for evaluation is O(1) for all x and m. Nearly full double and extended precision accuracies are achieved in their respective environments. The performance of the scheme is illustrated via several numerical examples.
Similar content being viewed by others
References
Abramowitz, M., Stegun, I.A. (eds.): Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables. US Government Printing Office, Washington (1964)
DiDonato, A., Morris, A.H.: Computation of the Incomplete Gamma Function Ratios and their Inverse. ACM TOMS 12.4, 377–393 (1986)
Feller, W.: An Introduction to Probability and Its Applications, 3rd Ed, vol. 1. Wiley, New York (1968)
Gautschi, W.: A computational procedure for incomplete gamma functions. ACM TOMS 5.4, 466–481 (1979)
Gil, A., Seura, J., Temme, N.M.: Efficient and accurate algorithms for the computation and inversion of the incomplete. SIAM J Sci. Comput. 34.6, A2965–A2981 (2013)
Gradshteyn, I.S., Ryzhik, I.M. In: Jeffrey, A., Zwillinger, D. (eds.) : Table of Integrals, Series, and Products. Academic Press, San Diego (2000)
Nemes, G.: exponential improvement for Hermite’s asymptotic expansion for the gamma function. Appl. Anal Error bounds Discrete Math. 7.1, 161–179 (2013)
Nemes, G.: An explicit formula for the coefficients in Laplace’s method. Constr. Approx. 38.3, 471–487 (2013)
Temme, N.M.: The asymptotic expansion of the incomplete gamma functions. SIAM J. Math. Anal. 10.4, 757–766 (1979)
Temme, N.M.: On the Computation of the Incomplete Gamma Functions for Large Values of the Parameters. Algorithms for Approximation, pp. 479–489. Clarendon Press, New York (1987)
Author information
Authors and Affiliations
Corresponding author
Additional information
Communicated by: Zydrunas Gimbutas
Rights and permissions
About this article
Cite this article
Greengard, P., Rokhlin, V. An algorithm for the evaluation of the incomplete gamma function. Adv Comput Math 45, 23–49 (2019). https://doi.org/10.1007/s10444-018-9604-x
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s10444-018-9604-x