牛顿向前差分插值

#include "iostream.h"
#include "math.h"
#define SIZE 50

double jc(int n)
{
double t=1,i;
for(i=1;i<=n;i++)
t=t*i;
return t;
}

void main()
{
double x0,x,y,f[SIZE][SIZE],h,t;
int n,i,j;
cout<<"\n请输入插值节点的个数:";
cin>>n;
cout<<"\n请输入步长:";
cin>>h;
cout<<"\n请输入X[0]:";
cin>>x0;
cout<<"\n请输入插值节点的纵坐标:";
for(i=0;i<=n-1;i++)
{
cout<<"\nx["<<i<<"]="<<(x0+i*h)<<",y["<<i<<"]=";
cin>>f[i][0];
}
cout<<"\n请输入插值点";
cin>>x;