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