Approximation of analytic functions: a method of enhanced convergence

We deal with a method of enhanced convergence for the approximation of analytic functions. This method introduces conformal transformations in the approximation problems, in order to help extract the values of a given analytic function from its Taylor expansion around a point. An instance of this method, based on the Euler transform, has long been known; recently we introduced more general versions of it in connection with certain problems in wave scattering. In ?2 we present a general discussion of this approach. As is known in the case of the Euler transform, conformal transformations can enlarge the region of convergence of power series and can enhance substantially the convergence rates inside the circles of convergence. We show that conformal maps can also produce a rather dramatic improvement in the conditioning of Pade approximation. This improvement, which we discuss theoretically for Stieltjes-type functions, is most notorious in cases of very poorly conditioned Pade problems. In many instances, an application of enhanced convergence in conjunction with Pade approximation leads to results which are many orders of magnitude more accurate than those obtained by either classical Pade approximants or the summation of a truncated enhanced series.

[1]  C. Isenberg Expansion of the Vibrational Spectrum at Low Frequencies , 1966 .

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

[3]  W. Gautschi On Generating Orthogonal Polynomials , 1982 .

[4]  J. M. Taylor,et al.  The condition of gram matrices and related problems , 1978, Proceedings of the Royal Society of Edinburgh: Section A Mathematics.

[5]  Fernando Reitich,et al.  Numerical solution of diffraction problems: a method of variation of boundaries , 1993 .

[6]  P. Morse,et al.  Methods of theoretical physics , 1955 .

[7]  G. A. Baker,et al.  THE THEORY AND APPLICATION OF THE PADE APPROXIMANT METHOD , 1964 .

[8]  G. Forsythe,et al.  Computer solution of linear algebraic systems , 1969 .

[9]  Albert Edrei,et al.  Sur les déterminants récurrents et les singularités d'une fonction donnée par son développement de Taylor , 1939 .

[10]  S. Gustafson On stable calculation of linear functionals , 1979 .

[11]  Fernando Reitich,et al.  Solution of a boundary value problem for the Helmholtz equation via variation of the boundary into the complex domain , 1992, Proceedings of the Royal Society of Edinburgh: Section A Mathematics.

[12]  Yudell L. Luke,et al.  Algorithms for the Computation of Mathematical Functions , 1977 .

[13]  J. Uhlig C. Forsythe and C. B. Moler, Computer Solution of Linear Algebraic Systems. (Series in Automatic Computation) XI + 148 S. Englewood Cliffs, N.J. 1967. Prentice-Hall, Inc. Preis geb. 54 s. net , 1972 .

[14]  Stanley Cabay,et al.  Algebraic Computations of Scaled Padé Fractions , 1986, SIAM J. Comput..

[15]  C. Isenberg Moment Calculations in Lattice Dynamics. I. fcc Lattice with Nearest-Neighbor Interactions , 1963 .

[16]  L. Mead,et al.  Maximum entropy in the problem of moments , 1984 .

[17]  G. A. Baker,et al.  An investigation of the applicability of the Padé approximant method , 1961 .

[18]  W. Gautschi Construction of Gauss-Christoffel quadrature formulas , 1968 .

[19]  R. Scraton A note on the summation of divergent power series , 1969, Mathematical Proceedings of the Cambridge Philosophical Society.

[20]  Yudell L. Luke,et al.  Computations of coefficients in the polynomials of Pade´approximations by solving systems of linear equations , 1980 .

[21]  M. Abramowitz,et al.  Mathematical functions and their approximations , 1975 .

[22]  P. Graves-Morris The numerical calculation of Padé approximants , 1979 .

[23]  G. Forsythe,et al.  Computer solution of linear algebraic systems , 1969 .