for k = 1:101 x(k) = (k-1)*0.01; y(k) = (k-1)*0.01;endfor k = 1:101 for m = 1:101 u(k,m) = sin(2*pi*x(m))*sin(2*pi*y(k))*exp(-(x(m)^2+y(k)^2)); endendpcolor(x,y,u)shading interphold onfor m = 1:9 pm = m*0.1; nm = m*(-0.1); contour(x,y,u,pm,'k-') contour(x,y,u,nm,'k--')endhold off

for ir = 1:21 for ith = 1:61 r = (ir-1)*0.05; th = (ith-1)*(pi/180); x(ir,ith) = r*cos(th); y(ir,ith) = r*sin(th); c(ir,ith) = r*th; endendpcolor(x,y,c)shading interphold oncontour(x,y,c,'k-')colorbaraxis equal tight

for k = 1:101 x(k) = (k-1)*0.01; y(k) = (k-1)*0.01;endfor p = 1:101 for q = 1:101 u(p,q) = 0; for n = 1:20 if (mod(n,2) == 0) an = 0; elseif (mod(n,4) == 1) an = 2/(n*n*pi*pi); else an = -2/(n*n*pi*pi); end an = (an*2)/sinh(n*pi); u(p,q) = u(p,q)+an*sinh(n*pi*(1-x(q)))*sin(n*pi*y(p)); end endendpcolor(x,y,u)shading interphold onfor m = 1:9 pm = m*0.05; contour(x,y,u,pm,'k-')endhold off

for ir = 1:21 for ith = 1:61 r = 0.5+(ir-1)*0.025; th = (ith-1)*(pi/30); x(ir,ith) = r*cos(th); y(ir,ith) = r*sin(th); c(ir,ith) = r*cos(3*th); endendpcolor(x,y,c)shading interphold oncontour(x,y,c,'k-')colorbaraxis equal tight