列主元消元法解线性方程组.docVIP

列主元消元法解线性方程组（已通过验证） PXEOA ?#define ? N ? 100 opww qP ?#includestdio.h Ki~Y0Mu ?#includemath.h |_NVsz R h ?main( ) {4 :1MC! ?{ S*}X7cRu ?int n, i,j,row=0,k=0,m=0; ]2sY%`Pb ?char s[2]; GTQi t@ ?? float a[N][N],b[N],sum=0.0,max=0.0,mik=0.0; zu x){ ?? float t1; @xWL]q ?? clrscr( ); HqF[Ojw ?? printf(please input fang cheng zu de jie shu zhi:\n); Den 1 rK4j ?? scanf(%d,k); rXN#ziX ?? n=k; ? // 适应阶数可变 // 9(.O2vl ?? printf(please input xi shu de zhi:\n); D4U,r 36 ?? for(i=1;i=n;i++) ? // 按行输入 // I}ir:/,L ?? ? for(j=1;j=n;j++) h.UP#@P ?? ? ? ? ? scanf(%f,a[j]); 8G0@lJ0 ?printf(\ncheck the xi shu:\n); @1YR`,W ?for(i=1;i=n;i++) xwE9j
?for(j=1;j=n;j++)// 形成系数距阵输入便于检查 // 3LXE|T?$ ?? ? ? ? ? ? { ? printf(%4.3f\t,a[j]); lkSe{e6bke ?if(!(++m%k))printf(\n); ./a` ?? ? ? ? ? ? } jTy0 .Q ?printf(输入方程组右边的值:\n); px[(?j6!R ?? for(i=1;i=n;i++) // 从上往下依次录入右边的值 // D:[0M ?? ? ? ? scanf(%f,b); ! ?c
)PZ ?? printf(\ncheck you bian de zhi:\n); W,f`P%SQ ?? for(j=1;j=n;j++) $IWV[3P+Y ?? ? ? printf(%4.3f\n,b[j]);// 检验输入工作完成 // Ev9p#%O ?// 下面是求每列的最大值并实现交换行，k控制消元步骤 // ?ER?w@-i ?for(k=1;kn;k++) y28NA^uk` ?? { ? max=a[k][k]; Sdib=,K.8L ?? ? { ? for(j=k+1;j=n;j++) 6Hg/fr ?? ? ? ? ? if(fabs(a[j][k])=fabs(max)) iN%JsQz ?? ? ? ? ? ? { ]*m8tnb F? ?? ? ? ? ? ? ? max=a[j][k]; o?c` ?? ? ? ? ? ? ? row=j;//记录最大元素的行标// \{w*]@)V=N ?? ? ? ? ? ? } aYSo]}m ?? ? ? ? if(fabs(max)==0) //检查方程组是否可解// of*LUmz ?{ ? printf(error\n); 9/UZzc=O ?exit(1); 55.sZZ ?} u{P/\qS` ?? ? ? ? else d9#;}g _ ?? ? ? ? ? ? for(i=k;i=n;i++) qz{iMl ?? ? ? ? ? ? { ? t1=a[k]; IseB mw ?? ? ? ? ? ? ? a[k]=a[row]; (7c104j  ?? ? ? ? ? ? ? a[row]=t1; -{qp^x-4 ?? ? ? ? ? ? } zzaoFki[ ?? ? ? ? t1=b[k]; :Fu6buq ?? ? ? ? b[k]=b[row]; {d^ AB)I ?? ? ? ? b[row]=t1; // 实现右边的交换 // )f+f=92 ?} vE *Pi9gh ?//下面是计算每列中的乘法因子并计算a[I][k],a[j],b[I]// @9;Obzfd ?? for(i=k+1;i=n;i++) }C$ /{5Z ?? ? { ? mik=a[k]/a[k][k]; kikWZI4 ??