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

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

作业：①

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

#includeiostream

usingnamespacestd;

//J法解线性方程

intmain(){

intm,n,i,j,times=0,mtimes;

doubles,sum,max;

cout请输入系数矩阵行数m、列数n：endl;

cinmn;

if(mn)

cout方程组无唯一解！endl;

double**A=newdouble*[m];

for(i=0;im;i++){

A[i]=newdouble[n];

}

double*B=newdouble[m];

for(i=0;im;i++){

sum=0;

for(j=0;jn;j++)

if(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;im;i++)

T[i]=X[i];

times++;

max=S[0];

for(i=1;im;i++)

if(S[i]max