summaryrefslogtreecommitdiff
path: root/162/CH8/EX8.3/example83.sce
blob: bcc08506fc0716c947bcdb039b25a13d934b47ad (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
//Example 8.3
clc;
// Define the polynomial
s=poly(0,"s");
p=3+10*s+5*s^2+s^3;
// Calculate the routh of above polynomial
r=routh_t(p);
disp(r,"Routh array=");
A=r(:,1);
c=0;
x1=0;
eps=0;
for i=1:4
   x1=A(i,1);
   if x1<0
     c=c+1;
   end
end
   if(c>=1) then
       printf("system is unstable");
   else 
       printf("system is stable");
   end
x=roots(p);