Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
94 5. INTERPOLACI ÓN<br />
Si se eligen 10 puntos equiespaciados se obtiene un polinomio como muestra la Figura 5.2.<br />
Si consideramos el error numérico, que es el que se obtiene como diferencia máxima entre f<br />
y el polinomio evaluados en una malla suficientemente fina (puntos equiespaciados con paso<br />
h = 0,01) se tiene un error de 0.4303. . . Tan sólo al considerar 25 puntos equiespaciados (tomados<br />
a intervalos de longitud 0.25) se obtiene un error numérico menor que 10 −6 . En este caso, en<br />
una figura como la anterior, los gráficos del polinomio y la función se confunden.<br />
4<br />
3<br />
2<br />
1<br />
0<br />
−1<br />
−2<br />
−3<br />
−4<br />
−5<br />
−6<br />
−1 −0.8 −0.6 −0.4 −0.2 0 0.2 0.4 0.6 0.8 1<br />
Figura 5.2. Interpolación de f(x) = cos(x) 3 en [−3, 3], (10 puntos equiespaciados)<br />
3. Forma de Newton<br />
La forma de Newton es conveniente para calcular el polinomio, en Pn, que interpola a una<br />
función f en x0, . . . , xn−1, xn una vez conocido el polinomio interpolador de f en x0, . . . , xn−1.<br />
La forma de Newton del polinomio interpolador puede verse como una generalización del polinomio<br />
de Taylor asociado a una función. En esta construcción aparecen las diferencias divididas<br />
que presentamos a continuación.<br />
Primera diferencia dividida<br />
Segunda diferencia dividida<br />
f[x0, x1] = f(x1) − f(x0)<br />
.<br />
x1 − x0<br />
f[x0, x1, x2] = f[x1, x2] − f[x0, x1]<br />
.<br />
x2 − x0<br />
Así sucesivamente se define la diferencia de orden k asociada a los puntos x0, . . . , xk,