用C程序接线性代数方程组（全主元消去法）

#include <stdio.h>
#include <math.h>
#define N 3/*方程阶数*/

void myprintf(double a[][N],double *b);/*打印运算矩阵a,b*/
void hdai(double a[][N],double *b);/*回代求解*/
void lzyxq(double a[][N],double *b);/*高斯变换*/
void lzy(double a[][N],double *b,int k);/*行变换*/

void main()
{
double a[N][N]={3,2,1,5,6,4,7,8,9};
double b[N]={5,6,7};

lzyxq(a,b);
myprintf(a,b);

hdai(a,b);

}

void myprintf(double a[][N],double *b)
{
int i,j;
printf("a is:\n");
for(i=0;i<N;i++)
{
for(j=0;j<N;j++)
{
printf("%4.3f ",a[i][j]);
}
printf("\n");
}
printf("b is:\n");
for(i=0;i<N;i++)
{
printf("%4.3f\n",b[i]);
}
}
void lzyxq(double a[][N],double *b)
{
int i,j,k;

double l;
for(k=0;k<N-1