path: root/257/CH9/EX9.20/example_9_20.sce

s=%s
sys=syslin('c',(s+1)*k/(s*(s-1)*(s^2+5*s+20)))
evans(sys)

//stability

[Ki,s]=kpure(sys)
disp(Ki)