code: def test_sho_R_nl():
      omega, r = symbols("omega r")
      l = symbols("l", integer=True)
      u = Function("u")
  
      # check that it obeys the Schrodinger equation
      for n in range(5):
          schreq = -diff(u(r), r, 2) / 2 + (
              (l * (l + 1)) / (2 * r**2)
              + omega**2 * r**2 / 2
              - E_nl(n, l, omega)
          ) * u(r)
          result = schreq.subs(u(r), r * R_nl(n, l, omega / 2, r))
          assert simplify(result.doit()) == 0
  
imports: from sympy.core import symbols, Rational, Function, diff
  from sympy.physics.sho import R_nl, E_nl
  from sympy.simplify.simplify import simplify
label_desc: Use image node_id 1 for calling a global function with example usage: test_sho_R_nl() without return types