C语言解方程组程序.doc

实用标准文案 精彩文档 #includestdio.h #define M 20 #define N 21 int k,l; void main() { int i,i1,i2,i3,j,j1,j2,m,n,n1,min,min1,imax,*p,(*q)[N],a[M][N]; char c; do { printf(请输入方程个数,不大于20个:); scanf(%d,k); if(kM) printf(error!); } while(kM);//如果大于20则重新输入 do { printf(请输入未知数个数，不大于20个:); scanf(%d,l); if(lN-1) printf(error!); } while(lN-1);//如果大于20则重新输入 for(i=0;ik;i++) { printf(请输入第%d个方程系数包括等号右边的常数按下回车键进入下一步（必须是整数且各方程的第一个系数为正整数）\n请输入,i+1); for(j=0;jl;j++) { printf(第%d个未知数系数:,j+1); scanf(%d,a[i][j]); } if(j==l) { printf(请输入第%d个方程等号右边的常数:,i+1); scanf(%d,a[i][j]); } p=a[0][0]; q=a; } p=a[0][0]; q=a; for(i=0;ik;i++,q++) { for(j=0,p=*q;j=l;j++,p++) printf(%d ,*p); printf(\n); } printf(下一步:\n); for(i=0;ik;i++)//将方程按照首系数的降序排列 { m=a[i][0]; for(i1=i;i1k;i1++)//交换第一未知数系数 if(a[i1][0]m) { m=a[i1][0]; a[i1][0]=a[i][0]; a[i][0]=m; for(j=1;j=l;j++)//交换其他系数 { n=a[i1][j]; a[i1][j]=a[i][j]; a[i][j]=n; } } } for(i2=0;i2k;i2++)//化简为最简方程 { imax=1;//初始化imax for(i3=1,min=2;i3=min;i3++)//min=2是为了保证第一次执行 { for(j2=0;j2=l;j2++)//寻找不为零的系数 if(a[i2][j2]!=0) { min=a[i2][j2]; break; } if(j2l)//如果所有系数都为零 min=1; else { for(j2=0;j2=l;j2++)//找到绝对值最小的系数(零除外) { if(a[i2][j2]0) min1=-a[i2][j2]; else min1=a[i2][j2]; if(min1minmin1!=0) min=min1; } for(j2=0,n1=0;j2=ln1==0;j2++)//找到最大公约数 { n1=a[i2][j2]%i3; if(n1==0j2==l) imax=i3; } } if(i3=min/2i3!=min)//减少循环次数 i3=min-1; } for(j2=0;j2=l;j2++)//除以最大公约数 a[i2][j2]/=imax; } p=a[0][0]; q=a; for(i=0;ik;i++,q++) { for(j=0,p=*q;j=l;j++,p++) printf(%d ,*p); printf(\n); } printf(下一步:\n); for(i=0;ik;i++)//矩阵上三角化 { for(i1=i+1;i1k;i1++) if(a[i1][i]!=0) { for(j=0,m=a[i][i],n=a[i1][i];j=l;j++) a[i1][j]=a[i1][j]*m-a[i][j]*n; for(i2=0;i2k;i2++)//化简为最简方程 { imax=1;//初始化imax for(i3=1,min=2;i3=min;i3++)//min=2是为了保证第一次执行 { for(j2=