//equation// s=%s; m=s^5+6*s^4+3*s^3+2*s^2+s+1 r=coeff(m) n=length(r) routh=[r([6,4,2]);r([5,3,1])] routh=[routh;-det(routh(1:2,1:2))/routh(2,1),-det(routh(1:2,2:3))/routh(2,2),0] routh=[routh;-det(routh(2:3,1:2))/routh(3,1),-det(routh(2:3,2:3))/routh(3,2),0] routh=[routh;-det(routh(3:4,1:2))/routh(4,1),-det(routh(3:4,2:3))/routh(4,2),0] routh=[routh;-det(routh(4:5,1:2))/routh(5,1),0,0] c=0; for i=1:n if (routh(i,1)<0) c=c+1; end end if(c>=1) printf("system is unstable") else ("system is stable") end