summaryrefslogtreecommitdiff
path: root/257/CH8/EX8.25/example_8_25.sce
blob: 0d03de8b80b0466a1ec50478b827d857b909ac4c (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
s=%s
P=s^5-s^4-2*s^3+2*s^2-8*s+8
routh=routh_t(P)
disp(routh)
r=coeff(P)
n=length(r)
c=0;
for i=1:n
if (routh(i,1)<0)
c=c+1;
 end 
  end
  if(c>=1)
    printf("there are %d roots on RHS",c)
  else printf("there are no roots on RHS")
  end
  
  disp("s^2 is")
  disp(roots(routh(2,:)))