Simple Search

Query:

Enter a query and click »Search«...

Format:

Display: entries per page entries

Zbl 0618.65009
Allasia, G.; Besenghi, R.
Numerical computation of Tricomi's psi function by the trapezoidal rule. (English)
[J] Computing 39, 271-279 (1987). ISSN 0010-485X; ISSN 1436-5057/e

The trapezoidal rule is applied to the numerical calculation of the integral representation of Tricomi's psi function $\Psi (a,c;x)=2x\sp{1- c}/\Gamma (a)\int \sp{+\infty}\sb{0}e\sp{-u\sp 2}(u\sp{2a-1}/(x+u\sp 2)\sp{1-c+a})du$ for $a,x\in {\bbfR}\sp +$, $c\in {\bbfR}$. The unexpectedly high accuracy is explained by means of a careful investigation in the complex field of the Euler-Maclaurin formula, and particularly of its remainder terms, considered as an extension of the trapezoidal rule. Since the same method has been used previously to evaluate the complementary incomplete gamma function [the authors, Numer. Math. 50, 419-428 (1987; Zbl 0593.65017)], the Euler gamma function and the digamma function, the present paper limits itself to quote the main features of the method and describes fully the peculiarities of this application. \par A simple and efficient numerical procedure for obtaining values of the psi function is given; moreover, to reduce the amount of calculation, an iterative algorithm for the evaluation of the trapezoidal rule, similar to Horner's scheme for polynomials, is suggested.

MSC 2000:: *65D20 Computation of special functions
33B15 Gamma-functions, etc.
65B15 Euler-Maclaurin formula

Keywords: trapezoidal rule; integral representation of Tricomi's psi function; high accuracy; Euler-Maclaurin formula; complementary incomplete gamma function; Euler gamma function; digamma function; iterative algorithm

Citations: Zbl 0593.65017

Comment on this Item PDF MathML XML ASCII DVI PS BibTeX Online Ordering Article Journal

	Scientific prize winners of the ICM 2010
	Overhang
	Lie groups, physics and geometry. An introduction for physicists, engineers and chemists.

Simple Search

Zentralblatt MATH Berlin [Germany]