Neville's Method
Neville's's method of polynomial interpolation.
neville(x, y, xs)
x, y |
x-, y-coordinates of data points defining the polynomial. |
xs |
single point to be interpolated. |
Straightforward implementation of Neville's method; not yet vectorized.
Interpolated value 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 <- numeric(51) for (i in 1:51) yy[i] <- neville(x, y, xx[i]) ## Not run: ezplot(fp, 0, 4) points(xx, yy) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.