def coherent_state(n, alpha):

    return exp(- Abs(alpha)**2/2)*(alpha**n)/sqrt(factorial(n))
Returns <n|alpha> for the coherent states of 1D harmonic oscillator.
dataset/ETHPy150Open sympy/sympy/sympy/physics/qho_1d.py/coherent_state