n = 0:20;
step = ones(size(n)); 
ramp = n; 
parabolic = n.^2; 
sinusoidal = sin(0.2 * pi * n); 

figure;
subplot(2,2,1); 
stem(n, step, 'r'); 
title('Step Function'); 
xlabel('n'); ylabel('Amplitude');

subplot(2,2,2); 
stem(n, ramp, 'g'); 
title('Ramp Function'); 
xlabel('n'); ylabel('Amplitude');

subplot(2,2,3); 
stem(n, parabolic, 'b'); 
title('Parabolic Function'); 
xlabel('n'); ylabel('Amplitude');

subplot(2,2,4); 
stem(n, sinusoidal, 'm'); 
title('Sinusoidal Function'); 
xlabel('n'); ylabel('Amplitude');