迭代法的程序.docVIP

#includeiostream#includemath.h#includeiomanipusing namespace std; #define kk 50 //定义最大方程元数int n,i,c,j,ll,hh,gg,mm;double A[kk][kk],x[kk][kk],b[kk],y[kk],a[kk],z[kk],m,nn,d,e=1,w,fff ; void main(){cout　　　**********************************************************************endl;cout　　　***本程序可以用雅可比迭代法,塞德尔迭代法,逐次超松弛法求解线性方程组***endl;cout　　　**********************************************************************endl;//*********************************数据的输入********************************/*bb:*/cout输入的方程元数endl;cinn;cout请输入方程系数矩阵:endl;for(i=0;in;i++)?? for(j=0;jn;j++)??? cinA[i][j];cout请输入右边向量:endl;for(i=0;in;i++)?? cinb[i];cout输入你想要的迭代精度(建议1e-5以上)！endl;cinfff;cout输入最大迭代次数（建议300次以上）！endl;cinmm; ////*******************************判断是否对角占优*************************//??? for(i=0;in;i++)// {//?? for(j=0;jn;j++)//?? {//??? nn=0;//??? if(i==j)//??? {//???? d=fabs(A[i][i]);//??? }//??? else //???? nn=nn+fabs(A[i][j]);//?? }//?? if(nnd)//?? {//??? cout该方程不对角占优，迭代不收敛endl;//??? cout是否继续？是（0），否（1）endl;//??? cinhh;//??? if(hh!=1)//???? goto bb;//??? else exit(1);//?? }// }//********************************计算出迭代矩阵*************************for(i=0;in;i++){?? b[i]=b[i]/A[i][i];???? for(j=0;jn;j++)?? {??? if(i==j)??? {???? x[i][i]=0;??? }??? else??? {???? x[i][j]=-A[i][j]/A[i][i];??? }?? } } //*******************************输出迭代矩阵*****************************cout计算出迭代矩阵为：endl;for(i=0;in;i++){?? for(j=0;jn;j++)??? coutx[i][j]??? ;?? coutb[i]?? ;?? coutendl; } //****************************迭代方法的选择*****************************cout请你选择迭代方法！endl;coutendl;coutendl;cout选用雅可比迭代法，请输入（1）！endl;coutendl;coutendl;cout选用塞德尔迭代法，请输入（2）！endl;coutendl;coutendl;cout选用逐次超松弛法，请输入（3）！endl;coutendl;coutendl;cinll;//*****************************赋迭代初值***********************************cout输入迭代初值endl;for(i=0;in;i++)?? ciny[i];int f=1;switch(ll){case 1:goto cc;?? break;case 2:goto aa;?? break;case 3:goto dd;}//***********