% Example 3.19
% Feb 1, 2020

numb = [ 0 0 100  ];
denb = [ 1 10.1 101 ];
[ z, p, k ] = tf2zp( numb, denb )

%z = [  ]
%p = [ -5.0500+8.6889j	-5.0500-8.6889j ]'
%K = 100

s = tf( 's' );
sysb = 100*s/(s^3 + 10.1*s^2 + 101*s );
t = 0:0.01:2;

figure(1)
y = step( sysb, t );
plot( t, y )
