summaryrefslogtreecommitdiff
path: root/104/CH8/EX8.14/8_14.sce
diff options
context:
space:
mode:
Diffstat (limited to '104/CH8/EX8.14/8_14.sce')
-rwxr-xr-x104/CH8/EX8.14/8_14.sce15
1 files changed, 15 insertions, 0 deletions
diff --git a/104/CH8/EX8.14/8_14.sce b/104/CH8/EX8.14/8_14.sce
new file mode 100755
index 000000000..93a3f4988
--- /dev/null
+++ b/104/CH8/EX8.14/8_14.sce
@@ -0,0 +1,15 @@
+//breakaway points
+s=%s
+sys=syslin('c',1/(s*(s^2+2*s+2)))
+evans(sys)
+syms s
+d=derivat(sys)
+n=numer(d)
+a=roots(n) //a=breakaway points
+disp(a,"breakaway points=")
+for i=1:2
+ K=-a(i,1)^2+2*a(i,1)+2
+ disp(a(i,1),"s=")
+ disp(K,"K=")
+end
+printf("if K is complex then point is not a break away point") \ No newline at end of file