您的当前位置：首页雅克比迭代法和高斯-赛德尔法解线性方程组(C++)

雅克比迭代法和高斯-赛德尔法解线性方程组(C++)

来源：化拓教育网

作业：①

分别用J法和G-S法求解下列方程，并讨论结果。

122x11111x21221x13

#include

using namespace std;

//J法解线性方程

int main(){

int m,n,i,j,times=0,mtimes;

double s,sum,max;

cout<<\"请输入系数矩阵行数m、列数n：\"<cin>>m>>n;

if(mcout<<\"方程组无唯一解！\"<double **A=new double *[m];

for (i=0; iA[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：\"<for(i=0;ifor(j=0;jcin>>A[i][j];

cout<<\"请输入常数向量B：\"<for(i=0;icin>>B[i];

cout<<\"请输入最大允许误差s:\"<cin>>s;

cout<<\"请输入最大迭代次数:\"<cin>>mtimes;

cout<<\"请输入一零级向量X:\"<for(i=0;icin>>X[i];

T[i]=X[i];//T[]存放上一次迭代结果

}

do{

for(i=0;isum=0;

for(j=0;jif(j!=i)

sum+=A[i][j]*T[j];

X[i]=(B[i]-sum)/A[i][i];

S[i]=(T[i]-X[i])*(T[i]-X[i]);

}

for(i=0;iT[i]=X[i];

times++;

max=S[0];

for(i=1;iif(S[i]>max)

max=S[i];

}while((max>s*s)&&(timescout<<\"该方程组的解为：\"<for(i=0;icout<cout<for (i=0; idelete[] A[i];

A[i]=NULL;

}

delete[] A; A=0;

delete[] B; B=0;

delete[] T; T=0;

delete[] X; X=0;

delete[] S; S=0;

return 0;

}