1 files changed, 16 insertions, 0 deletions
diff --git a/104/CH8/EX8.13/8_13.sce b/104/CH8/EX8.13/8_13.sce
new file mode 100755
index 000000000..d09a64a9b
--- /dev/null
+++ b/104/CH8/EX8.13/8_13.sce
@@ -0,0 +1,16 @@
+//breakaway points
+s=%s
+sys=syslin('c',1/(s*(s+4)*(s^2+4*s+20)))
+evans(sys)
+syms s
+d=derivat(sys)
+n=numer(d)
+a=roots(n)    //a=breakaway points
+disp(a,"breakaway points=")
+for i=1:3
+    K=-a(i,1)*(a(i,1)+4)*(a(i,1)^2+4*a(i,1)+20)
+        disp(a(i,1),"s=")
+        disp(K,"K=")
+end
+printf("if K is positive breakaway point lies on RL or else on CRL")
+    
+\ No newline at end of file