Привести систему к треугольному виду

#include <stdio.h>
#include <iomanip>
#include <stdlib.h>
#include <time.h>
#include <math.h>
#include <iostream>
 
using namespace std;
 
int **matrix1(int &ROW , int &COL)
{
    cout<<"Введите размер массива: ";cin>>ROW>>COL;"\n";
    int **matr1=new int *[ROW];
    for (int i=0; i<ROW; i++)
        matr1[i]=new int [COL];
    cout<<"Введите элементы массива :"<<endl;
    for ( int i = 0; i < ROW; i++ )
        for ( int j = 0; j < COL; j++ ) 
             cin >> matr1[i][j];
    cout<<"\tИсходная матрица: ""\n";        
    for ( int i = 0; i < ROW; i++ )
    {
        for ( int j = 0; j < COL; j++ ) 
            cout << setw(4) << matr1[i][j] << " ";
    }
    cout << endl;
    return matr1;
}
 
int **matrix2(int &ROW , int &COL)
{
    int i,j;
    cout<<"Введите размер массива: ";cin>>ROW>>COL;"\n";
    int **matr2=new int *[ROW];
    for (int i=0; i<ROW; i++)
        matr2[i]=new int [COL];
    cout<<"Введите элементы массива :"<<endl;
    for ( int i = 0; i < ROW; i++ )
        for ( int j = 0; j < COL; j++ ) 
             cin >> matr2[i][j]; 
    cout<<"\tИсходная матрица: ""\n";        
    for ( int i = 0; i < ROW; i++ )
    {
        for ( int j = 0; j < COL; j++ ) 
            cout << setw(4) << matr2[i][j] << " ";
    }
    cout << endl;
    return matr2;
}
 
void exercise_1 (int **matr1,int ROW,int COL)
{
    float r;
    for(int k=0;k<ROW;k++)
    {
        for(int j=k+1;j<ROW+2;j++)
        {
            r=matr1[j][k]/matr1[k][k];
            for(int i=k;i<ROW+1;i++)
            {
                matr1[j][i]=matr1[j][i]-r*matr1[k][i];
            }
        }
    }
    cout<<"\tТреугольный вид: ""\n";
    for(int i=0;i<ROW;i++)                       
    {
        for(int j=0;j<COL;j++)
            cout<<setw(4)<<matr1[i][j]<<" ";  
        cout<<"\n";
    }
}
 
void exercise_2 (int **matr2,int ROW,int COL)
{
    int q=0;
    float num,sred_arif;
    cout<<"Введите элементы массива:"<<endl;
    for ( int i = 0; i < ROW; i++ )
        for ( int j = 0; j < COL; j++ ) 
            cin >> matr2[i][j]; 
    cout<<"\tВторая матрица: ""\n";        
        for ( int i = 0; i < ROW; i++ )
        {
            for ( int j = 0; j < COL; j++ ) 
                cout << setw(4) << matr2[i][j] << " "; 
            cout << endl;
        }
    cout<<"Введите контрольное значение: ";cin>>num;
    for (int i=0; i<ROW; i++)
    {
        sred_arif=0;
        for (int j=0; j<COL; j++)
            sred_arif+= matr2[i][j];
        sred_arif/=ROW;
        if (sred_arif<num)
            q++;
    }
    cout<<"Количество строк: "<<q<<endl;
    cout<<"Среднее арифметическое элементов которых меньше заданного равно: "<<sred_arif;
}
 
int main()
{   
    int ROW,COL,**A,**B;
    A=matrix1(ROW,COL);
    exercise_1(A,ROW,COL);
    B=matrix2(ROW,COL);
    exercise_1(B,ROW,COL);
    return 0;
}

for(int j=k+1;j<ROW+2;j++) // <-- j в один прекрасный момент станет больше ROW…
{
    r=matr1[j][k]/matr1[k][k];
    for(int i=k;i<ROW+1;i++) // <-- здесь то же самое…
    {
        matr1[j][i]=matr1[j][i]-r*matr1[k][i];