Deducción: coeficientes de Newton por diferencias divididas

De dónde salen b0b_0 y b1b_1 al forzar que el polinomio pase por los puntos, y cómo la recursión genera todos los coeficientes de orden superior.

Forzar que pase por dos puntos

  1. Partimos de la forma lineal p1(x)=b0+b1(xx0)p_1(x)=b_0+b_1(x-x_0) e imponemos que pase por (x0,f(x0))(x_0,f(x_0)) y (x1,f(x1))(x_1,f(x_1)):

    {p1(x0)=f(x0)=b0p1(x1)=f(x1)=b0+b1(x1x0)\begin{cases} p_1(x_0)=f(x_0)=b_0\\ p_1(x_1)=f(x_1)=b_0+b_1(x_1-x_0) \end{cases}
  2. De la primera ecuación b0=f(x0)b_0=f(x_0). Sustituyendo en la segunda y despejando b1b_1:

    b1=f(x1)f(x0)x1x0=f[x1,x0]b_1=\frac{f(x_1)-f(x_0)}{x_1-x_0}=f[x_1,x_0]

A b1b_1 se le llama diferencia dividida de primer orden. Repetir el mismo argumento con tres puntos da la diferencia de segundo orden, y así se llega a la recursión general.

f[xn,,x0]=f[xn,,x1]f[xn1,,x0]xnx0f[x_n,\dots,x_0]=\frac{f[x_n,\dots,x_1]-f[x_{n-1},\dots,x_0]}{x_n-x_0}