open-loop transfer function:
G(s)=5*(0.0167s+1)/[s*(0.03s+1)(0.0025s+1)(0.001s+1)]
bode picture in matlab
>Solution :
num=5*[0.0167,1];
den=conv(conv([1,0],[0.03,1]),conv([0.0025,1],[0.001,1]));
G=tf(num,den);
w=logsapce(0,4,50);
bode(G,w);grid;
[Gm,Pm,Wcg,Wcp]=margin(G)