root locus -lead compensation

s = tf('s');
sys = (20*16.8)/(s*(s + 5)*(s + 20));
rlocus(sys)
axis([-11 0 -1 4])
Zeta = 0.5;
Wn = 4.8;
sgrid(Zeta,Wn)