gauss, który wydaje się działać

#include <iostream>
#include <vector>


using namespace std;

vector<double> mnozenie(vector<double> matylda, double liczba)
{
    for(int i=0; i<matylda.size(); i++)
    {
        matylda[i] = matylda[i]*liczba;
    }
    return matylda;
}

//matylda - wladyslaw
vector<double> odejmowanie(vector<double> matylda, vector<double> wladyslaw)
{
    for(int i=0; i<matylda.size(); i++)
    {
        matylda[i] = matylda[i]-wladyslaw[i];
    }
    return matylda;
}

int main()
{
    int t;
    cin >> t;

    for(int iterator=0; iterator<t; iterator++)
    {
        //tu wstaw operacje na macierzach
        int liczba_rownan;
        cin >> liczba_rownan;
        vector<vector<double> > macierz;

        for(int i=0; i<liczba_rownan; i++)
        {
            vector<double> usagi;
            for(int j=0; j<liczba_rownan+1; j++)
            {
                double sailor_moon;
                cin >> sailor_moon;
                usagi.push_back(sailor_moon);
            }
            macierz.push_back(usagi);
        }
        //jebac wyznacznik glowny
        //eliminacja gaussa bez cudow
        vector<vector<double> > eg = macierz;
        vector<double> x_eg;
        for(int i=0; i<eg.size(); i++)
        {
            x_eg.push_back(0);
        }

        for(int i=0; i<eg.size(); i++)
        {
            //szukamy najwiekszego w kolumnie i
            //tylko ostatniej kolumny nie ruszamy
            int najnaj=i;
            for(int j=i; j<eg.size(); j++)
            {
                if(eg[j][i] > eg[najnaj][i])
                    najnaj=j;
            }
            eg[i].swap(eg[najnaj]);
            //jej, swapniete

            //teraz jak juz wybralismy najwiekszy element
            //i jest on w eg[i][i] to ten wiersz dzielimy przez ten element
            //robiac w miejscu eg[i][i] jedynke
            eg[i] = mnozenie(eg[i], 1/eg[i][i]);

            //for(int j=i; j<eg.size(); j++)
            //{
            //  //if(j==najnaj)
            //  //  eg[j][i]=1;
            //  //else
            //  eg[j][i] = eg[j][i]/eg[i][i];
            //  //eg[j][i] -= eg[j][i]/eg[najnaj][i];
            //}
            //pieknie, mam jedynki. co teraz z nimi?
            // a no trzeba wiersze odejmowac. a to juz nie jest takie proste.

            for(int j=i+1; j<eg.size(); j++)
            {
                //matylda-wladyslaw;
                eg[j]=odejmowanie(eg[j], mnozenie(eg[i], eg[j][i]));
            }
            //super, mamy macierz trojkatna
        }

        cout << "\n";
        for(int i=0; i<eg.size(); i++)
        {
            for(int j=0; j<eg.size() + 1; j++)
            {
                cout << eg[i][j] << " ";
            }
            cout << "\n";
        }

        x_eg.back() = eg.back().back();
        //ostatniego x-a znamy juz bez kombinowania
        for(int i=eg.size()-2; i>=0; i--)
        {
            //;
            x_eg[i] = eg[i][eg.size()];
            //dajemy mu wartosc b, bedziemy od niej odejmowac
            for(int j=eg.size()-1; j>i; j--)
            {
                x_eg[i] -= x_eg[j] * eg[i][j];
            }
        }

        cout << "\n";
        for(int i=0; i<eg.size(); i++)
        {
            for(int j=0; j<eg.size() + 1; j++)
            {
                cout << eg[i][j] << " ";
            }
            cout << "\n";
        }
        cout << "\n";
        for(int i=0; i<x_eg.size(); i++)
            cout << x_eg[i] << " ";

        cout << "\n";
        //macierz.clear();
    }


    return 0;
}