作业:①
分别?/p>
J
法和
G
-
S
法求解下列方程,并讨论结果?/p>
1
2
3
1
2
2
1
1
1
1
1
2
2
1
1
x
x
x
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?/p>
?
#include<iostream>
using namespace std;
//J
法解线性方?/p>
int main(){
int m,n,i,j,times=0,mtimes;
double s,sum,max;
cout<<"
请输入系数矩阵行?/p>
m
、列?/p>
n
?/p>
"<<endl;
cin>>m>>n;
if(m<n)
cout<<"
方程组无唯一解!
"<<endl;
double **A=new double *[m];
for (i=0; i<m; i++) {
A[i] =new double [n];
}
double *B=new double [m];
double *X=new double [m];
double *T=new double [m];
double *S=new double [m];
cout<<"
请输入系数矩?/p>
A
?/p>
"<<endl;
for(i=0;i<m;i++)
for(j=0;j<n;j++)
cin>>A[i][j];
cout<<"
请输入常数向?/p>
B
?/p>
"<<endl;
for(i=0;i<m;i++)
cin>>B[i];
cout<<"
请输入最大允许误?/p>
s:"<<endl;
cin>>s;
cout<<"
请输入最大迭代次?/p>
:"<<endl;
cin>>mtimes;
cout<<"
请输入一零级向量
X:"<<endl;
for(i=0;i<m;i++){
cin>>X[i];
T[i]=X[i];//T[]
存放上一次迭代结?/p>
}