s=%s
GH=syslin('c',40/((s+4)*(s^2+2*s+2)))
nyquist(GH)