diff options
Diffstat (limited to '27/CH2/EX2.7.2/Example_2_7_2.sce')
-rwxr-xr-x | 27/CH2/EX2.7.2/Example_2_7_2.sce | 23 |
1 files changed, 23 insertions, 0 deletions
diff --git a/27/CH2/EX2.7.2/Example_2_7_2.sce b/27/CH2/EX2.7.2/Example_2_7_2.sce new file mode 100755 index 000000000..16e1b7ee2 --- /dev/null +++ b/27/CH2/EX2.7.2/Example_2_7_2.sce @@ -0,0 +1,23 @@ +//Example 2.7.2 Page 31
+//Non-Linear Dynamics and Chaos, First Indian Edition Print 2007
+//Steven H. Strogatz
+
+clear;
+clc;
+close;
+set(gca(),"auto_clear","off") //hold on
+
+//Given : x(dot) = f(x) = x - x^3.
+//On integrating we get --> V(x) = -(1/2)x^2 + (1/4)x^4 + C ; C=0.
+
+//Now plotting V(x) v/s x ; and observe fix points and their Stabilities
+
+for x = -2:0.05:2
+ V = -(1/2)*(x^2) + (1/4)*(x^4);
+ plot2d(x,V,style=-2)
+end
+plot2d(0,0,style=-3) //Just to show that the fixed point is UnStable.
+plot2d(-1,-1/4,style=-4) //Just to show that the fixed point is Stable.
+plot2d(1,-1/4,style=-4) //Just to show that the fixed point is Stable.
+set(gca(),"grid",[2,5])
+xtitle("Double-Well Potential Diagram (Bistable)","x-Axis(x)","y-Axis (V(x))")
\ No newline at end of file |