n=-10:0.01:10;
om=2*pi/10;
mw=pi^1/4*exp(i*om*n).*exp(-(n.^2)/2);
r=abs(mw
plot(n,mw)