code: def test_psi_nlm():
      r = S("r")
      phi = S("phi")
      theta = S("theta")
      assert Psi_nlm(1, 0, 0, r, phi, theta) == exp(-r) / sqrt(pi)
      assert (Psi_nlm(2, 1, -1, r, phi, theta)) == S.Half * exp(
          -r / (2)
      ) * r * (sin(theta) * exp(-I * phi) / (4 * sqrt(pi)))
      assert (
          Psi_nlm(3, 2, 1, r, phi, theta, 2)
          == -sqrt(2)
          * sin(theta)
          * exp(I * phi)
          * cos(theta)
          / (4 * sqrt(pi))
          * S(2)
          / 81
          * sqrt(2 * 2**3)
          * exp(-2 * r / (3))
          * (r * 2) ** 2
      )
  
file_path: sympy/sympy/physics/tests/test_hydrogen.py
imports: from sympy.core.numbers import I, Rational, oo, pi
  from sympy.core.singleton import S
  from sympy.core.symbol import symbols
  from sympy.functions.elementary.exponential import exp
  from sympy.functions.elementary.miscellaneous import sqrt
  from sympy.functions.elementary.trigonometric import cos, sin
  from sympy.integrals.integrals import integrate
  from sympy.simplify.simplify import simplify
  from sympy.physics.hydrogen import R_nl, E_nl, E_nl_dirac, Psi_nlm
  from sympy.testing.pytest import raises
label_desc: Use image node_id 4 for calling a global function with example usage: test_psi_nlm() without return types