#include
#include
#include

/* Interpolation with Lagrange method */

main()

{

double x[10],y[10];
double xx,yy;
double pn,li;
double deg=3.141592/180.;

int i,j;
int n,m;

n=10;
for (i=0;i x[i]=(double)i *10.0*deg; y[i]=sin(x[i]);
}

xx=45.0*deg;
yy=sin(xx);

pn=0.0;
for (i=0;i li=1;
for (j=0;j if(i != j){
li *= (xx-x[j])/(x[i]-x[j]);
}
}
pn+=y[i]*li;
}

printf(" 45 deg:%10.5f rad Pn:%7.5f True:%7.5f \n",xx,pn,yy);

return(0);
}