diff options
Diffstat (limited to '257/CH8/EX8.25/example_8_25.sce')
-rw-r--r-- | 257/CH8/EX8.25/example_8_25.sce | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/257/CH8/EX8.25/example_8_25.sce b/257/CH8/EX8.25/example_8_25.sce new file mode 100644 index 000000000..0d03de8b8 --- /dev/null +++ b/257/CH8/EX8.25/example_8_25.sce @@ -0,0 +1,23 @@ +s=%s
+P=s^5-s^4-2*s^3+2*s^2-8*s+8
+routh=routh_t(P)
+disp(routh)
+r=coeff(P)
+n=length(r)
+c=0;
+for i=1:n
+if (routh(i,1)<0)
+c=c+1;
+ end
+ end
+ if(c>=1)
+ printf("there are %d roots on RHS",c)
+ else printf("there are no roots on RHS")
+ end
+
+ disp("s^2 is")
+ disp(roots(routh(2,:)))
+
+
+
+
\ No newline at end of file |