You are here: MIMS > EPrints
MIMS EPrints

2006.139: The numerical stability of barycentric Lagrange interpolation

2006.139: Nicholas J Higham (2004) The numerical stability of barycentric Lagrange interpolation. IMA Journal of Numerical Analysis, 24. pp. 547-556. ISSN 1464-3642

Full text available as:

PDF (Author PDF of final accepted paper) - Archive staff only - Requires a PDF viewer such as GSview, Xpdf or Adobe Acrobat Reader
113 Kb

DOI: 10.1093/imanum/24.4.547


The Lagrange representation of the interpolating polynomial can be rewritten in two more computationally attractive forms: a modified Lagrange form and a barycentric form. We give an error analysis of the evaluation of the interpolating polynomial using these two forms. The modified Lagrange formula is shown to be backward stable. The barycentric formula has a less favourable error analysis, but is forward stable for any set of interpolating points with a small Lebesgue constant. Therefore the barycentric formula can be significantly less accurate than the modified Lagrange formula only for a poor choice of interpolating points. This analysis provides further weight to the argument of Berrut and Trefethen that barycentric Lagrange interpolation should be the polynomial interpolation method of choice.

Item Type:Article
Uncontrolled Keywords:Keywords: polynomial interpolation; Lagrange interpolation; barycentric formula; rounding error analysis; backward error; forward error; Lebesgue constant.
Subjects:MSC 2000 > 15 Linear and multilinear algebra; matrix theory
MSC 2000 > 65 Numerical analysis
MIMS number:2006.139
Deposited By:Miss Louise Stait
Deposited On:27 June 2006

Download Statistics: last 4 weeks
Repository Staff Only: edit this item