diff options
author | priyanka | 2015-06-24 15:03:17 +0530 |
---|---|---|
committer | priyanka | 2015-06-24 15:03:17 +0530 |
commit | b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b (patch) | |
tree | ab291cffc65280e58ac82470ba63fbcca7805165 /104/CH8/EX8.11/8_11.sce | |
download | Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.gz Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.tar.bz2 Scilab-TBC-Uploads-b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b.zip |
initial commit / add all books
Diffstat (limited to '104/CH8/EX8.11/8_11.sce')
-rwxr-xr-x | 104/CH8/EX8.11/8_11.sce | 16 |
1 files changed, 16 insertions, 0 deletions
diff --git a/104/CH8/EX8.11/8_11.sce b/104/CH8/EX8.11/8_11.sce new file mode 100755 index 000000000..d10efd0d1 --- /dev/null +++ b/104/CH8/EX8.11/8_11.sce @@ -0,0 +1,16 @@ +//breakaway points
+s=%s
+sys=syslin('c',(s+4)/(s*(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)*(a(i,1)+2)/(a(i,1)+4)
+ disp(a(i,1),"s=")
+ disp(K,"K=")
+end
+printf("if K is positive breakaway point lies on RL or else on CRL")
+
|