x= [-5:.01:5];


init= exp(-(x.^2));

hold on
plot(x,init,'g')

for t = .02:1:10.02
    
    u = exp(-(x.^2)./(1 +4*t))./ sqrt(1 + 4*t);
    plot(x,u)
    
end

hold off