C语言编写牛顿插值多项式

我编的程序 可运行答案不对，请问错在哪啊
#include <stdio.h>
void difference(float a[],float b[],int n)
{ float g[5];
int i,j;
for(i=1;i<n;i++)
{ g[0]=b[i];
for(j=0;j<i;j++)
g[j+1]=(g[j]-b[j])/(a[i]-a[j]);
b[i]=g[i];
}
return;
}
main()
{ int i,n=4;
float a[]={0,1,1.5,2},b[]={1,1.25,2.5,5.5};
float x=1.75,y;
difference(a,b,n);
y=b[4];
for(i=n-1;i>=0;i--)
y=y*(x-a[i])+b[i];
printf("y=%f",y);
}

牛顿的插值法用C语言怎么编写怎么编啊? 求，用c语言编写的一个程序，编写程序，因式分解多项式 c语言编写x3+x2+3x=0用牛顿迭代法求 C语言程序 牛顿迭代法 求c语言程序 一元多项式 多项式相加C语言实现（数据结构） 编写C语言程序 c语言编写应用程序 C语言怎样编写？ 用C语言编写：