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); %impedancia na admitanciu preratana
Yz13= 1/(3+4i);
Yz23= 1/(5);

Yadma= [(real(Yz12)+real(Yz13))  real(-Yz12)                   real(-Yz13);
        real(-Yz12)              (real(Yz12)+real(Yz23))       real(-Yz23);
        real(-Yz13)              real(-Yz23)                   (real(Yz13)+real(Yz23)); 
       ] %Re časť

Yadmr= [ (imag(Yz12)+imag(Yz13))  imag(-Yz12)                   imag(-Yz13);
        imag(-Yz12)              (imag(Yz12)+imag(Yz23))       imag(-Yz23);
        imag(-Yz13)              imag(-Yz23)                   (imag(Yz13)+imag(Yz23)); 
       ] %Im časť

Yadm= [[Yadma]+i*[Yadmr]]

Y21= Yadm(2:2,1:1);  
Y31= Yadm(3:3,1:1);  

Y22= Yadm(2:2,2:2);
Y23= Yadm(2:2,3:3);

Y32= Yadm(3:3,2:2);
Y33= Yadm(3:3,3:3);

g22= real(Y22);
g23= real(Y23);
b22= -imag(Y22);
b23= -imag(Y23);

g32= real(Y32);
g33= real(Y33);
b32= -imag(Y32);
b33= -imag(Y33);

g21= real(Y21);
g31= real(Y31);
b21= -imag(Y21);
b31= -imag(Y31);

% verzia s Y22+Y23= -Y21; Y32+Y33= -Y31;
X= [ g22  g23  b22  b23   ((sqrt(3)*real(I2))-(Ubil*g21));
     g32  g33  b32  b33   ((sqrt(3)*real(I3))-(Ubil*g31));
    -b22 -b23  g22  g23   ((sqrt(3)*imag(I2))+(Ubil*b21));
    -b32 -b33  g32  g33   ((sqrt(3)*imag(I3))+(Ubil*b31));
    ]

U= rref(X)

% verzia bez Y22+Y23= -Y21; Y32+Y33= -Y31;
T= [ g22  g23  b22  b23   ((sqrt(3)*real(I2))+(Ubil*g22)+(Ubil*g23));
     g32  g33  b32  b33   ((sqrt(3)*real(I3))+(Ubil*g32)+(Ubil*g33));
    -b22 -b23  g22  g23   ((sqrt(3)*imag(I2))-(Ubil*b22))-(Ubil*b23);
    -b32 -b33  g32  g33   ((sqrt(3)*imag(I3))-(Ubil*b32))-(Ubil*b33);
    ]

G= rref(T)