输入系数矩阵后

/* Note:Your choice is C IDE */
#include "stdio.h"
#include H 3
#include L 4
#include"math.h"
void main()
{int i,j,k,m,n;
	static float a[H][L],x[H],c[H],s1;
	for(i=0;i<=H-1;i++)
      for(j=0;j<=L-1;j++)
        scanf("%f",a[i][j]);
    for(i=0;i<=H-2;i++)
    {  for(j=i;j<=H-2;j++)
    	{k=i;
    	 if(fabs(a[k][i])<fabs(a[j+1][i]))
    	 k=j;}
        for(m=0;m<=H-2;m++)
        {s1=a[i][m];a[i][m]=a[k][m];a[k][m]=s1;	}
        for(m=0;m<=H-2;m++)
          for(n=0;n<=L-1;n++)
            a[m+1][n]=a[m+1][n]-a[m+1][n]*a[i][n]/a[i][i];
    }
    printf("消去后的上三角矩阵是:");
    for(m=0;m<=H-1;m++)
      {for(n=0;n<=L-1;n++)
      	printf("%f  ",a[m][n]);
      	printf("\n");
      }
      x[H-1]=a[H-1][L-1]/a[H-1][L-2];
      	for(i=H-2;i>=0;i--)
           {for(j=L-2;j>i;j--)
           	c[i]=c[i]+a[i][j]*x[j];
           	x[i]=(a[i][L-1]-c[i])/a[i][i];}
       printf("The answer is:\n");
       for(m=0,i=1;m<=H-1;m++,i++)
       printf("x[%d]=%f\n",i,x[m]);    
}