Hermite interpolation

Interpolation that imposes value and derivative at each node: the polynomial H2n+1H_{2n+1}, its construction from the Lagrange bases, the error, the practical route via divided differences with repeated nodes and a Bessel-function example.

What Hermite adds

So far we forced the polynomial through the points. If we also know ff's derivative at those points, we can impose it too: this controls the slope and often improves the approximation.

Hi(x)=[12(xxi)Li(xi)]Li2(x),H^i(x)=(xxi)Li2(x)H_{i}(x)=\bigl[1-2(x-x_i)L_i'(x_i)\bigr]L_i^2(x),\qquad \hat H_{i}(x)=(x-x_i)L_i^2(x)
Base blocks built from the Lagrange functions and their derivatives.

Procedure and error

Five steps
  1. Compute the Lagrange functions Li(x)L_i(x).

  2. Differentiate them to obtain Li(x)L_i'(x).

  3. Form the blocks Hi(x)H_i(x).

  4. Form the blocks H^i(x)\hat H_i(x).

  5. Sum everything to obtain H2n+1(x)H_{2n+1}(x).