1 files changed, 28 insertions, 0 deletions
diff --git a/1938/CH6/EX6.3/6_3.sce b/1938/CH6/EX6.3/6_3.sce
new file mode 100755
index 000000000..ab19695f2
--- /dev/null
+++ b/1938/CH6/EX6.3/6_3.sce
@@ -0,0 +1,28 @@
+clc,clear
+printf('Example 6.3\n\n')
+
+//note that a new function p2z has been defined below for direct representation of complex numbers in polar form
+function [FUN] = p2z(RRRR,Theeeta)
+ FUN = RRRR.*exp(%i*%pi*Theeeta/180.);
+endfunction
+
+Z1=complex(0,3) //impedance of alternator 1
+Z2=complex(0,4) //impedance of alternator 2
+Z=6 //load
+E1=p2z(220,0) //induced emf vector on no load
+E2=p2z(220,10)//induced emf vector on no load
+
+I1=((E1-E2)*Z+E1*Z2)/(Z*(Z1+Z2)+Z1*Z2)
+I2=((E2-E1)*Z+E2*Z1)/(Z*(Z1+Z2)+Z1*Z2)
+
+phi1=phasemag(I1)  //Phasemag returns the angle of complex number in degrees
+phi2=phasemag(I2)  //Phasemag returns the angle of complex number in degrees
+
+I=I1+I2
+V=I*Z //Terminal voltage
+printf('(i) Terminal voltage is %.1f volts at %.2f degrees\n',abs(V),phasemag(V))
+printf('(ii) Currents are %.2f A at %.2f degrees and %.2f A at %.2f degrees\n     Total current is %.2f A at %.2f degrees  ',abs(I1),phasemag(I1),abs(I2),phasemag(I2),abs(I),phasemag(I))
+
+P1=abs(V)*abs(I1)*cosd(phi1)
+P2=abs(V)*abs(I2)*cosd(phi2)
+printf('\n(iii)Power delivered is %.2f watts and %.2f watts',P1,P2)