Interpolació inversa

Suposem coneguts (xk,fk) les dades corresponents a una funció f(x) i volem trobar una aproximació del valor de x tal que f(x)=c, on c és un valor donat.

El que farem és resoldre l'equació x=g(c) on g és la funció inversa de f. Llavors, interpolarem aquesta funció g(y) i l'avaluarem a y=c, és a dir, si seguim el mètode de Newton posarem a la primera fila els valors fj i en la segona els valors xj i procedirem de la mateixa manera.

Exemple

Per exemple, suposem que volem calcular un zero de la funció f(x)=x315x+4 sabent que aquest està a prop de x=0.3. Llavors farem interpolació quadràtica, per exemple, de la inversa de f(x). Primer, doncs, avaluem la funció en tres punts prop de x=0.3:

x 0.2 0.3 0.4
f(x) 1.008 0.473 1.936

Ara escrivim la taula per calcular les diferències dividides de Newton, però intercanviant les columnes, obtenint els coeficients del polinomi interpolador:

1.008 0.2    
    0.0675  
0.473 0.3   0.00028963
    0.0684  
1.936 0.4    

D'aquesta forma el polinomi interpolació és:

P3(y)=0.2+0.0675(y1.008)+0.00028963(y1.008)(y0.473)=0.26790190900.067654y+0.00028963y2

Així una aproximació del zero de la funció és:

P3(0)=0.2679019090