From 7bc77cb1ed33745c720952c92b3b2747c5cbf2df Mon Sep 17 00:00:00 2001
From: prashantsinalkar
Date: Sat, 3 Feb 2018 11:01:52 +0530
Subject: Added new code

---
 3885/CH5/EX5.8/Ex5_8.sci | 28 ++++++++++++++++++++++++++++
 1 file changed, 28 insertions(+)
 create mode 100644 3885/CH5/EX5.8/Ex5_8.sci

(limited to '3885/CH5/EX5.8')

diff --git a/3885/CH5/EX5.8/Ex5_8.sci b/3885/CH5/EX5.8/Ex5_8.sci
new file mode 100644
index 000000000..9de51e2a7
--- /dev/null
+++ b/3885/CH5/EX5.8/Ex5_8.sci
@@ -0,0 +1,28 @@
+//control systems by Nagoor Kani A
+//Edition 3
+//Year of publication 2015
+//Scilab version 6.0.0
+//operating systems windows 10
+// Example 5.8
+
+clc;
+clear; 
+s=poly(0,'s')
+a=(s^6)+(s^5)+(3*s^4)+(3*s^3)+(3*s^2)+(2*s)+1
+b=coeff(a)
+n=length(b)
+R=routh_t(a)
+disp(R,'the routh array is;') 
+//characterstic polynomial can be expressed as product of auxillary polynomial and quotient polynomial
+//divide characterstic equation by auxilary polynomial to get quotient polynomial
+//routh table for quotient poly nomial
+a1=(s^4)+(s^3)+(2*s^2)+(2*s)+1 //quotient poly nomial
+b1=coeff(a1)
+n1=length(b1)
+R1=routh_t(a1)
+disp(R1,'the routh array for quotient poly nomial is;')
+ap=s^2+1
+r=roots(ap)
+disp(r,'the roots are')
+disp('the system is unstable')
+disp('two roota on imaginary axis ,two roots on right half of s plane and two roots lie in left half of s plane')
-- 
cgit