1 2 3 4 5 6 7 8
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)