function ps10_q5_anim

     nx = 100;
     nt = 200;
     tend = 10;
     m = 100;

     x = [0:1.0/nx:1];

     for t=0:tend/nt:tend
     S = zeros(1,length(x));
     for j=1:m

        S = S - 4.0/pi^2 * (-1)^j/(2*j-1)^2 * cos((2*j-1)*pi*t) * sin((2*j-1)*pi*x);

     end
     plot(x,S);
     axis([0,1,-0.5,0.5]);
     pause(0.1);
     end

