1 files changed, 44 insertions, 0 deletions
diff --git a/1445/CH2/EX2.47/ch2_ex_47.sce b/1445/CH2/EX2.47/ch2_ex_47.sce
new file mode 100644
index 000000000..807493965
--- /dev/null
+++ b/1445/CH2/EX2.47/ch2_ex_47.sce
@@ -0,0 +1,44 @@
+//CHAPTER 2- STEADY-STATE ANALYSIS OF SINGLE-PHASE A.C. CIRCUIT
+//Example 47
+
+disp("CHAPTER 2");
+disp("EXAMPLE 47");
+
+//VARIABLE INITIALIZATION
+r=10;                          //in Ohms
+xl=8.66;                       //in Ohms
+I=5-(%i*10);                   //in Amperes
+
+//SOLUTION
+z=r+(%i*(xl));
+//function to convert from rectangular form to polar form
+function [mag,angle]=rect2pol(x,y);
+mag=sqrt((x^2)+(y^2));         
+angle=atan(y/x)*(180/%pi);     //to convert the angle from radians to degrees
+endfunction;
+[z,angle_z]=rect2pol(real(z),imag(z)); 
+[I,angle_I]=rect2pol(real(I),imag(I)); 
+
+//solution(i)
+v=I*z;
+angle_v=angle_I+angle_z;
+disp(sprintf("(i) The applied voltage is %f V, %f degrees",v,angle_v));
+
+//solution (ii)
+phi=angle_I-angle_v;
+pf=cos(phi*(%pi/180));
+disp(sprintf("(ii) The power factor is %f (lagging)",pf));
+
+//solution(iii)
+s=v*I;
+angle_s=angle_v-angle_I;
+//function to convert from polar form to rectangular form
+function [x,y]=pol2rect(mag,angle);
+x=mag*cos(angle*(%pi/180));  //to convert the angle from degrees to radians
+y=mag*sin(angle*(%pi/180));
+endfunction;
+[p,q]=pol2rect(s,angle_s);
+disp(sprintf("(iii) The active power is %f W",p));
+disp(sprintf("      The reactive power is %f VAR",q));
+
+//END