Numerical computation of Tricomi's psi function by the trapezoidal rule
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.