#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=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=0) ;
/*下一个for是解方程组*/
for(i=0;i}
void menu()
{printf("
_ _ _ _ _
");
printf(" 1。
operation
");
printf(" 2。exit");
printf("
_ _ _ _ _
");
}
main()
{int i,j,choose;
double a[N][M]={0},answer[N];
clrscr();
while(1){
leep:menu();
scanf("%d",&choose);
switch(choose){
case 1:
printf("!!The size of Maxrix is %d * %d,each line enter %d element:
",N,M,M);
for(i=0;i {printf("Enter the Matrix's %d line:
",i);
for(j=0;j scanf("%lf",&a[i][j]);
}
printf("
the corss matrix is:
_ _ _ _ _
");
gauss(a,answer);
for(i=0;i {for(j=0;j printf("%-2lf ",a[i][j]);
putchar('
');
}
printf("_ _ _ _ _
the solve is:
");
for(i=0;i break;
case 2:
exit(0);break;
default:printf("input error:
");goto leep;
}
}
getch();
}
/*试验：
西安交通大学出版社出版的《计算方法》书上28页的例2。
1:
1 2 3 -4 -2
_ -3 -4 -12 13 5
A= 2 10 0 -3 10
4 14 9 -13 7
试验结果：x1=1,x2=2,x3=3,x4=4 */。
全部
