x0=0; xt=30; dt=.1;
L=(xt-x0)/dt+1;
n=10; 


figure(1)
for k=1:6
      t=5*k; 
      for i=1:L 
           summ=0;
           x(i)=x0+(i-1)*dt;
           for j=1:n 
                npi=j*pi; 
                a=2*(1-cos(npi))-npi^2*(1+cos(npi)); 
                b= exp((-npi^2*t)/900);
                c=sin(npi*x(i)/30);
                summ=summ+((a-b)/j^3)*b*c;
           end
           u(i)=30-x(i)+(60/pi^3)*summ;
      end
      plot(x,u);
      title(sprintf('u(x,%d)',t));
      m(k)=getframe();
end
movie(m)
movie2avi(m,'movieout.avi','compression','none');