Gaussian Time Development

Constants and definitions

The constants are defined below.

tau=1.0;
hbar=1.0;
m=1.0;
k0=0.6;
a=0.1;

Generating the movie

The command below uses the definitions above to generate the frames of the movie.

tmax=6.0;
step=(tmax/50);
xmin=-0.5;
xmax=3.0;
Animate[Plot[((Exp[-((((x - ((hbar*k0*time)/m)))^2)/(2*(a^2)*(1 + ((time^2)/(tau^2)))))])/(Sqrt[2*π*(a^2)*(1+((time^2)/(tau^2)))])),{x,-1,5},
PlotRange->{{xmin,xmax},{0,5}},
Frame->True,
FrameLabel->{"x","P  "}],
{time,0,tmax,step}]