clc

open=input('preiona enter')
//Variables de entrada//
fp=0.95
fprad=acos(fp)
S=10                             //Potencia aparente del sistema//
V=500                         //Voltaje de la línea//
f=50                              //Frecuencia del sistema//
l=400                          //Largo de la línea//
cir=1                             //Número de circuitos//
seccion=765                     //Sección del conductor en mm2//
n=4                               //Número de subconductores//
G=0                   //Conductancia de la línea//
d=0.4                            //Distancia entre conductores//
D12=14                      //Distancia entre conductores//
D23=14                      //Distancia entre conductores//
D13=28                           //Distancia entre conductores//
RL=(0.0402/n)                   //Resistencia de la línea por subconductor//
V2=complex(0.875,0)               //Voltaje en pu del consumo//
I2=complex(1.142857*cos(fprad),1.142857*sin(-fprad))          //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//
XT1=complex(0.00266666,0.007999999)
XT2=complex(0.03,0.12)

//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)
B=((sinh(Gamma*l)*Zc)/Zb)
C=((sinh(Gamma*l))/Zc)*Zb
D=A
T=[A,B;C,D]
F=1/sqrt(real(A)^2+(imag(A)^2))    //Efecto ferranti en la línea//

//Calculo de la matriz de entrada//
MT1=[1,XT1;0,1]
MT2=[1,XT2;0,1]
MC=[V2;I2]
MTG=MT1*T*MT2*MC


V1=MTG(1.1)
I1=MTG(2.1)

//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))-(sqrt(imag(V2)^2+real(V2)^2)))/(sqrt((imag(V1)^2)+(real(V1))^2)))*100
S1=V1*I1

disp(I1,'La corriente I1 es:')
disp(V1,'El voltaje V1 es:')
disp(Rend,'El rendimiento del sistema es:')
disp(Reg,'La regulación del sistema es:')
disp(S1,'La potencia aparente es:')
disp (fp1,'factor de potencia')
disp(r,'el radio')