blob: 2eb74881b263f93e59e8fbe3f7c4a56e842d7ec9 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
|
s=%s;
//F=s^5+s^4+2*s^3+2*s^2+3*s+15//
//replacing 's' by '1/z' we get F//
z=%z;
F=15*z^5+3*z^4+2*z^3+2*z^2+z+1;
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 printf("system is stable")
end
|
