#include <iostream>
#include <iomanip>
#include "Matrix.h"

using namespace std;

int main() {
    while (true) {
        int m, n;
        cout << "Enter the matrix dimensions m, n: ";
        cin >> m >> n;
        if (!cin)
            break;
            
        cout << fixed << setprecision(4);    

        Matrix a(m, n);
        cout << "Enter matrix elements:" << endl;
        cin >> a;
        cout << "Initial matrix:" << endl;
        cout << a;
        
        Matrix b = a;
        int r = b.gaussElimination();
        cout << "Row echelon form:" << endl;
        cout << b;
        cout << "rank = " << r << endl;
        if (a.rows() == a.cols()) {
            cout << "determinant = " << a.det() << endl;
        }
        //... to do...

    }
    return 0;
}