//5-1 G-S迭代法解线性方程组
【源程序】:
double temp; //判断是否主对角占优
for(i=0;itemp=0;for(j=0;jif(j!=i) temp+=fabs(A[i][j]);if(A[i][i]}cout<for(i=0;icin>>X0[i];while(1){ //使用G-S算法求解方程组
for(i=0;ifor(i=0;idouble temp1=0,temp2=0;for(j=0;jfor(j=i+1;jX0[i]=(b[i]-temp1-temp2)/A[i][i];}
double max=fabs(X0[0]-X1[0]); //迭代终止条件判断
for(i=0;iif(fabs(X0[i]-X1[i])>max)max=fabs(X0[i]-X1[i]);
if(max<=eps) break;
}
【运行结果】: