257/CH8/EX8.15/example_8_15.sce


1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

s=%s
F=2*s^5+s^4+6*s^3+s+1

disp(routh_t(F)) 
r=coeff(F)
routh=routh_t(F)
n=length(r)
c=0;
for i=1:n
if (routh(i,1)<0)
c=c+1;
 end 
  end
  if(c>=1)
    printf("system is unstable")
  else printf("there are no roots on RHS")
  end