clc
//Variables de entrada//

S=100                             //Potencia aparente del sistema//
V=110                             //Voltaje de la línea//
f=50                              //Frecuencia del sistema//
l=400                             //Largo de la línea//
cir=2                             //Número de circuitos//
seccion=765                       //Sección del conductor en mm2//
n=4                               //Número de subconductores//
G=31.5*10^-9                        //Conductancia de la línea//
d=0.45                            //Distancia entre conductores//
D12=4.70744                       //Distancia entre conductores//
D23=4.70744                       //Distancia entre conductores//
D13=9.2                           //Distancia entre conductores//
RL=0.0402/n                       //Resistencia de la línea por subconductor//
V2=1                              //Voltaje en pu del consumo//
fp=0.95
         //Corriente en pu del consumo//
r=(sqrt(seccion/%pi))*(0.001)     //Radio del conductor//
W=(2*%pi*f)                       //Omega//
E0=8.8542*10^(-12)                //Epsilon en vacío//

//Condición dependiente del n //
if n==1 then reql=(%e^-0.25)*r
   elseif n==2 then reql=sqrt((%e^-0.25)*r*d)
   elseif n==3 then reql=((%e^-0.25)*r*(d^2))^(1/3)
   elseif n==4 then reql=1.09*((((%e)^-0.25)*r*(d^3))^(1/4))
   else reql=r
end

if n==1 then reqc=r
elseif n==2 then reqc=sqrt(r*d)
elseif n==3 then reqc= (r*(d^2))^(1/3)
elseif n==4 then reqc=1.09*((r*(d^3))^(1/4))
else reqc=r
end

Deq=(D12*D13*D23)^(1/3)            //Distancia media entre conductores//
Zb=(V*V)/S                         //Cálculo de la impedancia base de la línea//

//Ecuaciones de inductancia//

L=(2*10^(-7))*(log(Deq/reql))*1000 //Inductancia en H/KM//
XL=W*L                             //Reactancia en ohm/KM//
Z=complex(RL,XL)                   //Impedancia de la línea//

//Ecuaciones de conductancia//

c=((2*%pi*E0)/(log(Deq/reqc)))*1000//Capacitancia en F/KM//
B=W*c                              //Admitancia en S/KM//
Y=complex(G,B)                    //Conductancia de la línea//

//Ecuaciones para la matriz//

Gamma=sqrt(Z*Y)                   //Constante de propagación//
Zc=sqrt(Z/Y)                      //ni idea qué es//

A=cosh(Gamma*l)
B1=((sinh(Gamma*l)*Zc)/Zb)/cir
C=((sinh(Gamma*l))/Zc)*Zb*cir
D=A
T=[A,B1;C,D]
F=1/(sqrt(real(A)^2+(imag(A)^2)))    //Efecto ferranti en la línea//
fprad=acos(fp)

//Cálculo de los parámetros del generador//

I2=complex(1*cos(fprad),-1*sin(fprad))
V2=1
V1=(A*V2)+(B1*I2)
I1=(C*V2)+(D*I2)
S=V1*I1

//Cálculo de la potencia del generador//

phi1=(atan((imag(V1))/(real(V1))))  //Ángulo phi del voltaje 1//
phi2=(atan((imag(I1))/(real(I1))))  //Ángulo phi de la corriente 1//

phi12=phi1-phi2                     //Diferencia de ángulos entre V1 e I1//
fp1=cos(phi12)                      //Factor de potencia del generador//

I1r=real(I1)                        //Valor real de la corriente 1//
I1i=-(imag(I1))                        //Valor imaginario de la corriente 1//

P1=real(V1*complex(I1r,I1i))

//Cálculos de la potencia de salida//

phi3=atan(imag(V2)/real(V2))       //Ángulo phi del voltaje 2//
phi4=atan(imag(I2)/real(I2))       //Ángulo phi de la corriente 2//
phi34=phi3-phi4                    //Diferencia de ángulos entre V2 e I2//
fp2=cos(phi34)                     //Factor de potencia del consumo//
I2r=real(I2)                       //Valor real de la corriente 2//
I2i=imag(I2)                       //Valor imaginario de la corriente //
P2=real(V2*complex(I2r,-I2i))      //Potencia real del consumo//

Rend=P2/P1                         //Rendimiento del sistema//

Reg=(((sqrt(imag(V1)^2+real(V1)^2))-1)/(sqrt((imag(V1)^2)+(real(V1))^2)))*100

disp(S,'La potencia del generador es:')
disp(Rend,'El rendimiento del sistema es:')
disp(Reg,'La regulación del sistema es:')
disp(F,'El efecto ferranti es:')