Aitken’s and Neville’s Interpolation

Numerical Methods

Source inspiration: (Mathew 2000-2019).

Animations

The animation below illustrates Neville’s algorithm for polynomial interpolation of \(f(x) = \cos(x)\) at 5 nodes \(x_0 = 0,\, \tfrac{\pi}{4},\, \tfrac{\pi}{2},\, \tfrac{3\pi}{4},\, \pi\). Each frame adds one node, building successively higher-degree interpolants \(P_{0,j}(x)\).

Julia source scripts that generated these animations are linked under each case.

Case 1 — Neville’s algorithm: degree 0 through 4 for \(\cos(x)\)

Behavior: Neville’s algorithm builds the degree-\(n\) interpolating polynomial by combining two polynomials of degree \(n-1\): \(P_{i,j}(x) = \tfrac{(x-x_i)P_{i+1,j} - (x-x_j)P_{i,j-1}}{x_j - x_i}\). Adding each node refines the polynomial. By degree 4 (using all 5 nodes) the interpolant closely matches \(\cos(x)\) on the whole interval.

Julia source

Neville algorithm building successive polynomial interpolants of cos(x) from degree 0 constant to degree 4 at 5 nodes at 0, pi/4, pi/2, 3pi/4, pi

Derivation Notes (Planned)

Short derivations will be added to explain the core equations and assumptions.

References

Mathew, John H. 2000-2019. Numerical Analysis - Numerical Methods Modules. https://web.archive.org/web/20190808102217/http://mathfaculty.fullerton.edu/mathews/n2003/NumericalUndergradMod.html.