не получается составить программу

/*
4,13x1-2,87x2-1,94x3+0,61x4= 0,32,
1,27x1+7,23x2-0,15x3+1,71x4=-4,16,
0,19x1+2,75x2+3,14x3-0,76x4= 2,33;
2,87x1+4,33x2-2,41x3-3,42x4= 2,79;
    Метод Гаусса — Жордана
1. Выбирают первый слева столбец матрицы, в котором есть хоть одно отличное от нуля значение.
2. Если самое верхнее число в этом столбце есть ноль, то меняют всю первую строку матрицы с другой строкой матрицы, где в этой колонке нет нуля.
3. Все элементы первой строки делят на верхний элемент выбранного столбца.
4. Из оставшихся строк вычитают первую строку, умноженную на первый элемент соответствующей строки, с целью получить первым элементом каждой строки (кроме первой) ноль.
5. Далее проводят такую же процедуру с матрицей, получающейся из исходной матрицы после вычёркивания первой строки и первого столбца.
6. После повторения этой процедуры n-1 раз получают верхнюю треугольную матрицу
7. Вычитаем из предпоследней строки последнюю строку, умноженную на соответствующий коэффициент, с тем, чтобы в предпоследней строке осталась только 1 на главной диагонали.
8. Повторяют предыдущий шаг для последующих строк. В итоге получают единичную матрицу и решение на месте свободного вектора (с ним необходимо проводить все те же преобразования).
9. Чтобы получить обратную матрицу, нужно применить все операции в том же порядке к единичной матрице.
*/
 
#include <iostream>
#include <iomanip>
using namespace std;
 
const int n = 4;
double arr[n][n+1] = { {4.13, -2.87, -1.94,  0.61,  0.32},
                       {1.27,  7.23, -0.15,  1.71, -4.16},
                       {0.19,  2.75,  3.14, -0.76,  2.33},
                       {2.87,  4.33, -2.41, -3.42,  2.79} };
 
int firstLeftColumn(int);   // Выбор первого слева столбца матрицы
int firstTopRow(int);       // Выбор перовой сверху строки не содержащей 0
void swapLines(int, int);   // Обмен строк
void calcArr(int);          // вычисление
void inputArr();            // вывод массив в консоль
 
void main()
{
    setlocale(LC_ALL, "Russian");
    inputArr();
    
    // 5. Далее проводят такую же процедуру с матрицей, получающейся из исходной матрицы
    // после вычёркивания первой строки и первого столбца.
    // 6. После повторения этой процедуры n − 1 раз получают верхнюю треугольную матрицу
    for (int i = 0; i < n; i++)
        calcArr(i); 
    
    // 7. Вычитаем из предпоследней строки последнюю строку, умноженную на соответствующий
    // коэффициент, с тем, чтобы в предпоследней строке осталась только 1 на главной диагонали.
    double temp;
    for (int k = n - 1; k > 0; k--) {
        for (int i = k - 1; i >= 0; i--) {  // цикл по предпоследним строкам
            temp = arr[i][k];
            for (int j = 0; j < n + 1; j++) // цикл по последним строкам
                arr[i][j] -= arr[k][j] * temp;
        }
    }
 
    // 8. Повторяют предыдущий шаг для последующих строк. В итоге получают единичную матрицу
    // и решение на месте свободного вектора (с ним необходимо проводить все те же преобразования).
 
    // 9. Чтобы получить обратную матрицу, нужно применить все операции в том же порядке  
    // к единичной матрице.
 
    inputArr();
    for (int i = 0; i < n; i++)
        cout << "x" << i + 1 << "=" << arr[i][n] << endl;
 
    system("pause");
}
 
int firstLeftColumn(int p)
{
    for (int j = p; j < n; j++)
        for (int i = p; i < n; i++)
            if(arr[i][j])
                return j;
}
 
int firstTopRow(int p, int column)
{
    for (int i = p + 1; i < n; i++)
        if(arr[i][column])
            return i;
}
 
void swapLines(int p, int line)
{
    double temp;
    for (int j = 0; j < n + 1; j++) {
            temp = arr[p][j];
            arr[p][j] = arr[line][j];
            arr[line][j] = temp;
        }
}
 
void calcArr(int p)
{   
    // 1. Выбрать первый слева столбец матрицы, в котором есть хоть одно отлично от нуля значение
    int numCol = firstLeftColumn(p);    // номер столбца
    
    // 2. Если самое верхнее число в этом столбце есть ноль, то меняют всю первую строку матрицы
    // с другой строкой матрицы, где в этой колонке нет нуля
    if (!arr[p][numCol])
        swapLines(p, firstTopRow(p, numCol));
 
    // 3. Все элементы первой строки делят на верхний элемент выбранного столбца
    double temp = arr[p][numCol];
    for (int j = 0; j < n + 1; j++)
        arr[p][j] /= temp;
    
    // 4. Из оставшихся строк вычитают первую строку, умноженную на первый элемент соответствующей
    // строки, с целью получить первым элементом каждой строки (кроме первой) ноль.
    for (int i = p + 1; i < n; i++) {
        temp = arr[i][p];
        for (int j = 0; j < n + 1; j++)
            arr[i][j] -= arr[p][j] * temp;
    }
}
 
void inputArr()
{
    cout << endl;
    for (int i = 0; i < n; i++) {
        for (int j = 0; j < n + 1; j++) {
            if (j == n)
                cout << " |";
            cout << setw(6) << arr[i][j];
        }
        cout << endl;
    }
}

#include <iostream>
#include <ctime>
#include <cstdlib>
#include <iomanip>
 
int** mult_matr(int** first, int** second, int n)
{
    int** Arr;
    Arr=new int*[n];
    for(int i=0; i<n; ++i)
        Arr[i]=new int[n];
    for(int i=0; i<n; ++i)
    {
        for(int j=0; j<n; ++j)
        {
            Arr[i][j]=0;
            for(int k=0; k<n; ++k)
            {
                Arr[i][j]+=first[i][k]*second[k][j];
            }
        }
    }
    return Arr;
}
 
void print(int** Arr, int n)
{
    for(int i=0; i<n; ++i)
    {
        for(int j=0; j<n; ++j)
            std::cout<<std::setw(5)<<Arr[i][j]<<' ';
        std::cout<<'\n';
    }
    std::cout<<'\n';
}
 
void fill_rand(int** Arr, int n)
{
    int first, last;
    std::cout<<"Enter ranges for random fill: ";
    std::cin>>first>>last;
    for(int i=0; i<n; ++i)
    {
        for(int j=0; j<n; ++j)
        {
            Arr[i][j]=first+rand()%(last-first);
        }
    }
}
 
int main()
{
    srand(static_cast<unsigned>(time(0)));
    int** Arr;
    int n;
    std::cout<<"Enter n: ";
    std::cin>>n;
    Arr=new int*[n];
    for(int i=0; i<n; ++i)
        Arr[i]=new int[n];
    fill_rand(Arr, n);
    std::cout<<"A\n";
    print(Arr, n);
    int** secArr=mult_matr(Arr, Arr, n);
    std::cout<<"A^2\n";
    print(secArr, n);
    int** thirdArr=mult_matr(secArr, Arr, n);
    std::cout<<"A^3\n";
    print(thirdArr, n);
    for(int i=0; i<n; ++i)
        delete[] Arr[i];
    delete[] Arr;
    for(int i=0; i<n; ++i)
        delete[] secArr[i];
    delete[] secArr;
    for(int i=0; i<n; ++i)
        delete[] thirdArr[i];
    delete[] thirdArr;
    return 0;
}