Lagrange and Newtons Interpolation
Lagrange's and Newton's method of polynomial interpolation.
newtonInterp(x, y, xs = c()) lagrangeInterp(x, y, xs)
x, y |
x-, y-coordinates of data points defining the polynomial. |
xs |
either empty, or a vector of points to be interpolated. |
Straightforward implementation of Lagrange's Newton's method
(vectorized in xs).
A vector of values at xs of the polynomial defined by x,y.
Each textbook on numerical analysis.
p <- Poly(c(1, 2, 3)) fp <- function(x) polyval(p, x) x <- 0:4; y <- fp(x) xx <- linspace(0, 4, 51) yy <- lagrangeInterp(x, y, xx) yy <- newtonInterp(x, y, xx) ## Not run: ezplot(fp, 0, 4) points(xx, yy) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.