summaryrefslogtreecommitdiff
path: root/1151/CH6/EX6.10/example10.sce
blob: 69ada3df07f987891dd7d811a265231270198465 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
s=%s;
m=s ^5+2*s ^4+24*s^3+48*s^2-25*s-50;
disp (m)
r= coeff (m)
n= length (r)
routh = routh_t (m)
disp (routh ," r o u t h s t a b u l a t i o n s=")
c =0;
for i =1: n
 if ( routh (i ,1) <0)
c=c+1;
 end
 end
if(c >=1)
printf (" sys t em i s u n s t a b l e ")
else printf (" sys t em i s ma r g i n a l l y s t a b l e ")
 end