B ˜‘[Óã@shddlmZmZddlmZmZmZddlmZm Z m Z m Z m Z ddl mZdd„Zdd„Zd d „Zd S) é)Úprint_functionÚdivision)ÚSÚpiÚRational)ÚhermiteÚsqrtÚexpÚ factorialÚAbs)ÚhbarcCs~tt||||gƒ\}}}}||t}|ttdƒdtdd|t|ƒƒ}|t| |ddƒt|t|ƒ|ƒS)aK Returns the wavefunction psi_{n} for the One-dimensional harmonic oscillator. ``n`` the "nodal" quantum number. Corresponds to the number of nodes in the wavefunction. n >= 0 ``x`` x coordinate ``m`` mass of the particle ``omega`` angular frequency of the oscillator Examples ======== >>> from sympy.physics.qho_1d import psi_n >>> from sympy import var >>> var("x m omega") (x, m, omega) >>> psi_n(0, x, m, omega) (m*omega)**(1/4)*exp(-m*omega*x**2/(2*hbar))/(hbar**(1/4)*pi**(1/4)) ééé)Úmaprr rrr r r)ÚnÚxÚmÚomegaZnuÚC©rú3lib/python3.7/site-packages/sympy/physics/qho_1d.pyÚpsi_ns ,rcCst||tddƒS)aú Returns the Energy of the One-dimensional harmonic oscillator ``n`` the "nodal" quantum number ``omega`` the harmonic oscillator angular frequency The unit of the returned value matches the unit of hw, since the energy is calculated as: E_n = hbar * omega*(n + 1/2) Examples ======== >>> from sympy.physics.qho_1d import E_n >>> from sympy import var >>> var("x omega") (x, omega) >>> E_n(x, omega) hbar*omega*(x + 1/2) r r)r r)rrrrrÚE_n+srcCs*tt|ƒd dƒ||tt|ƒƒS)zð Returns for the coherent states of 1D harmonic oscillator. See http://en.wikipedia.org/wiki/Coherent_states ``n`` the "nodal" quantum number ``alpha`` the eigen value of annihilation operator r)r r rr )rZalpharrrÚcoherent_stateGs rN)Z __future__rrZ sympy.corerrrZsympy.functionsrrr r r Zsympy.physics.quantum.constantsr rrrrrrrÚs  #