c语言实现矩阵行初等变换

#include <stdio.h>

//行与行交换
void row_exchange(int **matrix, int r1, int r2)
{
    int n = matrix[0][1];
    int a[n];
    for (int j1 = 0; j1 < n; j1++)
        a[j1] = matrix[r1][j1];
    for (int j1 = 0; j1 < n; j1++)
        matrix[r1][j1] = matrix[r2][j1];
    for (int j1 = 0; j1 < n; j1++)
        matrix[r2][j1] = a[j1];
}

//行加上n倍的行
void row_sum(int **matrix, int r1, int n, int r2)
{
    int j = matrix[0][1];
    for (int j1 = 0; j1 < j; j1++)
        matrix[r1][j1] += n * matrix[r2][j1];
}

//行减n倍的行
void row_sub(int **matrix, int r1, int n, int r2)
{
    int j = matrix[0][1];
    for (int j1 = 0; j1 < j; j1++)
        matrix[r1][j1] -= n * matrix[r2][j1];
}

//打印矩阵
void print(int **matrix)
{
    int i = matrix[0][0], j = matrix[0][1];
    for (int m = 1; m < i; m++)
    {
        for (int n = 0; n < j; n++)
            printf("%4d ", matrix[m][n]);
        printf("\n");
    }
}



int main()
{
    // m×n矩阵用m+1行n列二维数组表示，a[0][0]为数组行数，a[0][1]为数组列数
    int a[5][4] =
        {
            {5, 4, 0, 0},
            {1, 1, 1, 1},
            {1, 1, 1, 1},
            {1, 1, 1, 1},
            {3, 1, 1, 1}
        };
    int *p[a[0][1]];
    for (int i = 0; i < a[0][0]; i++)
        p[i] = a[i];
    
    print(p);
}