Computation of special functions by Padé approximants with orthogonal polynomial denominators

This paper deals with the computation of special functions of mathematics and physics in the complex domain using continued fraction (one-point or two-point Padé) approximants. We consider three families of continued fractions (Stieltjes fractions, real J-fractions and non-negative T-fractions) whose denominators are orthogonal polynomials or Laurent polynomials. Orthogonality of these denominators plays an important role in the analysis of errors due to numerical roundoff and truncation of infinite sequences of approximants. From the rigorous error bounds described one can determine the exact number of significant decimal digits contained in the approximation of a given function value. Results from computational experiments are given to illustrate the methods.

[1]  D. Owen Handbook of Mathematical Functions with Formulas , 1965 .

[2]  P. Henrici Applied and Computational Complex Analysis, Volume 2: Special Functions, Integral Transforms, Asymptotics, Continued Fractions , 1977 .

[3]  Henry C. Thacher,et al.  Applied and Computational Complex Analysis. , 1988 .

[4]  G. A. Baker Essentials of Padé approximants , 1975 .

[5]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[6]  Annie Cuyt,et al.  Rounding error analysis for forward continued fraction algorithms , 1985 .

[7]  William B. Jones,et al.  A survey of truncation error analysis for Padé and continued fraction approximants , 1993 .

[8]  William B. Jones,et al.  Numerical stability in evaluating continued fractions , 1974 .

[9]  C. Brezinski Procedures for estimating the error in Padé approximation , 1989 .

[10]  Haakon Waadeland,et al.  Continued fractions with applications , 1994 .

[11]  K. S. Kölbig,et al.  Errata: Milton Abramowitz and Irene A. Stegun, editors, Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables, National Bureau of Standards, Applied Mathematics Series, No. 55, U.S. Government Printing Office, Washington, D.C., 1994, and all known reprints , 1972 .

[12]  William B. Jones,et al.  Continued fractions in numerical analysis , 1988 .

[13]  W. J. Thron,et al.  Continued Fractions: Analytic Theory and Applications , 1984 .

[14]  Peter Henrici,et al.  Truncation error estimates for Stieltjes fractions , 1966 .

[15]  C. Brezinski Error Estimates in Padé Approximation , 1999 .

[16]  Walter Gautschi,et al.  Efficient computation of the complex error function , 1970 .

[17]  L. R. Shenton,et al.  Continued fractions in statistical applications , 1991 .