Modified the code

1 files changed, 0 insertions, 40 deletions

diff --git a/1340/CH6/EX6.4/6_4.sce b/1340/CH6/EX6.4/6_4.sce
deleted file mode 100755
index 0c494884e..000000000
--- a/1340/CH6/EX6.4/6_4.sce
+++ /dev/null
@@ -1,40 +0,0 @@
-clc;// in-built fuction routh_t can be used to generate the routh table

-s = poly(0,'s');

-tf = syslin('c',10/(s^5+7*s^4+6*s^3+42*s^2+8*s+56));

-deno = denom(tf);

-coef = coeff(deno);

-routh = [coef([6,4,2]);coef([5,3,1])];

-// we will get a row of all zeroes

-T = routh(2,:)/7;

-coef1 = coeff(T);

-// auxillary polynomial s^2+6*s+8 generation

-second = poly([coef1(3) 0 coef1(2) 0 coeff(1)],"s","coeff");disp(second);

-

-aux = derivat(second);//auxillary polynomial

-len = coeff(aux);

-routh = [routh;len(4) len(2) 0];

-disp(routh);

-t = routh(2:3,1:3);

-det1 = det(t(1:2,1:2))/t(2,1);

-det2 = -(t(1,1)*t(2,3)-t(2,1)*t(1,3))/t(2,1);

-routh = [routh;-det1 det2 0];

-t1 = routh(3:4,1:2);

-det3 = det(t1(1:2,1:2))/t1(2,1);

-routh = [routh;-det3 0 0];

-t2 = routh(4:5,1:2);

-det4 = det(t2(1:2,1:2))/t2(2,1);

-routh = [routh;-det4 0 0];

-disp(routh)

-c = 0;

-for k = 1:length(coef)

-    if(routh(k,1)<0)

-        c =c +1;

-        end

-end

-

-if(c>=1)

-    printf("system is unstable")

-else printf("system is stable,hence no poles in RHP")

-    end

-

-