% Example 3.26
% Feb 3, 2020

s = tf( 's' )

sysH = ( 2*s + 1 )/( s^2+2*s+5 );

t = 0:0.1:6;
y = impulse( sysH, t );

plot( t, y, 'b','LineWidth',2 )
grid
hold on;

plot( t, 2*exp(-t), ':r','LineWidth',2 )
plot( t, -2*exp(-t), ':r','LineWidth',2 )

