% Example for subplot
[x,y] = meshgrid(linspace(-pi,pi,100));
 
subplot(2,2,1);
title('Vibration modes of a rectangular membrane');
z = cos(x/2).*cos(y/2); 
surf(x,y,z,z)

subplot(2,2,2);
z = sin(x).*cos(y/2); 
surf(x,y,z,z)

subplot(2,2,3);
z = sin(y).*cos(x/2); 
surf(x,y,z,z)

subplot(2,2,4);
z = sin(x).*sin(y); 
surf(x,y,z,z)