Матрица. Исключение

std::vector<std::vector<double>> Matrix::inverse(std::vector<std::vector<double>> matrix)
{
    std::vector <std:: vector <double>> matrix1(a, std::vector<double>(2*b));
    double r;
    double g;
    for (int i = 0; i < a; i++)
    {
        for (int j = 0; j <2* b; j++)
        {
            matrix1[i][j]=matrix[i][j];
        }
    }
    for (int i = 0; i < a; i++)
    {
        for (int j = 0; j < 2 * b; j++)
        {
            if (i == (j - a))
                matrix1[i][j] = 1;
            else
                matrix1[i][j] = 0;
        }
}
    for (int i = 0; i < a; i++)
    {
        for (int j = 0; j < b; j++)
        {
            if (i != j)
            {
                r = matrix1[j][i] / matrix1[i][i];
                for (int k = 0; k < 2 * b; k++)
                {
                    matrix1[j][k] =- r * matrix1[i][k];
                }
            }
        }
    }
    for (int i = 0; i < a; i++)
    {
        g = matrix1[i][i];
        for (int j = 0; j < 2 * b; j++)
        {
            matrix1[i][j] /= g;
        }
    }
    for (int i = 0; i < a; i++)
    {
        for (int j = b; j < 2 * b; j++)
        {
            std::cout << std::fixed << std::setprecision(4) << matrix1[i][j] << "\t";
            //std::cout << "\n";
        }
        std::cout << "\n";
    }
    return matrix1;
}

std::vector<std::vector<double>> Matrix::razmer(std::vector < std::vector <double>> matrix) {
    a = matrix.size();
    b = matrix.size();
    return matrix;
}

for (int i = 0; i < a; i++)
    matrix1[i][i] += 1;

for (int i = 0; i < a; i++)
    {
        for (int j = b; j < 2 * b; j++) //  j=b  !!!!!!!!!!!!!!!!
        {
            if (i == (j - a))
                matrix1[i][j] = 1;
            else
                matrix1[i][j] = 0;
        }
}

std::vector<std::vector<double>> Matrix::inverse(std::vector<std::vector<double>> matrix)
{
    std::vector <std::vector <double>> matrix1(a, std::vector<double>(2 * b));
    float temp;
    for (int i = 0; i < a; i++)
    {
        for (int j = 0; j < b; j++)
        {
            matrix1[i][j] = matrix[i][j];
        }
    }
    for (int i = 0; i < a; i++)
    {
        for (int j = 0; j < 2 * b; j++)
        {
            if ((i+a) == j)
                matrix1[i][j] = 1;
        }
    }
    for (int i = a - 1; i > 0; i--) {
        if (matrix1[i - 1][0] < matrix1[i][0]) {
            for (int j = 0; j < 2 * b; j++) {
                temp = matrix1[i][j];
                matrix1[i][j] = matrix1[i - 1][j];
                matrix1[i - 1][j] = temp;
            }
        }
    }
    for (int i = 0; i < a; i++) {
        for (int j = 0; j < b; j++) {
            if (j != i) {
                temp = matrix1[j][i] / matrix1[i][i];
                for (int k = 0; k < 2 * b; k++) {
                    matrix1[j][k] -= matrix1[i][k] * temp;  }
            }
        }
}
    for (int i = 0; i < a; i++) {
        temp = matrix1[i][i];
 
        for (int j = 0; j < 2* b; j++) {
            matrix1[i][j] = matrix1[i][j]/temp;
        }
}
    for (int i = 0; i < a; i++) {
        for (int j = 0; j <  2*b; j++) {    //[COLOR="Red"]ЕСЛИ ЗДЕСЬ ОСТАВИТЬ просто b , то выводится только единичная [/COLOR]
            std::cout << std::fixed << std::setprecision(4) << matrix1[i][j] << "\t";
        }
        std::cout << std::endl;
 
    }
    return matrix1;

for (int j = 0; j < b; j++) { // для вывода левой половины

for (int j = b; j < 2*b; j++) { // для вывода правой половины