In [1]:

    

%pylab
from qutip import *


    

    




Using matplotlib backend: MacOSX
Populating the interactive namespace from numpy and matplotlib

In [74]:

    

N=100
qmin = -60
qmax = 5
pmin = -60
pmax = 30
qvec = linspace(qmin,qmax,200)
pvec = linspace(pmin,pmax,200)
qdistvac = qfunc(coherent_dm(N,0),qvec,pvec)
qdist1 = qfunc(coherent_dm(N,sqrt(10)*exp(1j*0.85*pi)),qvec,pvec)
qdist2 = qfunc(coherent_dm(N,sqrt(60)*exp(1j*0.75*pi)),qvec,pvec)


    

In [76]:

    

imshow(qdistvac+qdist1+qdist2,extent=(qmin,qmax,pmin,pmax),origin='lower')


    

    
Out[76]:





<matplotlib.image.AxesImage at 0x1122e0d50>

In [ ]: