double calcDet(vector<vector<double>> &Matrix) {
        double det = 0;
        if (Matrix.size() == 1) {
            return Matrix[0][0];
        }
        else if (Matrix.size() == 2) {
            det = Matrix[0][0] * Matrix[1][1] - Matrix[0][1] * Matrix[1][0];
            return det;
        }
        else {
            for (int p = 0; p < Matrix[0].size(); p++) {
                vector<vector<double>> TempMatrix;
                for (int i = 1; i < Matrix.size(); i++) {
                    vector<double> TempRow;
                    for (int j = 0; j < Matrix[i].size(); j++) {
                        if (j != p) {
                            TempRow.push_back(Matrix[i][j]);
                        }
                    }
                    if (!TempRow.empty())
                        TempMatrix.push_back(TempRow);
                }
                det = det + Matrix[0][p] * pow(-1, p) * calcDet(TempMatrix);
            }
            return det;
        }
    }

    double Matrix::det() const {
        if (matSizeRow != matSizeCol) {
            throw("matrix is not square");
        }
        vector<vector<double>> matvec(matSizeRow, vector<double>(matSizeRow));
        for (int i = 0; i < matSizeRow; i++) {
            for (int j = 0; j < matSizeRow; ++j) {
                matvec[i][j] = get(i, j);
            }
        }
        return calcDet(matvec);
    }