1 files changed, 34 insertions, 0 deletions
diff --git a/27/CH6/EX6.6.3/Example_6_6_3.sce b/27/CH6/EX6.6.3/Example_6_6_3.sce
new file mode 100755
index 000000000..86ba7b6c6
--- /dev/null
+++ b/27/CH6/EX6.6.3/Example_6_6_3.sce
@@ -0,0 +1,34 @@
+clear;
+clear;
+clc;
+close;
+//obtain the General Jacobian "A";
+
+x1= %pi/2; y1= %pi/2;
+x2= %pi/2; y2= -%pi/2;
+x3= -%pi/2; y3= %pi/2;
+x4= -%pi/2; y4= -%pi/2;
+
+A1 = [2*sin(x1) sin(y1); sin(x1) 2*sin(y1)]
+t1=det(A1)
+tau1=trace(A1)
+d1 = ((tau1)^2) - 4*t1
+A2 = [2*sin(x2) sin(y2); sin(x2) 2*sin(y2)]
+t2=det(A2)
+tau2=trace(A2)
+d2 = ((tau2)^2) - 4*t2
+A3 = [2*sin(x3) sin(y3); sin(x3) 2*sin(y3)]
+t3=det(A3)
+tau3=trace(A3)
+d3 = ((tau3)^2) - 4*t3
+A4 = [2*sin(x4) sin(y4); sin(x4) 2*sin(y4)]
+t4=det(A4)
+tau4=trace(A4)
+d4 = ((tau4)^2) - 4*t4
+
+
+disp("From the above information we come to following conclusion:")
+disp("1.    (pi/2,pi/2) --> Unstable Node.")
+disp("2.    (pi/2,-pi/2) --> Saddle.")
+disp("3.    (-pi/2,pi/2) --> Saddle.")
+disp("4.    (-pi/2,-pi/2) --> Stable Node.")
+\ No newline at end of file