Guass消元法和Guass列主元素消元法解线性方程组

#include<stdio.h>
#include<math.h>
#include<conio.h>
#define N 3
void Guass(double a[][N+1],double x[])
{
 double l[N];
 int i,j,k;
 l[0]=1.0;
 for(i=0;i<N-1;i++)
 {
  if(fabs(a[i][i])<0.000001)
  {
   printf("不能求解！/n");
   return;
  }
  for(j=i+1;j<N;j++)
   l[j]=a[j][i]/a[i][i];
  for(j=i+1;j<N;j++)
    for(k=i+1;k<N+1;k++)
     a[j][k]=a[j][k]-l[j]*a[i][k];
 }

 for(i=N-1;i>=0;i--)
 {
  x[i]=a[i][N];
  for(j=i+1;j<N;j++)x[i]=x[i]-a[i][j]*x[j];
  x[i]=x[i]/a[i][i];
 }
}
void main()
{
int i,j;
double x[N];
double s[N][N+1]=