> with(plots);
> u:=(x,t)->(40/Pi^2)*evalf(sum((1/n^2)*sin(n*Pi/2)*sin(n*Pi/10)*sin(n*Pi*x/10)*sin(n*Pi*t/10),n=1..100));
> animate(u(x,t),x=0..10,t=0..20,frames=75,color=red);
>