1 2 3 4 5 6
function y=nqst1(z) s=%s H=syslin('c',(s+2)/((s+1)*(s-1))) nyquist(H) show_margins(H,'nyquist') endfunction