1 2 3 4 5 6 7 8 9
function y=nqst10(z) s=%s H=syslin('c',1/(s*(0.2*s+1)*(0.05*s+1))) nyquist(H) show_margins(H,'nyquist') mtlb_axis([-1 1 -5 1]) gm=g_margin(H) pm=p_margin(H) endfunction