%% Nombre del Estudiante %%
clear all
clc
%% Zprim= matriz primitiva de impedancias 
Zprim=[0.4013+1.4133*i 0.0953+0.8515*i 0.0953+0.7266*i 0.0953+0.7524*i;...
    0.0953+0.8515*i 0.4013+1.4133*i 0.0953+0.7802*i 0.0953+0.7865*i;...
    0.0953+0.7266*i 0.0953+0.7802*i 0.4013+1.4133*i 0.0953+0.7674*i;...
    0.0953+0.7524*i 0.0953+0.7865*i 0.0953+0.7674*i 0.6873+1.5465*i];
%% Definición de sub-matrices de la matriz primitiva de impedancias
Zij= Zprim(1:3,1:3);
Zin= Zprim(1:3,4);
Znj= Zprim(4,1:3);
Znn= Zprim(4,4);
%% Determinar la Matriz de impedancias de fase Zabc en Ohm/milla

%% Matriz primitiva de coeficientes de potencial Prim en milla/microF
Pprim=[84.56 35.1522 23.7147 25.2469;...
    35.1522 84.56 28.6058 28.359;...
    23.7147 28.6058 84.56 26.6131;...
    25.2469 28.359 26.6131 85.6659];
%% Definición de sub-matrices de la matriz primitiva de coeficientes de
% potencial
Pij= Pprim(1:3,1:3);
Pin= Pprim(1:3,4);
Pnj= Pprim(4,1:3);
Pnn= Pprim(4,4);
%% Determinar la Matriz de coeficientes de potencial de fase
Pabc=
%% Determinar la Matriz de Capacitancias Cabc
Cabc= 
%% Matriz de admitancias Yabc microS/milla
Yabc= 
%% Determinar las matrices generalizadas a,b,c y d
a=
b=
c=
d=
%% Datos de tensiones y corrientes en la carga (Nodo m)
t=(cosd(120)+i*sind(120));
Vabcm=7199.56*[1;conj(t);t];
DI=cosd(-25.84)+i*sind(-25.84);
Iabcm=277.79*[DI;conj(t)*DI;t*DI];
%% Determianar las tensiones y corrientes en la fuente Nodo n
Vabcn=
Iabcn=