Обратная матрица

public Matrix Inv
{
  get
  {
    if (!IsSquare)
      throw new MatrixMustBeSquare("Матрица должна быть квадратной");

    if (this.Det == 0)
      throw new Exception("Матрица вырожденная. Ошибка Inv");

    Matrix res = new Matrix(this);
               
    res = res.Join(Eye(this.rows), RowOrColumn.Row);

    for (int j = 0; j < res.rows; j++)
    {
      if (res[j, j] == 0)
      {
        Matrix tm = res.GetRow(j);
        for (int i = j + 1; i < res.rows; i++)
          if (res[i, j] != 0)
          {
            res = res.SetRow(j, res.GetRow(i));
            res = res.SetRow(i, tm);
            break;
          }
      }

      for (int i = j + 1; i < res.rows; i++)
      {
        Matrix tm1 = res.GetRow(i) - res.GetRow(j) * res[i, j] / res[j, j];
        res = res.SetRow(i, tm1);
      }
    }

    for (int j = res.rows - 1; j >= 0; j--)
    {
      if (res[j, j] == 0)
      {
        Matrix tm = res.GetRow(j);
        for (int i = j + 1; i < res.rows; i++)
          if (res[i, j] != 0)
          {
            res = res.SetRow(j, res.GetRow(i));
            res = res.SetRow(i, tm);
            break;
          }
      }
      for (int i = j - 1; i >= 0; i--)
      {
        Matrix tm1 = res.GetRow(i) - res.GetRow(j) * res[i, j] / res[j, j];
        res = res.SetRow(i, tm1);
      }

      Matrix tm2 = res.GetRow(j) / res[j, j];
      res = res.SetRow(j, tm2);
    }

     return res.GetSubMatrix(0, res.rows - 1, res.columns / 2, res.columns - 1);
  }
}