Plotting the output of Part (c), Example 4.1.5, from Orfanidis

% Plot of y(n) = 4 u(n) - 4 u(n-25) - 3 (3/4)^n u(n) + 3 (3/4)^(n-25) u(n-25)

n=0:50;

un = [n >= 0];
unm25 = [(n-25) >= 0];
y = 4*un - 4*unm25 - 3*(3/4).^n.*un + 3*(3/4).^(n-25).*unm25;
stem(n,y);
title('y(n) = 4 u(n) - 4 u(n-25) - 3 (3/4)^n u(n) + 3 (3/4)^(n-25) u(n-25)')
xlabel('n')