Numerical computation of Tricomi's psi function by the trapezoidal rule

Computing - Tập 39 - Trang 271-279 - 1987
G. Allasia1, R. Besenghi1
1Dipartimento di Matematica, Università di Torino, Torino, Italy

Tóm tắt

The trapezoidal rule is applied to the numerical calculation of an integral representation of Tricomi's psi function Ψ(a, c; x) fora, x ε ℝ+ andc ε ℝ. The unexpectedly high accuracy is explained by means of a careful investigation of the remainder terms of the Euler-Maclaurin formula. A simple and efficient numerical procedure for obtaining values of the psi function is given.

Tài liệu tham khảo

Tricomi, F. G.: Funzioni Ipergeometriche Confluenti. Roma: Cremonese 1954. Slater, L. J.: Confluent Hypergeometric Functions. Cambridge: Cambridge Univ. Press 1960. Luke, Y. L.: The Special Functions and their Approximations, Vols. I, II. New York: Academic Press 1969. Davis, P. J., Rabinowitz, P.: Methods of Numerical Integration, 2nd ed. New York: Academic Press 1984. Allasia, G., Besenghi, R.: Numerical calculation of incomplete gamma functions by the trapezoidal rule. Num. Math.50, 419–428 (1987). Allasia, G., Besenghi, R.: Sul calcolo numerico delle funzioni gamma e digamma mediante la formula del trapezio. To appear on Boll Unione Mat. Italiana. Abramowitz, M., Stegun, I. A.: Handbook of Mathematical Functions. New York: Dover Publications 1970. Wimp, J.: On the computation of Tricomi's Ψ function. Computing13, 195–203 (1974). Temme, N. M.: The numerical computation of the confluent hypergeometric functionU(a, b; z). Num. Math.41, 63–82 (1983). Luke, Y. L.: Mathematical Functions and their Approximations. New York: Academic Press 1975.