Gauss消元法

#include "Stdio.h"
#include "Conio.h"
/*L是矩阵的行减1,从程序上看是最外层循环的次数
N 对应矩阵的行数，M对应矩阵的列数
可以通过改变L、N、M来控制矩的阶数 */
#define L 3
#define N 4
#define M 5
void gauss(double a[N][M],double x[N])
{int i,j,l,n,m,k=0;
double temp[N];
/*第一个do-while是将增广矩阵消成上三角形式*/
do{n=0;
for(l=k;l<L;l++)temp[n++]=a[l+1][k]/a[k][k];
for(m=0,i=k;i<N;i++,m++)
for(j=k;j<M;j++)a[i+1][j]-=temp[m]*a[k][j];
k++;
}while(k<N) ;
/*第二个do-while是将矩阵消成对角形式，并且重新给k赋值*/
k=L-1;
do{n=0;
for(l=k;l>=0;l--)temp[n++]=a[k-l][k+1]/a[k+1][k+1];
for(m=0,i=k;i>=0;i--,m++)
for(j=k;j<M;j++)a[k-i][j]-=temp[m]*a[k+1][j];
k--;
}while(k>=0) ;
/*下一个for是解方程组*/
for(i=0;i<N;i++)x[i]=a[i][N]/a[i][i];

}
void menu()
{printf("\n _ _ _ _ _\n");
printf(" 1.operation\n");
printf(" 2.exi