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