function ps10_q2_sol

     nx = 50;
     ny = 50;
     m = 200;

     x = [0:pi/nx:pi];
     y = [0:pi/ny:pi];

     S = zeros(length(y),1) * x;
     for j=1:m

        S = S - 2.0*(-1)^j/(j*sinh(j*pi)) * ...
            (sinh(j*(pi-y))'*sin(j*x) + sin(j*y)'*sinh(j*(pi-x)));

     end

     surf(x,y,S);

