牛顿插值法C语言程序

#includestdio.h#includemath.h#defineN6floatsub(floata[],floatb[],floatx,floate);voidmain(void){floatu[N]={100,121,144,169,196,225};floatv[N]={10,11,12,13,14,15};floatx,y,e,*p1,*p2;printf(InputnumberxE=:);scanf(%f%e,&x,&e);p1=u;p2=v;y=sub(p1,p2,x,e);printf(y=%f\n,y);}floatsub(float*pp1,float*pp2,floatx,floate){floata[N],b[N],t[N],y,y1,c;inti,k;for(i=0;iN;i++,pp1++){a[i]=*pp1;printf(%12.6f,a[i]);}printf(\n);for(i=0;iN;i++,*pp2++){b[i]=*pp2;printf(%12.6f,b[i]);}printf(\n);y1=b[0];y=0;c=1.0;for(k=1;kN;k++){for(i=k;iN;i++)t[i]=(b[i]-b[i-1])/(a[i]-a[i-k]);//t数组保存差商c=c*(x-a[k-1]);y1=y1+c*t[k];if(fabs(y-y1)e)break;y=y1;for(i=k;iN;i++)b[i]=t[i];}return(y);}