牛顿插值法

#include<stdio.h>
#include<stdlib.h>
#include<string.h>
float chashang(int x,float a[10],float c[10],int num);
float qiuzhi(int j,float a[10],float x1);
int h;
int main ()
{float a[10],b[10],c[10];//a用来储存x的值，b用来储存最终需要计算的各阶差商，
//c用来回计算的差商值 
 int i,j;int x;int num;//用来记录计算第一次差商;
  float x1,t=1,y=0; //所需要带入的求的值 
  float u;//用来记录x1的初值 
 memset(c,0,sizeof(c));
 printf("请输入多少个数");
 scanf("%d",&x);
 printf("请输入x:\n");
for(i=0;i<x;i++)
{scanf("%f",&a[i]);
}
 printf("请输入y:\n");
 for(i=0;i<x;i++)
scanf("%f",&c[i]);
h=x;
num=1;
b[0]=c[0];
for( ;x>0;x--)
c[x-1]=(c[x-1]-c[x-2])/(a[x-1]-a[x-2]);//求第一次差商 
b[1]=c[1];
j=2;// 
x=h;
 for(i=1;i<h;i++)
 {   num++;
	 chashang (x,a,c,num);
	 b[j]=c[j]; 
	 j++;
 } 
printf("输入所求x:\n"); 
scanf("%f",&x1);
u=x1;//此时num用来记录x1的初值 
 for(i=0;i<h;i++)
{
  for(j=0;j<i;j++) 
 {
  x1=qiuzhi(j,a,x1);
  t=x1*t;
  x1=u;
  if(j==0)
  t=t*b[i];
 }
 if(i==0)
 {t=b[i];
 }
 y=t+y;
 t=1;
}
printf("%f",y);
}
float qiuzhi(int j,float a[10],float x1)
{  x1=(x1-a[j]);
    return x1;
}
float chashang(int x,float a[10],float c[10],int num)
 { 
 
 for(int i=0;i<h-num;i++)
 {
	 c[x-1]=(c[x-1]-c[x-1-1])/(a[x-1]-a[x-1-num]);
	 x--;
 }

 }

太菜了别喷有点混乱