%% Nombre del Estudiante %%
clear all
clc
%% Para el segmento 1-2
Z12=[0.1414+0.5353*i 0.0361+0.3225*i 0.0361+0.2752*i;...
    0.0361+0.3225*i 0.1414+0.5353*i 0.0361+0.2955*i;...
    0.0361+0.2752*i 0.0361+0.2955*i 0.1414+0.5353*i];
a1=eye(3);
b1=Z12;
d1=a1;
c1=zeros(3,3);
%% Para el segmento 3-4
Z34=[0.1907+0.5035*i 0.0607+0.2302*i 0.0598+0.1751*i;...
    0.0607+0.2302*i 0.1939+0.4885*i 0.0614+0.1931*i;...
    0.0598+0.1751*i 0.0614+0.1931*i 0.1921+0.4970*i];
a2=eye(3);
b2=Z34;
d2=a2;
c2=zeros(3,3);
%% Para el transformador
Ztabc=[0.0288+0.1728*i 0 0;0 0.0288+0.1728*i 0; 0 0 0.0288+0.1728*i];
Zt=0.0288+0.1728*i;
at=[0 -3.4639 -1.7319;...
    -1.7319 0 -3.4639;...
    -3.4639 -1.7319 0];
bt1=[0 2*Zt Zt;Zt 0 2*Zt;2*Zt Zt 0];
bt=(-5.1958/3).*bt1;
ct=zeros(3,3);
dt=(1/5.1958).*[1 -1 0;0 1 -1;-1 0 1];
At=(1/5.1958).*[1 0 -1;-1 1 0;0 -1 1];
Bt=Ztabc;
%% Tensiones de fase en el bus infinito
t=(cosd(120)+i*sind(120));
En=7199.6*[1;conj(t);t];
%% Potencia de la carga
Sabc=1000*[750*exp(i*0.55481);1000*exp(i*0.45103);1250*exp(i*0.31756)];
%% Primera iteración
t=(cosd(120)+i*sind(120));
Vln4=2400*[1;conj(t);t];