¿Cómo programar el método de interpolación de Newton en lenguaje C?

#include lt;iostream.hgt;

#include lt;math.hgt;

void main()

{

char L;

hacer

{

doble M[100][100]; [100], y[100];

doble X=1, xx=0, w=1, N=0, P, R=1;

int n;

coutlt;lt;"Ingrese el orden de diferencia de medias requerido:";

cingt;gt;n;

for(int i=0;ilt ; = n; i )

{

coutlt;lt;"Por favor, introduzca el valor de x"lt;lt;ilt;lt;": "lt;lt;endl;

cingt;

cingt; y[i];

M[i][0]=x[i];

M[i][1] =y[i];

}

for( int j=2; jlt; =n 1; j )

{

for( i=1;ilt;=n;i)

{

M[i][j]=(M[ i][j-1]-M[ i-1][j-1])/(M[i][0]-M[i-j 1][0]);

}

}

for(i=1;ilt;=n;i)

{

coutlt;lt;"its"lt; lt;ilt;lt;"orden de diferencia de medias Para: "lt;lt;M[i][i 1]lt;lt;endl;

}

coutlt;lt; " Introduzca el valor de x: x=" ;

cingt;gt;xx;

for(i=0;ilt;n;i)

{

X* =xx-x[i];

N =M[i 1][i 2]*X;

P=M [0][1] N;

}

coutlt;lt;"su valor de función: y="lt;lt;Plt;lt;endl;

for(i=0;ilt;n ;i )

{

w*=xx-x[i];

R=fabs( M[n][n 1]*w) ;

}

coutlt;lt;"El error de truncamiento: R="lt;lt;Rlt;lt;endl;

coutlt;lt;endllt ;lt;"¿Quieres calcular otras interpolaciones? Presione 'y' en caso afirmativo; de lo contrario, presione 'n'"lt;lt;endl;

cingt;gt;L;

p> }mientras(L=='y');

}