% U50, Different PD Control
% 

s = tf( 's' )

sysG1 = (1)/(s^2);
sysG2 = (s+1)/(s^2);
sysG3 = (s+1)/(s^2*(s+4));
sysG4 = (s+1)/(s^2*(s+9));
sysG5 = (s+1)/(s^2*(s+12));

rlocus( sysG1, 'r' );
hold on;
rlocus( sysG2, 'm' );
rlocus( sysG3, 'b' );
rlocus( sysG4, 'c' );
rlocus( sysG5, 'g' );


axis('equal')

axis([-6, 2, -3, 3]) 

% axis([-13, 1, -5, 5]) 

set(gco,'LineWidth',4)
% For example