code: def test_energy():
      n, l, hw = symbols("n l hw")
      assert (
          simplify(E_nl(n, l, hw) - (2 * n + l + Rational(3, 2)) * hw)
          == 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 2 for calling a global function with example usage: test_energy() without return types