updated the code

1 files changed, 24 insertions, 23 deletions

diff --git a/854/CH1/EX1.2/Example1_2.sce b/854/CH1/EX1.2/Example1_2.sce
index 225fc4169..ac185a2e1 100755
--- a/854/CH1/EX1.2/Example1_2.sce
+++ b/854/CH1/EX1.2/Example1_2.sce
@@ -1,23 +1,24 @@
-//clear//

-//Caption: Program to find the phase angle between two vectors

-//Example1.2

-//page 11

-clc;

-Q = [4,5,2]; //point Q

-x = Q(1);

-y = Q(2);

-z = Q(3);

-G = [y,-2.5*x,3]; //vector field

-disp(G,'G(rQ) =')

-aN = [2/3,1/3,-2/3]; //unit vector- direction of Q

-G_dot_aN = dot(G,aN); //dot product of G and aN

-disp(G_dot_aN,'G.aN =')

-G_dot_aN_aN = G_dot_aN*aN;

-disp(G_dot_aN_aN,'(G.aN)aN=')

-teta_Ga = Phase_Angle(G,aN) //phase angle between G and unit vector aN

-disp(teta_Ga,'phase angle between G and unit vector aN in degrees =')

-//Result

-// G(rQ) =       5.  - 10.    3.  

-// G.aN =      - 2.  

-// (G.aN)aN =     - 1.3333333  - 0.6666667    1.3333333  

-// phase angle between G and unit vector aN in degrees =   99.956489  

+//Caption: Program to find the phase angle between two vectors
+//Example1.2
+//page 11
+clc;
+clear
+Q = [4,5,2]; //point Q
+x = Q(1);
+y = Q(2);
+z = Q(3);
+G = [y,-2.5*x,3]; //vector field
+disp(G,'G(rQ) =')
+aN = [2/3,1/3,-2/3]; //unit vector- direction of Q
+G_dot_aN = sum(G.*aN); //dot product of G and aN
+disp(G_dot_aN,'G.aN =')
+G_dot_aN_aN = G_dot_aN*aN;
+disp(G_dot_aN_aN,'(G.aN)aN=')
+teta_Ga = 2 * atand(norm(G*norm(aN) - norm(G)*aN) / norm(G * norm(aN) + norm(G) * aN))
+ //phase angle between G and unit vector aN
+disp(teta_Ga,'phase angle between G and unit vector aN in degrees =')
+//Result
+// G(rQ) =       5.  - 10.    3.  
+// G.aN =      - 2.  
+// (G.aN)aN =     - 1.3333333  - 0.6666667    1.3333333  
+// phase angle between G and unit vector aN in degrees =   99.956489
\ No newline at end of file