
x = [-pi:0.005:pi];
tt = [0.01:0.1:3];

for i=1:length(tt)
A = (1/sqrt(0.1*tt(i)*2*pi))*exp(-x.*x./(4*0.1*tt(i)));
plot(x,A)
axis([-pi,pi,0,4])
M(:,i) = getframe;
end
