%
% Example 5.20 and Figure 5.6
%
num=[1 1]
den=[1 -1/6 -1/6]
h=dimpulse(num,den)
d=length(h)
k=0:1:d-1
plot(k,h,'*')
grid
xlabel('discrete time k')
ylabel('system impulse response')
print -deps fig5_6.eps