blob: f3b367aeafa538415df5f79bb5ea5df2977591b4 (
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
F=s^6+2*s^5+5*s^4+8*s^3+8*s^2+8*s+4
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
disp("s^2 is")
k=(roots(routh(3,:)))
disp(k)
end
if(k(2,1)==k(3,1))
printf("repeated roots on imaginary axis. hence unstable system")
|
