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
|