x = 1:20;  y1 = x;  y2 = sin(x/3); y3 = x.^2;     
subplot(3, 1, 1)
plot(x, y1, 'p') 
xlabel('x'), ylabel('y1'), title('just some arbitrary plot')
subplot(3, 1, 2) 
