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