求助 这到底哪里错了

#include<stdio.h>
#define N 2
#define M 3
void xsj(float a[N][M],float b[N]);
void main()
{float a[N][M];
float x[N];
int i,j;
printf("请输入方程组的增广矩阵：\n");
for(i=0;i<=N;i++)
{for(j=0;j<=M;j++)
scanf("%f",&a[i][j]);}
xsj(a,x);
printf("方程组的解为：\n");
for(i=0;i<=N;i++)
printf("x[%d]=%f\n",i+1,x[i]);
}
void xsj(float a[N][M],float b[N])
{float sum=0;
int i,j;
b[0]=a[0][M]/a[0][0];
for(i=1;i<=N;i++)
{for(j=0;j<i;j++)
sum+=(b[j]*a[i][j]);
b[i]=(a[i][M]-sum)/a[i][i];
sum=0;}
}