Как через push_back заполнить вектор векторов?

typedef vector<vector<double>>Elements;//вектор векторов
class Matrix
{
    Elements M;
    int R, C;//строки столбцы
public:
    //по умолчанию незадаваемые эл-ты равны нулю
    Matrix(int r, int c);//пр¤моуг матрица
    Matrix(int n);//квадр м-ца
    Matrix(int n, double d);
    ~Matrix() {};
    int GetR();
    int GetC();
    void rand(int mod);//зап случ знач по модулю mod
 
    Matrix inverse();//обрат
 
    
    vector<double> & operator[](int r);//доступ к элементу
 
    friend ostream & operator<< (ostream &os, Matrix &A); //вывод м-цы
 
};
 
 
Matrix Matrix::inverse()
{
    Matrix Ones(R, 1.0);
    Matrix A = *this;
    for (int i = 0; i < R; i++) //по всем строкам
    {
        int maxi = i;
        for (int s = i + 1; s < R; s++)
            maxi = abs(A[maxi][i]) < abs(A[s][i]) ? s : maxi; //Если maxi = abs(A[maxi][i]) < abs(A[s][i]) истинно, результатом операции будет s, в противном случае maxi
        swap(A[i], A[maxi]); 
        swap(Ones[i], Ones[maxi]);
        if (abs(A[i][i]) < eps)
            throw string("Matrix is not invertable");
 
        double buf = A[i][i];
        for (int ii = 0; ii < R; ii++)//то по всей строке
        {
            A[i][ii] /= buf;//делаем 1
            Ones[i][ii] /= buf;
        }
        for (int j = 0; j < R; j++)
        {
            if (i == j) continue;
            double buf = A[j][i];
            for (int jk = 0; jk < R; jk++)//то по всей строке
            {
                A[j][jk] -= A[i][jk] * buf;//вычитаем строку
                Ones[j][jk] -= Ones[i][jk] * buf;
            }
        }
    }
    return Ones;
}

int cFoo(
    int m,          // Число строк матрицы
    int n,          // Число столбцов матрицы
    double *a,      // Адрес массива элементов матрицы
    double eps      // Точность вычислений
    ) {
    int i, j, k, l;
    double r;
 
    i = 0; j = 0;
    while (i < m && j < n) {
        r = 0.0;
        for (k = i; k < m; ++k) {
            if (fabs(a[k*n + j]) > r) {
                l = k;     
                r = fabs(a[k*n + j]); 
            }
        }
        
        // etc
    }
 
    return i; // Возвращаем число ненулевых строк
}