def matrix_element(self, j, m, jp, mp):
        result = hbar*sqrt(j*(j + S.One) - mp*(mp + S.One))
        result *= KroneckerDelta(m, mp + 1)
        result *= KroneckerDelta(j, jp)
        return result
dataset/ETHPy150Open sympy/sympy/sympy/physics/quantum/spin.py/JplusOp.matrix_element/original