close all
clear all
clc

I2= (1/(sqrt(3)))*(0.5+1.5i); %odb. prúd v uzle c.2
I3= (1/(sqrt(3)))*(0.8+0.2i);
I4= (1/(sqrt(3)))*(1.9+1.6i);

U1=400; %U1=Ub
Ubil=400;

Yz12= 1/(2+9i); %impedancie vetiev
Yz14= 1/(1+5i);
Yz23= 1/(1+8i);
Yz34= 1/(5+4i);

Yadm= [ Yz12+Yz14  -Yz12          0      -Yz14;
        -Yz12      Yz12+Yz23     -Yz23     0 ;
          0          -Yz23        Yz23+Yz34   -Yz34; 
         -Yz14      0              -Yz34      Yz14+Yz34;
       ]
   
Y21= Yadm(2:2,1:1);  
Y31= Yadm(3:3,1:1);  
Y41= Yadm(4,1);

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

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

Y42= Yadm(4,2);
Y43= Yadm(4,3);
Y44= Yadm(4,4);




g22= real(Y22);
g23= real(Y23);
g24= real(Y24);

b22= -imag(Y22);
b23= -imag(Y23);
b24= -imag(Y24);

g32= real(Y32);
g33= real(Y33);
g34= real(Y34);

b32= -imag(Y32);
b33= -imag(Y33);
b34= -imag(Y34);


g42= real(Y42);
g43= real(Y43);
g44= real(Y44);

b42= -imag(Y42);
b43= -imag(Y43);
b44= -imag(Y44);

g21= real(Y21);
g31= real(Y31);
g41= real(Y41);

b21= -imag(Y21);
b31= -imag(Y31);
b41= -imag(Y41);    

% verzia s Y22+Y23= -Y21; Y32+Y33= -Y31;
X= [ g22  g23 g24  b22  b23 b24  ((sqrt(3)*real(I2))-(Ubil*g21));
     g32  g33 g34  b32  b33 b34   ((sqrt(3)*real(I3))-(Ubil*g31));
     g42  g43 g44  b42  b43 b44   ((sqrt(3)*real(I4))-(Ubil*g41));
    -b22 -b23 -b24 g22  g23 g24   ((sqrt(3)*imag(I2))+(Ubil*b21));
    -b32 -b33 -b34 g32  g33 g34  ((sqrt(3)*imag(I3))+(Ubil*b31));
     -b42 -b43 -b44 g42  g43 g44  ((sqrt(3)*imag(I4))+(Ubil*b41));
    ]

U= rref(X)


%NE= [Y22  Y23 Y24 ((sqrt(3)*(I2))-(Y21*Ubil));
 %    Y32  Y33 Y34 ((sqrt(3)*(I3))-(Y31*Ubil));
  %   Y42   Y43 Y44 ((sqrt(3)*(I4))-(Y41*Ubil));
%     ]
%rref(NE)

E= [Y22  Y23 Y24 ((sqrt(3)*(I2))+(Y22*Ubil)+(Y23*Ubil)+(Y24*Ubil));
    Y32  Y33 Y34 ((sqrt(3)*(I3))+(Y32*Ubil)+(Y33*Ubil)+(Y34*Ubil));
     Y42  Y43 Y44 ((sqrt(3)*(I4))+(Y42*Ubil)+(Y43*Ubil)+(Y44*Ubil));
    ]
rref(E)