Определить матрицы по методу Гаусса

#include <math.h>  
#include <clocale>
#include <iostream>
#include <fstream>
#include <string>
 
using namespace std;
 
void delete_array(double **array, int row_count) 
{
    for (int r = 0; r < row_count; r++)
        delete[] array[r];
    delete[] array;
}
 
int get_matrix(ifstream &file, int size, double** arr, string fileName)
{
    for (int i = 0; i < size; i++)
    {
        for (int j = 0; j < size; j++)
            file >> arr[i][j];
    }
    file.close();
    return 0;
}
 
int gaussMethod(int m, int n, double **a, double eps)
{
    int i, j, k, l;
    double r;
 
    i = 0; j = 0;
    while (i < m && j < n)
    {
        // Ищем максимальный элемент в j-м столбце,
        // начиная с i-й строки
        r = 0.0;
        for (k = i; k < m; ++k)
        {
            if (fabs(*a[k*n + j]) > r)
            {
                l = k;      // Запомним номер строки
                r = fabs(*a[k*n + j]); // и макс. эл-т
            }
        }
        if (r <= eps)
        {
            // Все элементы j-го столбца по абсолютной
            // величине не превосходят eps.
            // Обнулим столбец, начиная с i-й строки
            for (k = i; k < m; ++k)
                *a[k*n + j] = 0.0;
 
            j++;      // Увеличим индекс столбца
        }
 
        if (l != i)
        {
            // Меняем местами i-ю и l-ю строки
            for (k = j; k < n; ++k)
            {
                r = *a[i*n + k];
                *a[i*n + k] = *a[l*n + k];
                *a[l*n + k] = (-r); // Меняем знак строки
            }
        }
 
        // Обнуляем j-й столбец, начиная со строки i+1,
        // применяя элем. преобразования второго рода
        for (k = i + 1; k < m; ++k)
        {
            r = (-(*a[k*n + j]) / (*a[i*n + j]));
 
            // К k-й строке прибавляем i-ю, умноженную на r
            *a[k*n + j] = 0.0;
            for (l = j + 1; l < n; ++l)
                *a[k*n + l] += r * (*a[i*n + l]);
 
        }
 
        i++; j++;   // Переходим к следующему минору
    }
 
    return i; 
}
 
int main() 
{
    setlocale(LC_ALL, "rus");
    int m, n, i, j, rank;
    double **a;
    double eps, det;
 
    string fileName;
    cout << "Введите имя файла: ";
    cin >> fileName;
    ifstream file(fileName);
 
    file >> m >> n;
    a = new double*[n];
    for (int i = 0; i < n; i++)
        a[i] = new double[n];
 
    get_matrix(file, n, a, fileName);
    cout << "[Матрица]" << endl;
    cout << endl;
    cout << "Размерность матрицы матрицы: " << n << "x" << n << endl;
    for (int i = 0; i < n; i++)
    {
        for (int j = 0; j < n; j++)
        {
            cout << a[i][j] << "\t";
        }
        cout << endl;
    }
 
    file >> eps;
 
    rank = gaussMethod(m, n, a, eps);
 
    cout << "Ранг " << rank;
 
    det = 1.0;
    for (i = 0; i < m; ++i) 
        det *= *a[i*n + i];
        
    cout << det << endl;
 
 
    delete_array(a, n);
    return 0;  
}

int gaussMethod(int m, int n, double **a, double eps)
{
    int i, j, k, l;
    double r;
 
    i = 0; j = 0;
    while (i < m && j < n)
    {
        // Ищем максимальный элемент в j-м столбце,
        // начиная с i-й строки
        r = 0.0;
        for (k = i; k < m; ++k)
        {
            if (fabs(a[k][j]) > r)
            {
                l = k;      // Запомним номер строки
                r = fabs(a[k][j]); // и макс. эл-т
            }
        }
        if (r <= eps)
        {
            // Все элементы j-го столбца по абсолютной
            // величине не превосходят eps.
            // Обнулим столбец, начиная с i-й строки
            for (k = i; k < m; ++k)
                a[k][j] = 0.0;
 
            j++;      // Увеличим индекс столбца
        }
 
        if (l != i)
        {
            // Меняем местами i-ю и l-ю строки
            for (k = j; k < n; ++k)
            {
                r = a[i][k];
                a[i][k] = a[l][k];
                a[l][k] = (-r); // Меняем знак строки
            }
        }
 
        // Обнуляем j-й столбец, начиная со строки i+1,
        // применяя элем. преобразования второго рода
        for (k = i + 1; k < m; ++k)
        {
            r = (-(a[k][j]) / (a[i][j]));
 
            // К k-й строке прибавляем i-ю, умноженную на r
            a[k][j] = 0.0;
            for (l = j + 1; l < n; ++l)
                a[k][l] += r * (a[i][l]);
 
        }
 
        i++; j++;   // Переходим к следующему минору
    }
 
    return i; 
}