#include <iostream>
#include <math.h>
#include <stdlib.h>
using namespace std;

int main()
{
    int i, j;
    // создаем массив

    const int n = 4, m = 5;
    float matrix[n][m] = {{-2, -9, -3, 7, -26}, {-7, 2, 2, 5, -25}, {-6, 2, 0, 0, -16}, {0, -3, 8, -3, -5}};

    //выводим массив
    cout << "matrix: " << endl;
    for (i = 0; i < n; i++)
    {
        for (j = 0; j < m; j++)
            cout << matrix[i][j] << " ";
        cout << endl;
    }
    cout << endl;

    //Метод Гаусса
    //Прямой ход, приведение к верхнетреугольному виду
    float tmp, xx[m];
    int k;

    for (i = 0; i < n; i++)
    {
        tmp = matrix[i][i];
        for (j = n; j >= i; j--)
            matrix[i][j] /= tmp;
        for (j = i + 1; j < n; j++)
        {
            tmp = matrix[j][i];
            for (k = n; k >= i; k--)
                matrix[j][k] -= tmp * matrix[i][k];
        }
    }
    /*обратный ход*/
    xx[n - 1] = matrix[n - 1][n];
    for (i = n - 2; i >= 0; i--)
    {
        xx[i] = matrix[i][n];
        for (j = i + 1; j < n; j++)
            xx[i] -= matrix[i][j] * xx[j];
    }

    //Выводим решения
    for (i = 0; i < n; i++)
        cout << xx[i] << " ";
    cout << endl;

    return 0;
}