dwz

function dwz(r,l,c,u)

for q=1:3
z(q)=r(q)+(l(q)-c(q))*i;
endfor
z

z23=(z(2)*z(3))/(z(2)+z(3))
zz=z(1)+z23
I(1)=u/zz
u23=I(1)*z23
I(2)=u23/z(2)
I(3)=u23/z(3)

for q=1:3
ir(q)=findz(I(q));
endfor
ir

Ur=ir.*r
Ul=ir.*l
Uc=ir.*c

Sfull=I(1)*u
for q=1:3
P(q)=ir(q)^2*r(q);
Q(q)=ir(q)^2*(l(q)-c(q));
endfor
SumP=sum(P)
SumQ=sum(Q)
endfunction

function ret=findz(c)
ret=sqrt(real(c)^2+imag(c)^2);
endfunction