c计算4x4阶矩阵

#include<stdio.h>
int main(){
    float a[4][4];
    int i,j;
    float p=1;
    printf("请输入4*4阶数组：\n");
    for(i=0;i<4;i++){
        for(j=0;j<4;j++){
            scanf("%f",&a[i][j]);
        }
    }

    /*for(i=0;i<4;i++){
        for(j=0;j<4;j++){
            a[i][j]=i+j+1;
                    printf("%f",a[i][j]);
        }
        printf("\n");
    }*/

    printf("\n");

//计算
    for(i=0;i<3;i++){
        float k=a[i+1][0]/a[0][0];
        for(j=0;j<4;j++){
            a[i+1][j]=a[i+1][j]-(k)*a[0][j];
        }
    }
    for(i=1;i<3;i++){
        float k=a[i+1][1]/a[1][1];
        for(j=1;j<4;j++){
            a[i+1][j]=a[i+1][j]-(k)*a[1][j];
        }
    }
    for(i=0;i<4;i++){
        for(j=0;j<4;j++){
            printf("%f",a[i][j]);
        }
        printf("\n");
    }
    for(i=0;i<4;i++){
        p=p*a[i][i];
    }
    printf("计算得：%f\n",p);
    return 0;
}