> with(plots):

> animate(sum(40/((2*n-1)*Pi)*sin((2*n-1)*Pi*x)*
exp(-(2*n-1)^2*Pi^2*t),n=1..20),
x=0..1,t=0..0.6,color=green,thickness=3,frames=50);

>

>