From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001
From: priyanka
Date: Wed, 24 Jun 2015 15:03:17 +0530
Subject: initial commit / add all books

---
 29/CH12/EX12.17/exa12_17.sce | 27 +++++++++++++++++++++++++++
 1 file changed, 27 insertions(+)
 create mode 100755 29/CH12/EX12.17/exa12_17.sce

(limited to '29/CH12/EX12.17/exa12_17.sce')

diff --git a/29/CH12/EX12.17/exa12_17.sce b/29/CH12/EX12.17/exa12_17.sce
new file mode 100755
index 000000000..7e9cd520d
--- /dev/null
+++ b/29/CH12/EX12.17/exa12_17.sce
@@ -0,0 +1,27 @@
+//caption:stability_using_Routh-hurwitz_criterion
+//example 12.17
+//page 530
+s=%s;
+syms  K
+G=sym('K/((s^3+3*s+2)*(s^2+6*s+24))');
+H=1;
+CH=((s^3+3*s+2)*(s^2+6*s+24)+K)
+disp('=0',CH,"characterstics_eq,CH=")
+c0=coeffs(CH,'s',0);
+c1=coeffs(CH,'s',1);
+c2=coeffs(CH,'s',2);
+c3=coeffs(CH,'s',3);
+c4=coeffs(CH,'s',4);
+b=[c0 c1 c2 c3 c4 ]
+routh=[b([5,3,1]);b([4,2]),0]
+routh=[routh;-det(routh(1:2,1:2))/routh(2,1),routh(1,3),0]
+routh(3,1)=simple(routh(3,1))
+t=routh(2:3,1:2)
+l=simple(-det(t)/t(2,1))
+routh=[routh;l,0,0]
+//routh=[routh;K,0,0]
+disp(routh,"routh=")
+disp("for given system to be stable,following condition should be satisfied");
+disp("78.84-0.259K>0")
+disp("which gives limiting value of K")
+disp("K<288.9");
-- 
cgit