高斯消元法选编.docVIP

#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;lL;l++)temp[n++]=a[l+1][k]/a[k][k];??for(m=0,i=k;iN;i++,m++)??for(j=k;jM;j++)a[i+1][j]-=temp[m]*a[k][j];??k++;??}while(kN)?;??/*第二个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;jM;j++)a[k-i][j]-=temp[m]*a[k+1][j];??k--;??}wh