“…Moreover, Higham has shown in [22] that the barycentric Lagrange formula (3.2) is stable for any set of interpolating points with a small Lebesgue constant, such as Chebyshev and Gauss-Legendre points, and ought to be the standard method of polynomial interpolation [4]. For an exhaustive discussion of the theoretical advantage of the barycentric Lagrange formula, we refer the reader to [4,10,16,22]. If f (x) is continuous and of bounded variation on [−1, 1], then the Lagrange interpolation polynomial q n (x) which interpolates f (x) at the Gauss-Legendre points converges uniformly to f (x) on [−1, 1] [38].…”