Задача: привести систему к треугольному виду

[limaxY]

Условие задачи:
Коэфиценты системы линейных уравнений заданы в виде прямоугольной матрицы. С помощью допустимых
преобразований привести систему к треугольному виду.

Найти количество строк, среднее арифметическое элементов которых меньше заданной величины.

Задача решена, но не могу привести никак к треугольному виду, помогите пожалуйста:

#include <iostream>
#include <stdio.h>
using namespace std; 
 
int main() // головной модуль
{
        setlocale (LC_ALL,"Rus");
        int a[10][10]; 
        int m, n;
        int i, j; 
        int b; 
        cout << "Введите количество строк: ";
        cin >> m;
        cout << "Введите количество столбцов: ";
        cin >> n;
        cout << "Введите элементы: \n";
        
        for(i = 0; i < m; i++)
                for(j = 0; j < n; j++)
                        
        cin >> a[i][j];
        
        cout << "\nМатрица: ";  
        for(i = 0; i < m; i++)
        {
                cout << "\n";           
                for(j = 0; j < n; j++)
                        {
                                cout << a[i][j];
                                cout << " \t"; 
                        }
        };
                
                cout << "\n\nВведите величину: ";
                cin >> b;
                int s = 0, kol = 0;
                
                for (i = 0; i < m; i++)
                        {
                                s = 0;
                                for(j = 0; j < n; j++)
                                        {
                                                s += a[i][j];
                        
                                        }
                                s = s / n;
                                if(s < b) kol++;
                        }
                cout << "\nКоличество строк: " << kol << " \n";
}

[KobolD]

тебе можно воспользоваться прямым ходом из метода гауса.
матрица коэфициентов a[][], вектор свободных членов b[], колличество строк n.

        // Прямой ход метода 
        double M;
        for(int i = 0; i<n-1; i++)
        {
                for(int j = i; j<n-1; j++)
                {
                        M = -a[j+1][i]/a[i][i];
                        for(int k = i; k<n; k++)
                        {
                                a[j+1][k] += a[i][k]*M; 
                                if(a[j+1][k] < 0.000001 & a[j+1][k]> -0.000001)
                                        a[j+1][k] = 0;
                        }
                        b[j+1] += b[i]*M;
                }
        }