1 2 3 4 5 6 7 8
s=%s sys=syslin('c',(k)/(s*(s+1)*(s+4))) evans(sys) //breakaway points disp("break away points are") disp(roots(numer(derivat(sys))))