function psi_x=moller_squeeze_coherent_wavefunction_hbar(qz,pz,zeta1,x,hbar)
format long
%zeta1=1.0+0.0i
%alpha1=1.0+1.0i;
%x=[-2:0.01:2];
%%%%%
disp('Moller Wavefunction Representation of the Squeezed Coherent_State')
r=abs(zeta1);
theta=angle(zeta1);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
ces_p=cosh(r)+exp(1i*theta)*sinh(r);
ces_n=cosh(r)-exp(1i*theta)*sinh(r);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
N1=(1/(sqrt(sqrt(pi*hbar))))*(1/sqrt(ces_p))
A1=(-1/(2*hbar))*(ces_n/ces_p)
Bx2=(x-qz).^2;
D=(+1i/hbar)*pz*(x-qz/2);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
psi_x=N1*exp(A1.*Bx2+D);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%figure,plot(x,(psi_x.*conj(psi_x)),'.')
%figure,plot(x,abs(fftshift(fft(psi_x))).^2)
%vecnorm=trapz(x,psi_x.*conj(psi_x))
%psi_x=psi_x./sqrt(vecnorm);
%%%%This wavefunction is seems to be OK !!!
end