close all
clear all
clc

%Zadanie
U1= 115;
delta1= 0;

U2= 115.415;
delta2= 0.1352;

U3= 109.712;
delta3= -2.1553;

Y12= (-0.02 + i*0.04);
Y13= (-0.0133 + i*0.0267);
Y23= (-0.0138 + i*0.0345);
g12= real(Y12);
g13= real(Y13);
g23= real(Y23);
b12= -imag(Y12);
b13= -imag(Y13);
b23= -imag(Y23);

g11= -(g12+g13);
g22= -(g12+g23);
g33= -(g13+g23);
b11= -(b12+b13);
b22= -(b12+b23);
b33= -(b13+b23);

g21= g12;
g31= g13;
g32= g23;

G= [g11 g12 g13
    g21 g22 g23
    g31 g32 g33];

b21= b12;
b31= b13;
b32= b23;

B= [b11 b12 b13
    b21 b22 b23
    b31 b32 b33];

U1ar= U1*exp(i*delta1*pi/180);
U1a= real(U1ar);
U1r= imag(U1ar);

U2ar= U2*exp(i*delta2*pi/180);
U2a= real(U2ar);
U2r= imag(U2ar);

U3ar= U3*exp(i*delta3*pi/180);
U3a= real(U3ar);
U3r= imag(U3ar);

Ua= [U1a
     U2a
     U3a]
 
Ur= [U1r
     U2r
     U3r]

I12a = (1/sqrt(3))*((U1a-U2a)*(-g12)+(U1r-U2r)*(-b12))
I12r = (1/sqrt(3))*((-(U1a-U2a)*(-b12))+(U1r-U2r)*(-g12))
I21a = -I12a
I21r = -I12r

I13a = (1/sqrt(3))*((U1a-U3a)*(-g13)+(U1r-U3r)*(-b13))
I13r = (1/sqrt(3))*((-(U1a-U3a)*(-b13))+(U1r-U3r)*(-g13))
I31a = -I13a
I31r = -I13r

I23a = (1/sqrt(3))*((U2a-U3a)*(-g23)+(U2r-U3r)*(-b23))
I23r = (1/sqrt(3))*((-(U2a-U3a)*(-b23))+(U2r-U3r)*(-g23))
I32a = -I23a
I32r = -I23r

P12= sqrt(3)*(U1a*I12a+U1r*I12r)
Q12= sqrt(3)*(U1r*I12a-U1a*I12r)
P21= sqrt(3)*(U2a*I21a+U2r*I21r)
Q21= sqrt(3)*(U2r*I21a-U2a*I21r)

P13= sqrt(3)*(U1a*I13a+U1r*I13r)
Q13= sqrt(3)*(U1r*I13a-U1a*I13r)
P31= sqrt(3)*(U3a*I31a+U3r*I31r)
Q31= sqrt(3)*(U3r*I31a-U3a*I31r)

P23= sqrt(3)*(U2a*I23a+U2r*I23r)
Q23= sqrt(3)*(U2r*I23a-U2a*I23r)
P32= sqrt(3)*(U3a*I32a+U3r*I32r)
Q32= sqrt(3)*(U3r*I32a-U3a*I32r)

deltaP12 = P12+P21
deltaP13 = P13+P31
deltaP23 = P23+P32

deltaQ12 = Q12+Q21
deltaQ13 = Q13+Q31
deltaQ23 = Q23+Q32

deltaPc = deltaP12+deltaP13+deltaP23
deltaQc = deltaQ12+deltaQ13+deltaQ23