def E_n(n, omega):
 return ((hbar * omega) * (n + Rational(1, 2)))

def E_n ( n , omega ) : return ( ( hbar * omega ) * ( n + Rational ( 1 , 2 ) ) )
returns the energy psi_{n} for a 1d potential hole with infinity borders n the "principal" quantum number .