s=%s; sys1=syslin('c',k*(s+1)/((s-1)*(s))) nyquist(sys1)