用牛顿插值多项式求函数近似值的算法，能用C语言编程实现

#include<stdio.h>
main()
{
int i,j,k,m,z=0;
double sum=0,w=1,x,b[5][6],cc[2][4];
for(i=0;i<5;i++)
{
printf("请输入x[%d]、y[%d]:  ",i,i);
scanf("%lf%lf",&b[i][0],&b[i][1]);
}
for(j=2,k=1;j<6;j++,k++)
{
for(i=j-1;i<5;i++)
{
b[i][j]=(b[i-1][j-1]-b[i][j-1])/(b[i-k][0]-b[i][0]);
}
}

printf("enter x:");
scanf("%lf",&x);
for(m=0;m<4;m++,z=0)
{
cc[0][m]=b[m+1][m+2];
do{
w*=(x-b[z][0]);
}while(z++!=m);
cc[1][m]=w;
}
sum=b[0][1];
for(m=0;m<4;m++)
sum+=(cc[0][m]*cc[1][m]);
printf("\n差值y为: %lf\n",sum);

}