C语言解线性方程组.doc

//求线性方程组的解,化成除对角线外其余的元素都为0#include<>intmain(){inti,j,n,m,r;floata[50][50],k,det=1;A:printf("输入方程组个数m=");scanf("%d",&m);printf("输入未知数个数n=");scanf("%d",&n);if(m==n){printf("输入线性方程组的系数和常数组成的矩阵\n"); //输入矩阵for(i=0;i<m;i++){ for(j=0;j<n+1;j++) scanf("%f",&a[i][j]);} printf("您输入的矩阵为:\n");for(i=0;i<m;i++){ for(j=0;j<n+1;j++) printf("%\t",a[i][j]); printf("\n"); } //用极小数代替矩阵中的0元素 for(i=0;i<m;i++) { for(j=0;j<n;j++) { if(a[i][j]==) a[i][j]=1e-10; } } //求解第一步化为下上角矩阵 for(i=1;i<m;i++) { for(j=0;j<i;j++) { k=-a[i][j]/a[j][j]; //printf("%\t",k); for(r=0;r<n+1;r++) { a[i][r]=k*a[j][r]+a[i][r]; } } } //求行列式 for(i=0;i<m;i++) det*=a[i][i]; printf("行列式:det=%\n",det); if(det!=0) { printf("上三角矩阵如下:\n"); for(i=0;i<m;i++) { for(j=0;j<n+1;j++) printf("%\t",a[i][j]); printf("\n"); } //求解第二步化为上三角矩阵 printf("消去右上角元素:\n"); for(i=m-2;i>=0;i--) { for(j=n-1;j>=i+1;j--) { k=-a[i][j]/a[j][j]; for(r=0;r<n+1;r++) { a[i][r]=k*a[j][r]+a[i][r]; } } } for(i=0;i<m;i++) { for(j=0;j<n+1;j++) { printf("