s=%s
sys=syslin('c',(k)*(s+4)/((s+2)*(s^2+2*s+2)))
evans(sys)