% Puntos P=[x1 x2 x3; y1 y2 y3]
P=[0 1 0;0 0 1];


% delta = [P1-P2 P3-P2] = [delta12 delta13]
delta=[P(:,2)-P(:,1) P(:,3)-P(:,1)]
delta=delta(:,1)+i*delta(:,2)

% teta = [fi12 chi12 teta12; fi13 chi13 delta13];
teta=(pi/180)*[30 60 30;60 90 45]

b=delta
A=[exp(j*teta(1,1))-1 exp(j*teta(1,3))-1;exp(j*teta(2,1))-1 exp(j*teta(2,3))-1];

Z23=A\b
z2=Z23(1,:);
z3=Z23(2,:);
z1=(P(1,1)+j*P(2,1))-z2-z3;


A=[exp(j*teta(1,3))-1 exp(j*teta(1,2))-1;exp(j*teta(2,3))-1 exp(j*teta(2,2))-1];
Z54=A\b
z4=Z54(1,:);
z5=Z54(2,:);

z6=(P(1,1)+j*P(2,1))-z5-z4-z1;

[z1;z2;z3;z4;z5;z6]

%%%%%
fill(real([z1+z2 z1+z2+z3 z1+z6+z5+z4 z1+z6+z5]),imag([z1+z2 z1+z2+z3 z1+z6+z5+z4 z1+z6+z5]),'c')
hold on
plot(z1,'ro');
plot(z1,'r*');
plot(z1+z6,'ro');
plot(z1+z6,'r*');


plot(z1+z2,'o');
plot(z1+z6+z5,'o');

plot([z1 z1+z2 z1+z6+z5 z1+z6])
plot([z1+z2 z1+z2+z3 z1+z6+z5+z4 z1+z6+z5])
plot(P(1,:),P(2,:),'.')
axis equal