新手上路，请多多帮助

#include <stdio.h>
#include <stdlib.h>
#define M 4
#define N 4
#define P 1

int main()
{
    float A[M][N],B[N][P],c[M][N+P];
    int i,j;
    float e;
    float x1,x2,x3,x4;
    printf("请输入%d*%d阶矩阵:\n",M,N);
    for(i=0;i<M;i++)
        for(j=0;j<N;j++)
    scanf("%f",&A[i][j]);
    printf("请输入%d*%d阶矩阵:\n",N,P);
    for(i=0;i<N;i++)
        for(j=0;j<P;j++)
    scanf("%f",&B[i][j]);
    for(i=0;i<M;i++)
    {
      for(j=0;j<N;j++)
        c[i][j]=A[i][j];
      for(j=4;j==4;j++)
        c[i][j]=B[i][0];

    }

    if(c[0][0]=0)
    {
        for(i=1;i<M&&c[0][0]==0;i++)
        for(j=0;j<N+P;j++)
          {
             e=c[0][j];
             c[i][j]=c[0][j];
             c[0][j]=c[i][j];
          }
    }


    if(c[0][0]!=0)
    {
         for(i=0,j<M;i<=j;j++)
        {
            c[i][j]=c[i][j]/c[0][0];
        }
        for(i=1,j<N+P;i<=j;j++)
        {
            c[i][j]=c[i][j]-c[i][j]*c[0][j];
        }
    }

    if(c[1][1]=0)
    {
        for(i=2;i<M&&c[1][1]==0;i++)
        for(j=1;j<N+P;j++)
          {
             e=c[1][j];
             c[i][j]=c[1][j];
             c[1][j]=c[i][j];
          }

    }

    if(c[1][1]!=0)
    {
        for(i=1,j<N+P;i<=j;j++)
        {
            c[i][j]=c[i][j]/c[1][1];
        }
        for(i=2,j<N+P;i<=j;j++)
        {
            c[i][j]=c[i][j]-c[i][j]*c[1][j];
        }
    }

    if(c[2][2]=0)
    {
        for(i=3;i<M&&c[2][2]==0;i++)
        for(j=2;j<N+P;j++)
          {
             e=c[2][j];
             c[i][j]=c[2][j];
             c[2][j]=c[i][j];
          }
    }

    if(c[2][2]!=0)
    {
       for(i=2,j<N+P;i<=j;j++)
        {
            c[i][j]=c[i][j]/c[2][2];
        }
        for(i=2,j<N+P;i<=j;j++)
        {
            c[i][j]=c[i][j]-c[i][j]*c[2][j];
        }
    }


    if(c[3][3]!=0)
    {
        x4=c[3][4]/c[3][3];
        x3=(c[2][4]-c[2][3]*x4)/c[2][2];
        x2=(c[1][4]-c[1][3]*x4-c[1][2]*x3)/c[1][1];
        x1=(c[0][4]-c[0][3]*x4-c[0][2]*x3-c[0][1]*x2)/c[0][0];
    }

    printf("%f\n%f\n%f\n%f\n",x1,x2,x3,x4);
    return 0;
}
试了这个程序，可是好像不太对啊，求大神改改或者来个新程序呗~