Программа, приводящая матрицу к ступенчатому виду

#include <stdio.h>
#include <stdlib.h>
#include <locale.h>
#include <math.h>
 
void detFind(int n, int m, double *a, double eps);
 
int main()
{
    setlocale(LC_ALL, "Rus");
    double eps = 0.0000001; //Погрешность
    double *a;
    int n = 3; //Кол-во строк
    int m = 3; //Кол-во столбцов
    a = (double *) malloc(n * m * sizeof(double));
 
    //Вводим массив
    for(int k = 0; k < n; k++)
    {
        for(int l = 0; l < m; l++)
        {
            printf("Введите элемент массива a[%d][%d]: ", k+1, l+1);
            scanf("%lf", &a[k*m + l]);
        }
    }
    system("cls");
    //Выводим массив
    printf("Исходный массив:\n");
    for(int k = 0; k < n; k++)
    {
        for(int l = 0; l < m; l++)
        {
            printf("%lf\t", a[k*m + l]);
        }
        printf("\n");
    }
 
    detFind(n, m, a, eps);
        //Отладка
    printf("\nМассив, приведённый к треугольному виду:\n");
    for(int z = 0; z < n; z++)
    {
        for(int y = 0; y < m; y++)
        {
            printf("%lf\t", a[z*m + y]);
        }
        printf("\n");
    }
    return 0;
}
 
void detFind(int n, int m, double *a, double eps)
{
    int i, j, k, l;
    double max = 0.0; int maxN;
    double r;
    i = j = 0;
 
    while(i < n && j < m)
    {
        //Поиск max по модулю в j-столбике
        for(k = i; k < n; k++)
        {
            if(fabs(a[k*m + j]) > fabs(max))
            {
                max = a[k*m + j]; //Запоминаем max значение
                maxN = k; //Запоминаем строку с max значением
            }
        }
        //Если max < eps => то все элементы столбца = 0 и мы переходим к следующей итерации, т.е. столбику j+1
        if(fabs(max) < eps)
        {
            for(k = j; k < n; k++)
            {
                a[k*m + j] = 0.0;
            }
            j++;
            continue;
        }
 
        //Меняем местами i-строчку с maxN
        for(l = 0; l < m; l++)
        {
            r = a[i*m + l]; //Буферная переменная равна элементы из i строчки
            a[i*m + l] = a[maxN*m + l]; //элемент i-строчки равен maxN
            a[maxN*m + l] = -r; //элемент maxN-строчки равен знач буферной пер, т.е. знач i строчки, с противоп знаком
        }
        //Обнуляем столбец j
        for(k = i + 1; k < n; k++)
        {
            r = -(a[k*m + j]/a[i*m + j]);
            a[k*m + j] = 0.0;
            for(l = j + 1; l < m; l++)
            {
                a[k*m + l] += r*a[i*m + l];
            }
        }
 
    i++;
    j++;
    }
 
return;
}