% % Example 4.20 and Figure 4.2 % num=[1 3] den=[1 3 2] figure (1) impulse(num,den) print -deps fig4_2.eps t=0:0.1:3 y=impulse(num,den,t) figure (2) plot(t,y) xlabel('time') ylabel('impulse response')