Cabay, Stan; Meleshko, Ron A weakly stable algorithm for Padé approximants and the inversion of Hankel matrices. (English) Zbl 0779.65009 SIAM J. Matrix Anal. Appl. 14, No. 3, 735-765 (1993). Authors’ summary: A new algorithm, algorithm NPADE, is presented for numerically computing Padé approximants in a weakly stable fashion. By this, it is meant that if the problem is well conditioned, then algorithm NPADE produces a good solution. No restrictions are imposed on the problem being solved. Except in certain pathological cases, the cost of the algorithm is \(O(n^ 2)\), where \(n\) is the maximum degree of the polynomials comprising the Padé approximant.The operation of algorithm NPADE is controlled by a single parameter. Bounds are obtained for the computed solution and it is seen that they are a function of this parameter. Experimental results show that the bounds, while crude, reflect the actual behavior of the error. In addition, it is shown how better bounds can easily be obtained a posteriori. As another application of algorithm NPADE, it is shown that it can be used to compute stably, in a weak sense, the inverse of a Hankel or Toeplitz matrix. Reviewer: D.Braess (Bochum) Cited in 2 ReviewsCited in 26 Documents MSC: 65D15 Algorithms for approximation of functions 65F05 Direct numerical methods for linear systems and matrix inversion 41A21 Padé approximation Keywords:inversion of Hankel matrices; inversion of Toeplitz matrices; algorithm NPADE; Padé approximants; Experimental results PDFBibTeX XMLCite \textit{S. Cabay} and \textit{R. Meleshko}, SIAM J. Matrix Anal. Appl. 14, No. 3, 735--765 (1993; Zbl 0779.65009) Full Text: DOI