作业:①
分别用 J 法和 G-S 法求解以下方程,并讨论结果。
1
2
3
1 2 2 1
1 1 1 1
2 2 1 1
x
x
x
-
æ öæ ö æ ö
ç ÷ç ÷ ç ÷
=
ç ÷ç ÷ ç ÷
ç ÷ç ÷ ç ÷
è øè ø è ø
#include<iostream>
using namespace std;
//J 法解线性方程
int main(){
int m,n,i,j,times=0,mtimes;
double s,sum,max;
cout<<"请输入系数矩阵行数 m、列数 n:"<<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<<"请输入系数矩阵 A:"<<endl;
for(i=0;i<m;i++)
for(j=0;j<n;j++)
cin>>A[i][j];
cout<<"请输入常数向量 B:"<<endl;
for(i=0;i<m;i++)
cin>>B[i];
cout<<"请输入最大允许误差 s:"<<endl;
cin>>s;
cout<<"请输入最大迭代次数:"<<endl;
cin>>mtimes;
cout<<"请输入一零级向量 X:"<<endl;
for(i=0;i<m;i++){
cin>>X[i];
T[i]=X[i];//T[]存放上一次迭代结果
}
- 1
- 2
- 3
前往页