Example 4.1.9 from Orfanidis

h = (0.95).^(0:24);
n=0:120;
x=[n==0] + 2*[n==40] + 2*[n==70] + [n==80];
y = conv(h,x);
stem(n,x)
hold on
plot(n,y(1:length(n)),'--');
axis([0 120 0 3])
xlabel('time samples')
legend('o','input','--','output');