def _diagonal_solve(self, rhs):
        return self._new(rhs.rows, rhs.cols, lambda i, j: rhs[i, j] / self[i, i])
Helper function of function diagonal_solve,
dataset/ETHPy150Open sympy/sympy/sympy/matrices/dense.py/DenseMatrix._diagonal_solve