blob: 36729e0ff166383af73ea88b3b0cb98f0b8dd423 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
|
//terms in a row become infinite//
s=%s;
//P=s^5+2*s^4+3*s^3+6*s^2+2*s+1;//
//replace 's' by '1/z' to get F//
z=%z;
F=z^5+2*z^4+6*z^3+3*z^2+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
|
