close all
clear all
clc

I2= (1/(sqrt(3)))*(-0.2+0.3i);
I3= (1/(sqrt(3)))*(-0.4+0.3i);

U1=115;
Ubil=115;

Yz12= 1/(4+3i);
Yz13= 1/(3+4i);
Yz23= 1/(5);

Yadm= [Yz12+Yz13 -Yz12 -Yz13;
       -Yz12 Yz12+Yz23 -Yz23;
       -Yz13 -Yz23 Yz23+Yz13;
       ];
   
Y11= Yadm(1,1);
Y12= Yadm(1,2);
Y13= Yadm(1,3);
Y21= Yadm(2,1);
Y22= Yadm(2,2);
Y23= Yadm(2,3);
Y31= Yadm(3,1);
Y32= Yadm(3,2);
Y33= Yadm(3,3);
   
% 0 iteracia
U3i= 112 + i*0; %treba tipnit hodnotu
U2i= 113 + i*0;

% 1 iteracia
U2= inv(Y22)*(sqrt(3)*I2 - Y21*U1 - Y23*U3i)
U3= inv(Y33)*(sqrt(3)*I3 - Y31*U1 - Y32*U2i)
real(U3-U3i)
real(U2-U2i)
% 2 iteracia
U3i= U3;
U2i= U2;
   
U2= inv(Y22)*(sqrt(3)*I2 - Y21*U1 - Y23*U3i)
U3= inv(Y33)*(sqrt(3)*I3 - Y31*U1 - Y32*U2i)
real(U3-U3i)
real(U2-U2i)
% 3 iteracia
U3i= U3;
U2i= U2;
   
U2= inv(Y22)*(sqrt(3)*I2 - Y21*U1 - Y23*U3i)
U3= inv(Y33)*(sqrt(3)*I3 - Y31*U1 - Y32*U2i)
real(U3-U3i)
real(U2-U2i)

% 4 iteracia (pridavam dovtedy kym nesplnim presnost a odmocnina z 3 bud je alebo nie)
%+ tabulka k, Ua2, Ur2