code: def test_wavefunction():
      Psi = {
          0: (nu / pi) ** Rational(1, 4) * exp(-nu * x**2 / 2),
          1: (nu / pi) ** Rational(1, 4)
          * sqrt(2 * nu)
          * x
          * exp(-nu * x**2 / 2),
          2: (nu / pi) ** Rational(1, 4)
          * (2 * nu * x**2 - 1)
          / sqrt(2)
          * exp(-nu * x**2 / 2),
          3: (nu / pi) ** Rational(1, 4)
          * sqrt(nu / 3)
          * (2 * nu * x**3 - 3 * x)
          * exp(-nu * x**2 / 2),
      }
      for n in Psi:
          assert simplify(psi_n(n, x, m, omega) - Psi[n]) == 0
  
file_path: sympy/sympy/physics/tests/test_qho_1d.py
imports: from sympy.core.numbers import Rational, oo, pi
  from sympy.core.singleton import S
  from sympy.core.symbol import Symbol
  from sympy.functions.elementary.exponential import exp
  from sympy.functions.elementary.miscellaneous import sqrt
  from sympy.integrals.integrals import integrate
  from sympy.simplify.simplify import simplify
  from sympy.abc import omega, m, x
  from sympy.physics.qho_1d import psi_n, E_n, coherent_state
  from sympy.physics.quantum.constants import hbar
label_desc: Use image node_id 1 for calling a global function with example usage: test_wavefunction() without return types