ssadfasdfsdaf.txt

牛顿插值计算牛顿插值计算牛顿插值计算牛顿插值计算牛顿插值计算牛顿插值计算牛顿插值计算牛顿插值计算牛顿插值计算牛顿插值计算

#include <iostream.h> 
#include <math.h> 
void main() 
{ 
char L; 
do 
{ 
double M[100][100]; 
double x[100],y[100]; 
double X=1,xx=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++) 
{ 
cout<<"其"<<i<<"阶均差为："<<M[i][i+1]<<endl; 
} 
cout<<"请输入x的值：x="; 
cin>>xx; 
for(i=0;i<n;i++) 
{ 
X*=xx-x[i]; 
N+=M[i+1][i+2]*X; 
P=M[0][1]+N; 
} 
cout<<"其函数值：y="<<P<<endl; 
for(i=0;i<n;i++) 
{ 
w*=xx-x[i]; 
R=fabs(M[n][n+1]*w); 
} 
cout<<"其截断误差：R="<<R<<endl; 
cout<<endl<<"还想算其它插值吗？是请按'y'否则按'n'"<<endl; 
cin>>L; 
}while(L=='y'); 
}