牛顿n阶插值

main()
{
  int i,j;
  float x[9];
  float cc[9][9];
  float m,n,xx,w,q,t,s;
  int a;
  printf("Input new num:/n");
  scanf("%d",&a);
  a=a+1;
  printf("/n");
  printf("Input new x:/n");
  for(i=0;i<a;i++)
   {
    scanf("%f",&x[i]);
    
   }
  printf("Input new y value:/n");
  for(i=0;i<a;i++)
  {
    scanf("%f",&cc[i][0]);
   
   }
 
  for(j=0;j<a;j++)
     {printf("/n");
      printf("%f",cc[j][0]);
     }
for(j=1;j<a;j++)
    for(i=j;i<a;i++)
         {m=cc[i][j-1]-cc[i-1][j-1];
          n=x[i]-x[i-j];
          cc[i][j]=m/n;
          }

  for(i=1;i<5;i++)
       {  printf("/n");
          printf("%f",cc[i][i]);
       }
  printf("/n");
  printf("Input a new value/n");
  scanf("%f",&xx);
    s=cc[0][0];
    t=cc[0][0];
  for(i=1;i<a;i++)
        {  w=xx-x[i-1];
           q=cc[i][i]/cc[i-1][i-1];
           t=t*q*w;
           s=s+t;
        }
  printf("/n");
  printf("%f",s);
 
 }