拉格朗日插值多项式(3)

拉格朗日插值多项式

printf("Your data:\n"); for(i=0;i<M;i++)

{scanf("%lf",&z[i].x); scanf("%lf",&z[i].y); }

printf("Input x:\n"); scanf("%lf",&x); for(k=0;k<M;k++) {p=1;

for(j=0;j<M;j++) if(j!=k)

p=p*(x-z[j].x)/(z[k].x-z[j].x); sum=sum+p*z[k].y; }

printf("The result is %lf\n",sum); getch(); }