数值方法（牛顿差值法）

#include<iostream.h>
#include<math.h>
void main()
{
char L;
do{
double M[100][100];
double x[100],y[100];
double f=1,t=0,w=1,N=0,p,R=1;
int n;
cout<<"请输入所求均差阶数:";
cin>>n;
for(int i=0;i<=n;i++)
{
cout<<"请输入x"<<i<<"的值:"<<endl;
cin>>x[i];
cout<<"请输入y"<<i<<"的值:"<<endl;
cin>>y[i];
M[i][0]=x[i];
M[i][1]=y[i];
}
for(int j=2;j<=n+1;j++)
{
for(i=1;i<=n;i++)
{
M[i][j]=(M[i][j-1]-M[i-1][j-1])/(M[i][0]-M[i-j+1][0]);
}
}
for(i=1;i<=n;i++)
{