【文件属性】:
文件名称:数值计算方法实验报告
文件大小:554KB
文件格式:DOC
更新时间:2013-12-19 14:47:19
计算方法
很好的自己编的都运行出来了
#include
#include
const int m=10;
void main()
{
int choice=1;
while(choice==1)
{
double a[m][m],b[m],e,x[m],y[m],w,se,max;
int n,i,j,N,k;
cout<<"\n请输入方程的个数n:";
cin>>n;cout<>a[i][j];
}
cout<<"\n请输入各个方程等号右边的常数项。\n";
for(i=1;i<=n;i++)
{
cin>>b[i];
}
cout<<"\n请输入最大迭代次数:";
cin>>N;
cout<<"\n请输入最大偏差:";
cin>>e;
for(i=1;i<=n;i++)
{
x[i]=0;
y[i]=x[i];
}
k=0;
while(k!=N)
{
k++;
for(i=1;i<=n;i++)
{
w=0;
for(j=1;j<=n;j++)
{
if(j!=i)
w=w+a[i][j]*y[j];
}
y[i]=(b[i]-w)/double(a[i][i]);
}
max=fabs(x[1]-y[1]);
for(i=1;i<=n;i++)
{
se=fabs(x[i]-y[i]);
if(se>max)
max=se;
}
if(max>choice;
}
}