function ps10_q5_sol

     nx = 50;
     nt = 50;
     tend = 1.2;
     m = 1000;
     
     x = [0:1.0/nx:1];
     t = [0:tend/nt:tend];

     S = zeros(length(t),1) * 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

     surf(x,t,S);

