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)