close all
clear all
clc

Ub= 115;
U1= Ub;
epsilon= 0.01;
 
Z12= (3 + i*4);
S2z= (-44 - i*33);

Y= 1/Z12;

Yadm= [ Y  -Y;
       -Y  Y; ]
   
Y11= Yadm(1:1,1:1);
Y12= Yadm(1:1,2:2);
Y21= Yadm(2:2,1:1);
Y22= Yadm(2:2,2:2);

g22= real(Y22);
b22= -imag(Y22);

g21= real(Y21);
b21= -imag(Y21);

% 0 iteracia
U2= 110-0.0i;      %+i*0;
Ua2= real(U2);
Ur2= imag(U2);

P= g22*(Ua2^2+Ur2^2)+g21*Ua2*Ub-b21*Ur2*Ub
Q= b22*(Ua2^2+Ur2^2)+g21*Ur2*Ub+b21*Ua2*Ub

deltaP2= real(S2z)-P
deltaQ2= imag(S2z)-Q

a= 2*g22*Ua2 + g21*Ub
b= 2*g22*Ur2 - b21*Ub

c= 2*b22*Ua2 + b21*Ub
d= 2*b22*Ur2 + g21*Ub

U= [a b deltaP2;
    c d deltaQ2];

O= rref(U)   %prvý člen hovorí o Ua (Aktivnej) a druhý o Ur (jalovej zložke napatia U)

% 1 iteracia
Ua= real(U2)+O(1,3)
Ur= imag(U2)+O(2,3)

P= g22*(Ua^2+Ur^2)+g21*Ua*Ub-b21*Ur*Ub;
Q= b22*(Ua^2+Ur^2)+g21*Ur*Ub+b21*Ua*Ub;

deltaP2= real(S2z)-P   
deltaQ2= imag(S2z)-Q

a= 2*g22*Ua + g21*Ub;
b= 2*g22*Ur - b21*Ub;

c= 2*b22*Ua + b21*Ub;
d= 2*b22*Ur + g21*Ub;

U= [a b deltaP2;
    c d deltaQ2];

O= rref(U)

% 2 iteracia
Uai= Ua+O(1,3)
Uri= Ur+O(2,3)

P= g22*(Uai^2+Uri^2)+g21*Uai*Ub-b21*Uri*Ub;
Q= b22*(Uai^2+Uri^2)+g21*Uri*Ub+b21*Uai*Ub;

deltaP2= real(S2z)-P
deltaQ2= imag(S2z)-Q

a= 2*g22*Uai + g21*Ub;
b= 2*g22*Uri - b21*Ub;

c= 2*b22*Uai + b21*Ub;
d= 2*b22*Uri + g21*Ub;

U= [a b deltaP2;
    c d deltaQ2];

O= rref(U)

Uanove= Uai+O(1,3)
Urnove= Uri+O(2,3)