1499/CH3/EX3.8/q8.sce


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

s=%s;
r=%r;

p=s^5+s^4+3*s^3+9*s^2+16*s+10
r=routh_t(p)
m=coeff(p)
l=length(m)
c=0;
for i=1:l
if (r(i,1)<0)
c=c+1;
end    
end
if(c>=1)
printf("System is unstable")
else("Sysem is stable")
end