埃特金逐步线性插值

代码如下：#includestdio.h#includemath.h#includemalloc.h#defineN5main(){inti,j;floatt,e,b;floatp[N][N],x[N],y[N];printf(请输入插值点t与拟合精度e：\n);scanf(%f%f,&t,&e);printf(请输入横坐标值：\n);for(i=0;iN;i++)scanf(%f,&x[i]);printf(请输入纵坐标值：\n);for(i=0;iN;i++)scanf(%f,&y[i]);for(i=0;iN;i++)p[i][0]=y[i];for(i=1;iN;i++)p[i][1]=(t-x[i])*y[0]/(x[0]-x[i])+(t-x[0])*y[i]/(x[i]-x[0]);for(i=2;iN;i++){for(j=2;ji+1;j++)p[i][j]=p[j-1][j-1]+(t-x[j-1])*(p[j-1][j-1]-p[i][j-1])/(x[j-1]-x[i]);}for(i=0;iN;i++){for(j=0;jN;j++)printf(%8.6f,p[i][j]);printf(\n);}for(j=1;jN;j++){for(i=j+1;iN;i++)if(fabs(p[i][j]-p[i-1][j])e)printf(%f,floor(p[i-1][j]/e)*e);}}1.插值点为0.512，精度为0.0012.插值点为0.512，精度为0.00013.插值点为0.462，精度为0.0014.插值点为0.462，精度为0.0001