blob: 84b2438de95b603005eda7e200aa9fc13f490ab2 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
|
s=%s
// C=s^4+10*s^3+36*s^2+70*s+75 characteristic equation//
F=(s-2)^4+10*(s-2)^3+36*(s-2)^2+70*(s-2)+75 //shifting the origin with respect to s=-2//
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("there are %d roots to the right of s=-2 line",c+1)
else printf("system is stable")
end
|