E=(6/sqrt(2))*exp(j*pi/4)
Iz1=(6/sqrt(2))*exp(-j*pi/4)*1e-3
r1=2500;
r2=1000;
r3=500;
r4=2000;
c=0.0000002;
l=0.2;
f=800;
w=2*pi*f;
Iz2=E/r3;

Y(1,1)=((1/r2)+(j*w*c)+(1/r1));
Y(1,2)=-j*w*c;
Y(1,3)=-1/r1;
Y(2,1)= -j*w*c;
Y(2,2)=((j*w*c)+(1/r4)+(1/(j*w*l)));
Y(2,3)=-1/r4;
Y(3,1)=-1/r1;
Y(3,2)=-1/r4;
Y(3,3)=((1/r3)+(1/r4)+(1/r1));

Im=[-Iz1;0;Iz2+Iz1];
Vm=Y\Im
I=(Vm(2)-Vm(3))*(1/r4)*1e3
Im=I*sqrt(2)
M=abs(Im)
K=angle(Im)*180/pi
abs(E)
abs(Iz1)