% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG54ic2\n');
fprintf('Esquema Linea Doble SR, Bifurcando a X,Y\n');
fprintf('5 Tramos\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('Eleccion: kE,L1,L2,L3,L4,L5,Sccs,Sccx,Sccr,Sccy\n');
fprintf('          Kx0,Kr0,Ky0,k0=zs0/zr0\n');
fprintf('Compensacion 21: Si, No\n');
fprintf('Falta Doble. Para cada Tramo: m=p=0...1\n');
fprintf('Para cada Tramo: Barra m=1, Barra p=2\n');
fprintf('Analisis: Extremos Sm, Sp\n');
%fprintf('Salida Grafica 1: Extremos m: Todas las Unidades (pu)\n');
%fprintf('                : Extremos p: Todas las Unidades (pu)\n');
%fprintf('Salida Grafica 2: Extremos m: Todas las Unidades (pua)\n');
%fprintf('                : Extremos p: Todas las Unidades (pua)\n');
%fprintf('Salida Grafica 3: Analisis Extremo m Izq, 67N\n');
%fprintf('                : Analisis Extremo m Der, 67N\n');
%fprintf('                : Analisis Extremo p Izq, 67N\n');
%fprintf('                : Analisis Extremo p Der, 67N\n');
%fprintf('Salida Grafica 4: Analisis Extremo m Izq, R-X\n');
%fprintf('                : Analisis Extremo m Der, R-X\n');
%fprintf('                : Analisis Extremo p Izq, R-X\n');
%fprintf('                : Analisis Extremo p Der, R-X\n');
fprintf('\n');
% ----------------------------------------------
%PG54ic2
datos=input('datos: standard=1, por input=2: ');
switch(datos)
case{1},
Un=220;
kE=1;
L1=100;
L2=100;
L3=100;
L4=100;
L5=100;
Sccs=5000;
Sccx=5000;
Sccr=5000;
Sccy=5000;
Kx0=1;
Kr0=1;
Ky0=1;
compens21=2;
% ----------------------------------------------
case{2},
Un=input('U(Fase-Fase) kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L1=input('Longitud Tramo 1 km:  ');
L2=input('Longitud Tramo 2 km:  ');
L3=input('Longitud Tramo 3 km:  ');
L4=input('Longitud Tramo 4 km:  ');
L5=input('Longitud Tramo 5 km:  ');
Sccs=input('Sccs(MVA): ');
Sccx=input('Sccx(MVA): ');
Sccr=input('Sccr(MVA): ');
Sccy=input('Sccy(MVA): ');
Kx0=input('Kx0: ');
Kr0=input('Kr0: ');
Ky0=input('Ky0: ');
compens21=input('Compensacion 21: Si=1, No=2 :');
end
% ----------------------------------------------
k0=input('k0=zs0/zr0: ');
tramo=input('Tramo en Falta, m=p: T1=1,T2=2,T3=3,T4=4,T5=5 :');

faltam=input('Falta en m: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :');
Rfm=input('Resistencia de Falta en m. Rf(Ohm) = ');
if faltam==7|faltam==8|faltam==9|faltam==11
Rgm=input('Resistencia de Falta en m. Rg(Ohm) = ');
else 
   Rgm=0;
end

faltap=input('Falta en p: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :');
Rfp=input('Resistencia de Falta en p. Rf(Ohm) = ');
if faltap==7|faltap==8|faltap==9|faltap==11
Rgp=input('Resistencia de Falta en p. Rg(Ohm) = ');
else 
   Rgp=0;
end
%--------------------------------------------------------------
pfaltam=1;
pfaltap=2;
%--------------------------------------------------------------
En=Un/sqrt(3);
U=kE*Un; 
E=kE*En;
if Un==400
   zl1u=0.0346+0.3208i; zl2u=zl1u; zl0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i;
elseif Un==220
   zl1u=0.0829+0.4135i; zl2u=zl1u; zl0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i;
elseif Un==132
   zl1u=0.1342+0.4042i; zl2u=zl1u; zl0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i;
elseif Un==66
   zl1u=0.2156+0.3952i; zl2u=zl1u; zl0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i; 
end
k=(zl0u-zl1u)/(3*zl1u);
zs1=(Un^2/Sccs)*i;
zs2=zs1;
zx1=(Un^2/Sccx)*i;
zx2=zx1;
zr1=(Un^2/Sccr)*i;
zr2=zr1;
zy1=(Un^2/Sccy)*i;
zy2=zy1;
zx0=Kx0*zx1; zr0=Kr0*zr1; zy0=Ky0*zy1;
zs0=k0*zr0; 
Km=zlM0u/zl1u;

switch(tramo)
   
case{1}, % tramo 1: Extremos de Linea SX,SY
Lma=L1+L4;
Lpa=L1+L2+L5;
zLma=abs(zl1u*Lma);
zLpa=abs(zl1u*Lpa);

Lt=L1;
Latm=0;
Lttm=L4;
Latp=0;
Lttp=L2+L5;

em1='Extremo S';
em2='Extremo X';
ep1='Extremo S';
ep2='Extremo Y';

m1=sprintf('Sm');
m2=sprintf('Xm');
p1=sprintf('Sp');
p2=sprintf('Yp');

case{2}, % tramo 2: Extremos de Linea XR,SY
Lma=L4+L2+L3;
Lpa=L1+L2+L5;
zLma=abs(zl1u*Lma);
zLpa=abs(zl1u*Lpa);

Lt=L2;
Latm=L4;
Lttm=L3;
Latp=L1;
Lttp=L5;

em1='Extremo X';
em2='Extremo R';
ep1='Extremo S';
ep2='Extremo Y';

m1=sprintf('Xm');
m2=sprintf('Rm');
p1=sprintf('Sp');
p2=sprintf('Yp');

case{3}, % tramo 3: Extremos de Linea XR,YR
Lma=L4+L2+L3;
Lpa=L3+L5;
zLma=abs(zl1u*Lma);
zLpa=abs(zl1u*Lpa);

Lt=L3;
Latm=L4+L2;
Lttm=0;
Latp=L5;
Lttp=0;

em1='Extremo X';
em2='Extremo R';
ep1='Extremo Y';
ep2='Extremo R';

m1=sprintf('Xm');
m2=sprintf('Rm');
p1=sprintf('Yp');
p2=sprintf('Rp');

case{4}, % tramo 4: Extremos de Linea SX,RX
Lma=L1+L4;
Lpa=L4+L2+L3;
zLma=abs(zl1u*Lma);
zLpa=abs(zl1u*Lpa);

Lt=L4;
Latm=L1;
Lttm=0;
Latp=L3+L2;
Lttp=0;

em1='Extremo S';
em2='Extremo X';
ep1='Extremo R';
ep2='Extremo X';

m1=sprintf('Sm');
m2=sprintf('Xm');
p1=sprintf('Rp');
p2=sprintf('Xp');

case{5}, % tramo 5:  SY,RY
Lma=L1+L2+L5;
Lpa=L3+L5;
zLma=abs(zl1u*Lma);
zLpa=abs(zl1u*Lpa);

Lt=L5;
Latm=L1+L2;
Lttm=0;
Latp=L3;
Lttp=0;

em1='Extremo S';
em2='Extremo Y';
ep1='Extremo R';
ep2='Extremo Y';

m1=sprintf('Sm');
m2=sprintf('Ym');
p1=sprintf('Rp');
p2=sprintf('Yp');

end

if faltam==1 cfaltam='ag'; elementom='AG';elseif faltam==2 cfaltam='bg'; elementom='BG';
elseif faltam==3 cfaltam='cg'; elementom='CG';elseif faltam==4 cfaltam='ab'; elementom='AB';
elseif faltam==5 cfaltam='bc'; elementom='BC';elseif faltam==6 cfaltam='ca'; elementom='CA';
elseif faltam==7 cfaltam='abg'; elementom='AB';elseif faltam==8 cfaltam='bcg'; elementom='BC';
elseif faltam==9 cfaltam='cag'; elementom='CA';elseif faltam==10 cfaltam='abc'; elementom='BC';
elseif faltam==11 cfaltam='abcg'; elementom='BC';
end

if faltap==1 cfaltap='ag'; elementop='AG';elseif faltap==2 cfaltap='bg'; elementop='BG';
elseif faltap==3 cfaltap='cg'; elementop='CG';elseif faltap==4 cfaltap='ab'; elementop='AB';
elseif faltap==5 cfaltap='bc'; elementop='BC';elseif faltap==6 cfaltap='ca'; elementop='CA';
elseif faltap==7 cfaltap='abg'; elementop='AB';elseif faltap==8 cfaltap='bcg'; elementop='BC';
elseif faltap==9 cfaltap='cag'; elementop='CA';elseif faltap==10 cfaltap='abc'; elementop='BC';
elseif faltap==11 cfaltap='abcg'; elementop='BC';
end

fprintf('\n');

fila=1;
i=1;
m=0.02;
while m<1
p=m;

switch(tramo)
   
case{1}, % tramo 1 
[Minc]=subMincT1(1);
[z0]=subz0T1(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,m,p);
[z1]=subz1T1(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,m,p);
[z2]=subz2T1(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,m,p);
case{2}, % tramo 2
[Minc]=subMincT2(1);
[z0]=subz0T2(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,m,p);
[z1]=subz1T2(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,m,p);
[z2]=subz2T2(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,m,p);
case{3}, % tramo 3
[Minc]=subMincT3(1);
[z0]=subz0T3(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,m,p);
[z1]=subz1T3(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,m,p);
[z2]=subz2T3(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,m,p);
case{4}, % tramo 4
[Minc]=subMincT4(1);
[z0]=subz0T4(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,m,p);
[z1]=subz1T4(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,m,p);
[z2]=subz2T4(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,m,p);
case{5}, % tramo 5
[Minc]=subMincT5(1);
[z0]=subz0T5(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,m,p);
[z1]=subz1T5(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,m,p);
[z2]=subz2T5(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,m,p);
end

% YBUS(0,1,2)
[Y0]=subY0(Minc,z0);
[Y1]=subY1(Minc,z1);
[Y2]=subY2(Minc,z2);

% ZBUS(0,1,2)
Z0=inv(Y0);
Z1=inv(Y1);
Z2=inv(Y2);

% CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2,Ip0,Ip1,Ip2,Vp0,Vp1,Vp2]=...
   mag4mp(Z0,Z1,Z2,E,faltam,faltap,Rfm,Rgm,Rfp,Rgp,pfaltam,pfaltap);

% IBUS(0,1,2)
I0=zeros(6,1);
I1=zeros(6,1);
I2=zeros(6,1);

I0(pfaltam,1)=-Im0; 
I1(pfaltam,1)=-Im1; 
I2(pfaltam,1)=-Im2; 

I0(pfaltap,1)=-Ip0;
I1(pfaltap,1)=-Ip1;
I2(pfaltap,1)=-Ip2;

% VBUS(0,1,2)
V0=Z0*I0;
V1=Z1*I1+E;
V2=Z2*I2;

% Ie(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc,V0,z0);
[Ie1]=subIe1(Minc,V1,z1);
[Ie2]=subIe2(Minc,V2,z2);

% V(a,b,c), I(a,b,c) EN LOS PUNTOS m,p DE FALTA
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
[Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2);

% ASIGNACIONES V(0,1,2) EN BARRAS
Vs0=V0(3);Vx0=V0(4);Vr0=V0(5);Vy0=V0(6);
Vs1=V1(3);Vx1=V1(4);Vr1=V1(5);Vy1=V1(6);
Vs2=V2(3);Vx2=V2(4);Vr2=V2(5);Vy2=V2(6);

% V(a,b,c) EN LAS BARRAS
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vxa,Vxb,Vxc]=abc(Vx0,Vx1,Vx2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Vya,Vyb,Vyc]=abc(Vy0,Vy1,Vy2);

switch(tramo)
   
case{1}, %tramo1
   
%case{1}, % SX,SY
%Lma=L1+L4;
%Lpa=L1+L2+L5;

% ASIGNACIONES I(0,1,2) EN LOS ELEMENTOS
Ism0=Ie0(5);Ixm0=Ie0(7);Isp0=Ie0(6);Iyp0=Ie0(10);
Ism1=Ie1(5);Ixm1=Ie1(7);Isp1=Ie1(6);Iyp1=Ie1(10);
Ism2=Ie2(5);Ixm2=Ie2(7);Isp2=Ie2(6);Iyp2=Ie2(10);
Ixr0=Ie0(8);
Iyr0=-Ie0(9);

% I(a,b,c) EN LOS ELEMENTOS
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2);
[Ixma,Ixmb,Ixmc]=abc(Ixm0,Ixm1,Ixm2);
[Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2);
[Iypa,Iypb,Iypc]=abc(Iyp0,Iyp1,Iyp2);

switch(compens21)
case{1},
% COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
  PDistComp(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k,Km,Isp0);
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
   PDistComp(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k,Km,Ixr0);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
   PDistComp(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k,Km,Ism0);
[Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG]=...
  PDistComp(Vya,Vyb,Vyc,Iypa,Iypb,Iypc,Iyp0,k,Km,Iyr0);

case{2},
% NO COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
  PDist(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k);
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
   PDist(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
   PDist(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k);
[Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG]=...
  PDist(Vya,Vyb,Vyc,Iypa,Iypb,Iypc,Iyp0,k);
end

% MEDICIONES 67N 
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0);
[Vxpol,Ixmop,AngDirXm]=PDire(Vx0,Ixm0);
[Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0);
[Vypol,Iypop,AngDirYp]=PDire(Vy0,Iyp0);

% Recuerdo: tramo1 SX,SY

%case{1}, % SX,SY
%Lma=L1+L4;
%Lpa=L1+L2+L5;

zLSm=abs(zl1u*m*L1);
zLXm=abs(zl1u*(L4+(1-m)*L1));
zLSp=abs(zl1u*p*L1);
zLYp=abs(zl1u*(L5+L2+(1-p)*L1));

[alcanceSm_ZABpu,alcanceSm_ZBCpu,alcanceSm_ZCApu,alcanceSm_ZAGpu,alcanceSm_ZBGpu,alcanceSm_ZCGpu]=...
   subalcance(Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG,zl1u,zLSm);
OrdenadaSm_ZABpu(1,i)=alcanceSm_ZABpu;
OrdenadaSm_ZBCpu(1,i)=alcanceSm_ZBCpu;
OrdenadaSm_ZCApu(1,i)=alcanceSm_ZCApu;
OrdenadaSm_ZAGpu(1,i)=alcanceSm_ZAGpu;
OrdenadaSm_ZBGpu(1,i)=alcanceSm_ZBGpu;
OrdenadaSm_ZCGpu(1,i)=alcanceSm_ZCGpu;

[alcanceXm_ZABpu,alcanceXm_ZBCpu,alcanceXm_ZCApu,alcanceXm_ZAGpu,alcanceXm_ZBGpu,alcanceXm_ZCGpu]=...
subalcance(Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG,zl1u,zLXm);
OrdenadaXm_ZABpu(1,i)=alcanceXm_ZABpu;
OrdenadaXm_ZBCpu(1,i)=alcanceXm_ZBCpu;
OrdenadaXm_ZCApu(1,i)=alcanceXm_ZCApu;
OrdenadaXm_ZAGpu(1,i)=alcanceXm_ZAGpu;
OrdenadaXm_ZBGpu(1,i)=alcanceXm_ZBGpu;

OrdenadaXm_ZCGpu(1,i)=alcanceXm_ZCGpu;

[alcanceSp_ZABpu,alcanceSp_ZBCpu,alcanceSp_ZCApu,alcanceSp_ZAGpu,alcanceSp_ZBGpu,alcanceSp_ZCGpu]=...
subalcance(Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG,zl1u,zLSp);
OrdenadaSp_ZABpu(1,i)=alcanceSp_ZABpu;
OrdenadaSp_ZBCpu(1,i)=alcanceSp_ZBCpu;
OrdenadaSp_ZCApu(1,i)=alcanceSp_ZCApu;
OrdenadaSp_ZAGpu(1,i)=alcanceSp_ZAGpu;
OrdenadaSp_ZBGpu(1,i)=alcanceSp_ZBGpu;
OrdenadaSp_ZCGpu(1,i)=alcanceSp_ZCGpu;

[alcanceYp_ZABpu,alcanceYp_ZBCpu,alcanceYp_ZCApu,alcanceYp_ZAGpu,alcanceYp_ZBGpu,alcanceYp_ZCGpu]=...
subalcance(Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG,zl1u,zLYp);
OrdenadaYp_ZABpu(1,i)=alcanceYp_ZABpu;
OrdenadaYp_ZBCpu(1,i)=alcanceYp_ZBCpu;
OrdenadaYp_ZCApu(1,i)=alcanceYp_ZCApu;
OrdenadaYp_ZAGpu(1,i)=alcanceYp_ZAGpu;
OrdenadaYp_ZBGpu(1,i)=alcanceYp_ZBGpu;
OrdenadaYp_ZCGpu(1,i)=alcanceYp_ZCGpu;

Abscisa(fila,i)=m;

[alcanceSm_ZABpua,alcanceSm_ZBCpua,alcanceSm_ZCApua,alcanceSm_ZAGpua,alcanceSm_ZBGpua,alcanceSm_ZCGpua]=...
subalcance(Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG,zl1u,zLma);
OrdenadaSm_ZABpua(1,i)=alcanceSm_ZABpua;
OrdenadaSm_ZBCpua(1,i)=alcanceSm_ZBCpua;
OrdenadaSm_ZCApua(1,i)=alcanceSm_ZCApua;
OrdenadaSm_ZAGpua(1,i)=alcanceSm_ZAGpua;
OrdenadaSm_ZBGpua(1,i)=alcanceSm_ZBGpua;
OrdenadaSm_ZCGpua(1,i)=alcanceSm_ZCGpua;

[alcanceXm_ZABpua,alcanceXm_ZBCpua,alcanceXm_ZCApua,alcanceXm_ZAGpua,alcanceXm_ZBGpua,alcanceXm_ZCGpua]=...
subalcance(Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG,zl1u,zLma);
OrdenadaXm_ZABpua(1,i)=alcanceXm_ZABpua;
OrdenadaXm_ZBCpua(1,i)=alcanceXm_ZBCpua;
OrdenadaXm_ZCApua(1,i)=alcanceXm_ZCApua;
OrdenadaXm_ZAGpua(1,i)=alcanceXm_ZAGpua;
OrdenadaXm_ZBGpua(1,i)=alcanceXm_ZBGpua;
OrdenadaXm_ZCGpua(1,i)=alcanceXm_ZCGpua;

[alcanceSp_ZABpua,alcanceSp_ZBCpua,alcanceSp_ZCApua,alcanceSp_ZAGpua,alcanceSp_ZBGpua,alcanceSp_ZCGpua]=...
subalcance(Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG,zl1u,zLpa);
OrdenadaSp_ZABpua(1,i)=alcanceSp_ZABpua;
OrdenadaSp_ZBCpua(1,i)=alcanceSp_ZBCpua;
OrdenadaSp_ZCApua(1,i)=alcanceSp_ZCApua;
OrdenadaSp_ZAGpua(1,i)=alcanceSp_ZAGpua;
OrdenadaSp_ZBGpua(1,i)=alcanceSp_ZBGpua;
OrdenadaSp_ZCGpua(1,i)=alcanceSp_ZCGpua;

[alcanceYp_ZABpua,alcanceYp_ZBCpua,alcanceYp_ZCApua,alcanceYp_ZAGpua,alcanceYp_ZBGpua,alcanceYp_ZCGpua]=...
subalcance(Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG,zl1u,zLpa);
OrdenadaYp_ZABpua(1,i)=alcanceYp_ZABpua;
OrdenadaYp_ZBCpua(1,i)=alcanceYp_ZBCpua;
OrdenadaYp_ZCApua(1,i)=alcanceYp_ZCApua;
OrdenadaYp_ZAGpua(1,i)=alcanceYp_ZAGpua;
OrdenadaYp_ZBGpua(1,i)=alcanceYp_ZBGpua;
OrdenadaYp_ZCGpua(1,i)=alcanceYp_ZCGpua;

%***************************************
%Para salida 67N
%***************************************
Vpol_m_izq=Vspol;
Iop_m_izq=Ismop;
AngDir_m_izq=AngDirSm;

Vpol_m_der=Vxpol;
Iop_m_der=Ixmop;
AngDir_m_der=AngDirXm;

cadVpol_m_izq='modulo Vspol(kV)';
cadIop_m_izq='modulo Ismop(kA)';
cadAngDir_m_izq='AngDirSm(º)';

cadVpol_m_der='modulo Vxpol(kV)';
cadIop_m_der='modulo Ixmop(kV)';
cadAngDir_m_der='AngDirXm(º)';

Vpol_p_izq=Vspol;
Iop_p_izq=Ispop;
AngDir_p_izq=AngDirSp;

Vpol_p_der=Vypol;
Iop_p_der=Iypop;
AngDir_p_der=AngDirYp;

cadVpol_p_izq='modulo Vspol(kV)';
cadIop_p_izq='modulo Ispop(kA)';
cadAngDir_p_izq='AngDirSp(º)';

cadVpol_p_der='modulo Vypol(kV)';
cadIop_p_der='modulo Iypop(kV)';
cadAngDir_p_der='AngDirYp(º)';

%***************************************
% Para Salida 4 R-X (Tramo 1)
%***************************************
ZAB_m1=Sm_ZAB;
ZBC_m1=Sm_ZBC;
ZCA_m1=Sm_ZCA;
ZAG_m1=Sm_ZAG;
ZBG_m1=Sm_ZBG;
ZCG_m1=Sm_ZCG;

ZAB_m2=Xm_ZAB;
ZBC_m2=Xm_ZBC;
ZCA_m2=Xm_ZCA;
ZAG_m2=Xm_ZAG;
ZBG_m2=Xm_ZBG;
ZCG_m2=Xm_ZCG;

ZAB_p1=Sp_ZAB;
ZBC_p1=Sp_ZBC;
ZCA_p1=Sp_ZCA;
ZAG_p1=Sp_ZAG;
ZBG_p1=Sp_ZBG;
ZCG_p1=Sp_ZCG;

ZAB_p2=Yp_ZAB;
ZBC_p2=Yp_ZBC;
ZCA_p2=Yp_ZCA;
ZAG_p2=Yp_ZAG;
ZBG_p2=Yp_ZBG;
ZCG_p2=Yp_ZCG;

%***************************************

case{2},% tramo2
   
%case{2}, % XR,SY
%Lma=L4+L2+L3;
%Lpa=L1+L2+L5;

Ixm0=Ie0(7);Irm0=Ie0(8);Isp0=Ie0(6);Iyp0=Ie0(10);
Ixm1=Ie1(7);Irm1=Ie1(8);Isp1=Ie1(6);Iyp1=Ie1(10);
Ixm2=Ie2(7);Irm2=Ie2(8);Isp2=Ie2(6);Iyp2=Ie2(10);
Ixs0=-Ie0(5);
Iry0=Ie0(9);
Isx0=Ie0(5);
Iyr0=-Ie0(9);

[Ixma,Ixmb,Ixmc]=abc(Ixm0,Ixm1,Ixm2);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
[Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2);
[Iypa,Iypb,Iypc]=abc(Iyp0,Iyp1,Iyp2);

switch(compens21)
case{1},
% COMPENSADA
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
  PDistComp(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k,Km,Ixs0);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
   PDistComp(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k,Km,Iry0);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
   PDistComp(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k,Km,Isx0);
[Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG]=...
  PDistComp(Vya,Vyb,Vyc,Iypa,Iypb,Iypc,Iyp0,k,Km,Iyr0);
case{2},
% NO COMPENSADA
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
  PDist(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
   PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
   PDist(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k);
[Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG]=...
  PDist(Vya,Vyb,Vyc,Iypa,Iypb,Iypc,Iyp0,k);
end

% MEDICION 21 EN LOS EXTREMOS DE LINEA
%
% MEDICION 67N EN LOS EXTREMOS DE LINEA
[Vxpol,Ixmop,AngDirXm]=PDire(Vx0,Ixm0);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);
[Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0);
[Vypol,Iypop,AngDirYp]=PDire(Vy0,Iyp0);

% Recuerdo: tramo 2

zLXm=abs(zl1u*(L4+m*L2));
zLRm=abs(zl1u*(L3+(1-m)*L2));
zLSp=abs(zl1u*(L1+p*L2));
zLYp=abs(zl1u*(L5+(1-p)*L2));

[alcanceXm_ZABpu,alcanceXm_ZBCpu,alcanceXm_ZCApu,alcanceXm_ZAGpu,alcanceXm_ZBGpu,alcanceXm_ZCGpu]=...
   subalcance(Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG,zl1u,zLXm);
OrdenadaXm_ZABpu(1,i)=alcanceXm_ZABpu;
OrdenadaXm_ZBCpu(1,i)=alcanceXm_ZBCpu;
OrdenadaXm_ZCApu(1,i)=alcanceXm_ZCApu;
OrdenadaXm_ZAGpu(1,i)=alcanceXm_ZAGpu;
OrdenadaXm_ZBGpu(1,i)=alcanceXm_ZBGpu;

OrdenadaXm_ZCGpu(1,i)=alcanceXm_ZCGpu;

[alcanceRm_ZABpu,alcanceRm_ZBCpu,alcanceRm_ZCApu,alcanceRm_ZAGpu,alcanceRm_ZBGpu,alcanceRm_ZCGpu]=...
   subalcance(Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG,zl1u,zLRm);
OrdenadaRm_ZABpu(1,i)=alcanceRm_ZABpu;
OrdenadaRm_ZBCpu(1,i)=alcanceRm_ZBCpu;
OrdenadaRm_ZCApu(1,i)=alcanceRm_ZCApu;
OrdenadaRm_ZAGpu(1,i)=alcanceRm_ZAGpu;
OrdenadaRm_ZBGpu(1,i)=alcanceRm_ZBGpu;
OrdenadaRm_ZCGpu(1,i)=alcanceRm_ZCGpu;

[alcanceSp_ZABpu,alcanceSp_ZBCpu,alcanceSp_ZCApu,alcanceSp_ZAGpu,alcanceSp_ZBGpu,alcanceSp_ZCGpu]=...
   subalcance(Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG,zl1u,zLSp);
OrdenadaSp_ZABpu(1,i)=alcanceSp_ZABpu;
OrdenadaSp_ZBCpu(1,i)=alcanceSp_ZBCpu;
OrdenadaSp_ZCApu(1,i)=alcanceSp_ZCApu;
OrdenadaSp_ZAGpu(1,i)=alcanceSp_ZAGpu;
OrdenadaSp_ZBGpu(1,i)=alcanceSp_ZBGpu;
OrdenadaSp_ZCGpu(1,i)=alcanceSp_ZCGpu;

[alcanceYp_ZABpu,alcanceYp_ZBCpu,alcanceYp_ZCApu,alcanceYp_ZAGpu,alcanceYp_ZBGpu,alcanceYp_ZCGpu]=...
   subalcance(Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG,zl1u,zLYp);
OrdenadaYp_ZABpu(1,i)=alcanceYp_ZABpu;
OrdenadaYp_ZBCpu(1,i)=alcanceYp_ZBCpu;
OrdenadaYp_ZCApu(1,i)=alcanceYp_ZCApu;
OrdenadaYp_ZAGpu(1,i)=alcanceYp_ZAGpu;
OrdenadaYp_ZBGpu(1,i)=alcanceYp_ZBGpu;
OrdenadaYp_ZCGpu(1,i)=alcanceYp_ZCGpu;

Abscisa(fila,i)=m;

[alcanceXm_ZABpua,alcanceXm_ZBCpua,alcanceXm_ZCApua,alcanceXm_ZAGpua,alcanceXm_ZBGpua,alcanceXm_ZCGpua]=...
   subalcance(Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG,zl1u,zLma);
OrdenadaXm_ZABpua(1,i)=alcanceXm_ZABpua;
OrdenadaXm_ZBCpua(1,i)=alcanceXm_ZBCpua;
OrdenadaXm_ZCApua(1,i)=alcanceXm_ZCApua;
OrdenadaXm_ZAGpua(1,i)=alcanceXm_ZAGpua;
OrdenadaXm_ZBGpua(1,i)=alcanceXm_ZBGpua;
OrdenadaXm_ZCGpua(1,i)=alcanceXm_ZCGpua;

[alcanceRm_ZABpua,alcanceRm_ZBCpua,alcanceRm_ZCApua,alcanceRm_ZAGpua,alcanceRm_ZBGpua,alcanceRm_ZCGpua]=...
   subalcance(Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG,zl1u,zLma);
OrdenadaRm_ZABpua(1,i)=alcanceRm_ZABpua;
OrdenadaRm_ZBCpua(1,i)=alcanceRm_ZBCpua;
OrdenadaRm_ZCApua(1,i)=alcanceRm_ZCApua;
OrdenadaRm_ZAGpua(1,i)=alcanceRm_ZAGpua;
OrdenadaRm_ZBGpua(1,i)=alcanceRm_ZBGpua;
OrdenadaRm_ZCGpua(1,i)=alcanceRm_ZCGpua;

[alcanceSp_ZABpua,alcanceSp_ZBCpua,alcanceSp_ZCApua,alcanceSp_ZAGpua,alcanceSp_ZBGpua,alcanceSp_ZCGpua]=...
   subalcance(Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG,zl1u,zLpa);
OrdenadaSp_ZABpua(1,i)=alcanceSp_ZABpua;
OrdenadaSp_ZBCpua(1,i)=alcanceSp_ZBCpua;
OrdenadaSp_ZCApua(1,i)=alcanceSp_ZCApua;
OrdenadaSp_ZAGpua(1,i)=alcanceSp_ZAGpua;
OrdenadaSp_ZBGpua(1,i)=alcanceSp_ZBGpua;
OrdenadaSp_ZCGpua(1,i)=alcanceSp_ZCGpua;

[alcanceYp_ZABpua,alcanceYp_ZBCpua,alcanceYp_ZCApua,alcanceYp_ZAGpua,alcanceYp_ZBGpua,alcanceYp_ZCGpua]=...
   subalcance(Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG,zl1u,zLpa);
OrdenadaYp_ZABpua(1,i)=alcanceYp_ZABpua;
OrdenadaYp_ZBCpua(1,i)=alcanceYp_ZBCpua;
OrdenadaYp_ZCApua(1,i)=alcanceYp_ZCApua;
OrdenadaYp_ZAGpua(1,i)=alcanceYp_ZAGpua;
OrdenadaYp_ZBGpua(1,i)=alcanceYp_ZBGpua;
OrdenadaYp_ZCGpua(1,i)=alcanceYp_ZCGpua;

%***************************************
%Para salida 67N
%***************************************
Vpol_m_izq=Vxpol;
Iop_m_izq=Ixmop;
AngDir_m_izq=AngDirXm;

Vpol_m_der=Vrpol;
Iop_m_der=Irmop;
AngDir_m_der=AngDirRm;

cadVpol_m_izq='modulo Vxpol(kV)';
cadIop_m_izq='modulo Ixmop(kA)';
cadAngDir_m_izq='AngDirXm(º)';

cadVpol_m_der='modulo Vrpol(kV)';
cadIop_m_der='modulo Irmop(kV)';
cadAngDir_m_der='AngDirRm(º)';

Vpol_p_izq=Vspol;
Iop_p_izq=Ispop;
AngDir_p_izq=AngDirSp;

Vpol_p_der=Vypol;
Iop_p_der=Iypop;
AngDir_p_der=AngDirYp;

cadVpol_p_izq='modulo Vspol(kV)';
cadIop_p_izq='modulo Ispop(kA)';
cadAngDir_p_izq='AngDirSp(º)';

cadVpol_p_der='modulo Vypol(kV)';
cadIop_p_der='modulo Iypop(kV)';
cadAngDir_p_der='AngDirYp(º)';

%***************************************
% Para Salida 4 R-X (Tramo 2)
%***************************************
ZAB_m1=Xm_ZAB;
ZBC_m1=Xm_ZBC;
ZCA_m1=Xm_ZCA;
ZAG_m1=Xm_ZAG;
ZBG_m1=Xm_ZBG;
ZCG_m1=Xm_ZCG;

ZAB_m2=Rm_ZAB;
ZBC_m2=Rm_ZBC;
ZCA_m2=Rm_ZCA;
ZAG_m2=Rm_ZAG;
ZBG_m2=Rm_ZBG;
ZCG_m2=Rm_ZCG;

ZAB_p1=Sp_ZAB;
ZBC_p1=Sp_ZBC;
ZCA_p1=Sp_ZCA;
ZAG_p1=Sp_ZAG;
ZBG_p1=Sp_ZBG;
ZCG_p1=Sp_ZCG;

ZAB_p2=Yp_ZAB;
ZBC_p2=Yp_ZBC;
ZCA_p2=Yp_ZCA;
ZAG_p2=Yp_ZAG;
ZBG_p2=Yp_ZBG;
ZCG_p2=Yp_ZCG;

%***************************************


case{3}, % tramo3
   
%case{3}, % XR,YR
%Lma=L4+L2+L3;
%Lpa=L3+L5;
   
Ixm0=Ie0(7);Irm0=Ie0(8);Iyp0=Ie0(9);Irp0=Ie0(10);
Ixm1=Ie1(7);Irm1=Ie1(8);Iyp1=Ie1(9);Irp1=Ie1(10);
Ixm2=Ie2(7);Irm2=Ie2(8);Iyp2=Ie2(9);Irp2=Ie2(10);
Ixs0=-Ie0(5);
Iys0=-Ie0(6);
   

[Ixma,Ixmb,Ixmc]=abc(Ixm0,Ixm1,Ixm2);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
[Iypa,Iypb,Iypc]=abc(Iyp0,Iyp1,Iyp2);
[Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2);

switch(compens21)
case{1},
% COMPENSADA
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
  PDistComp(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k,Km,Ixs0);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
   PDistComp(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k,Km,Irp0);
[Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG]=...
  PDistComp(Vya,Vyb,Vyc,Iypa,Iypb,Iypc,Iyp0,k,Km,Iys0);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
   PDistComp(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k,Km,Irm0);
   
case{2},
% NO COMPENSADA
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
  PDist(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
   PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
[Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG]=...
  PDist(Vya,Vyb,Vyc,Iypa,Iypb,Iypc,Iyp0,k);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
   PDist(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k);
end
%
% MEDICION 67N EN LOS EXTREMOS DE LINEA
[Vxpol,Ixmop,AngDirXm]=PDire(Vx0,Ixm0);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);
[Vypol,Iypop,AngDirYp]=PDire(Vy0,Iyp0);
[Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0);

% Recuerdo: tramo3

%case{3}, % XR,YR
%Lma=L4+L2+L3;
%Lpa=L3+L5;

zLXm=abs(zl1u*(L4+L2+m*L3));
zLRm=abs(zl1u*(1-m)*L3);
zLYp=abs(zl1u*(L5+p*L3));
zLRp=abs(zl1u*(1-p)*L3);

[alcanceXm_ZABpu,alcanceXm_ZBCpu,alcanceXm_ZCApu,alcanceXm_ZAGpu,alcanceXm_ZBGpu,alcanceXm_ZCGpu]=...
subalcance(Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG,zl1u,zLXm);
OrdenadaXm_ZABpu(1,i)=alcanceXm_ZABpu;
OrdenadaXm_ZBCpu(1,i)=alcanceXm_ZBCpu;
OrdenadaXm_ZCApu(1,i)=alcanceXm_ZCApu;
OrdenadaXm_ZAGpu(1,i)=alcanceXm_ZAGpu;
OrdenadaXm_ZBGpu(1,i)=alcanceXm_ZBGpu;

OrdenadaXm_ZCGpu(1,i)=alcanceXm_ZCGpu;

[alcanceRm_ZABpu,alcanceRm_ZBCpu,alcanceRm_ZCApu,alcanceRm_ZAGpu,alcanceRm_ZBGpu,alcanceRm_ZCGpu]=...
   subalcance(Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG,zl1u,zLRm);
OrdenadaRm_ZABpu(1,i)=alcanceRm_ZABpu;
OrdenadaRm_ZBCpu(1,i)=alcanceRm_ZBCpu;
OrdenadaRm_ZCApu(1,i)=alcanceRm_ZCApu;
OrdenadaRm_ZAGpu(1,i)=alcanceRm_ZAGpu;
OrdenadaRm_ZBGpu(1,i)=alcanceRm_ZBGpu;
OrdenadaRm_ZCGpu(1,i)=alcanceRm_ZCGpu;

[alcanceYp_ZABpu,alcanceYp_ZBCpu,alcanceYp_ZCApu,alcanceYp_ZAGpu,alcanceYp_ZBGpu,alcanceYp_ZCGpu]=...
   subalcance(Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG,zl1u,zLYp);
OrdenadaYp_ZABpu(1,i)=alcanceYp_ZABpu;
OrdenadaYp_ZBCpu(1,i)=alcanceYp_ZBCpu;
OrdenadaYp_ZCApu(1,i)=alcanceYp_ZCApu;
OrdenadaYp_ZAGpu(1,i)=alcanceYp_ZAGpu;
OrdenadaYp_ZBGpu(1,i)=alcanceYp_ZBGpu;
OrdenadaYp_ZCGpu(1,i)=alcanceYp_ZCGpu;

[alcanceRp_ZABpu,alcanceRp_ZBCpu,alcanceRp_ZCApu,alcanceRp_ZAGpu,alcanceRp_ZBGpu,alcanceRp_ZCGpu]=...
   subalcance(Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG,zl1u,zLRp);
OrdenadaRp_ZABpu(1,i)=alcanceRp_ZABpu;
OrdenadaRp_ZBCpu(1,i)=alcanceRp_ZBCpu;
OrdenadaRp_ZCApu(1,i)=alcanceRp_ZCApu;
OrdenadaRp_ZAGpu(1,i)=alcanceRp_ZAGpu;
OrdenadaRp_ZBGpu(1,i)=alcanceRp_ZBGpu;
OrdenadaRp_ZCGpu(1,i)=alcanceRp_ZCGpu;


Abscisa(fila,i)=m;


[alcanceXm_ZABpua,alcanceXm_ZBCpua,alcanceXm_ZCApua,alcanceXm_ZAGpua,alcanceXm_ZBGpua,alcanceXm_ZCGpua]=...
   subalcance(Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG,zl1u,zLma);
OrdenadaXm_ZABpua(1,i)=alcanceXm_ZABpua;
OrdenadaXm_ZBCpua(1,i)=alcanceXm_ZBCpua;
OrdenadaXm_ZCApua(1,i)=alcanceXm_ZCApua;
OrdenadaXm_ZAGpua(1,i)=alcanceXm_ZAGpua;
OrdenadaXm_ZBGpua(1,i)=alcanceXm_ZBGpua;

OrdenadaXm_ZCGpua(1,i)=alcanceXm_ZCGpua;

[alcanceRm_ZABpua,alcanceRm_ZBCpua,alcanceRm_ZCApua,alcanceRm_ZAGpua,alcanceRm_ZBGpua,alcanceRm_ZCGpua]=...
   subalcance(Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG,zl1u,zLma);
OrdenadaRm_ZABpua(1,i)=alcanceRm_ZABpua;
OrdenadaRm_ZBCpua(1,i)=alcanceRm_ZBCpua;
OrdenadaRm_ZCApua(1,i)=alcanceRm_ZCApua;
OrdenadaRm_ZAGpua(1,i)=alcanceRm_ZAGpua;
OrdenadaRm_ZBGpua(1,i)=alcanceRm_ZBGpua;
OrdenadaRm_ZCGpua(1,i)=alcanceRm_ZCGpua;

[alcanceYp_ZABpua,alcanceYp_ZBCpua,alcanceYp_ZCApua,alcanceYp_ZAGpua,alcanceYp_ZBGpua,alcanceYp_ZCGpua]=...
   subalcance(Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG,zl1u,zLpa);
OrdenadaYp_ZABpua(1,i)=alcanceYp_ZABpua;
OrdenadaYp_ZBCpua(1,i)=alcanceYp_ZBCpua;
OrdenadaYp_ZCApua(1,i)=alcanceYp_ZCApua;
OrdenadaYp_ZAGpua(1,i)=alcanceYp_ZAGpua;
OrdenadaYp_ZBGpua(1,i)=alcanceYp_ZBGpua;
OrdenadaYp_ZCGpua(1,i)=alcanceYp_ZCGpua;

[alcanceRp_ZABpua,alcanceRp_ZBCpua,alcanceRp_ZCApua,alcanceRp_ZAGpua,alcanceRp_ZBGpua,alcanceRp_ZCGpua]=...
   subalcance(Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG,zl1u,zLpa);
OrdenadaRp_ZABpua(1,i)=alcanceRp_ZABpua;
OrdenadaRp_ZBCpua(1,i)=alcanceRp_ZBCpua;
OrdenadaRp_ZCApua(1,i)=alcanceRp_ZCApua;
OrdenadaRp_ZAGpua(1,i)=alcanceRp_ZAGpua;
OrdenadaRp_ZBGpua(1,i)=alcanceRp_ZBGpua;
OrdenadaRp_ZCGpua(1,i)=alcanceRp_ZCGpua;

%***************************************
%Para salida 67N
%***************************************
Vpol_m_izq=Vxpol;
Iop_m_izq=Ixmop;
AngDir_m_izq=AngDirXm;

Vpol_m_der=Vrpol;
Iop_m_der=Irmop;
AngDir_m_der=AngDirRm;

cadVpol_m_izq='modulo Vxpol(kV)';
cadIop_m_izq='modulo Ixmop(kA)';
cadAngDir_m_izq='AngDirXm(º)';

cadVpol_m_der='modulo Vrpol(kV)';
cadIop_m_der='modulo Irmop(kV)';
cadAngDir_m_der='AngDirRm(º)';

Vpol_p_izq=Vypol;
Iop_p_izq=Iypop;
AngDir_p_izq=AngDirYp;

Vpol_p_der=Vrpol;
Iop_p_der=Irpop;
AngDir_p_der=AngDirRp;

cadVpol_p_izq='modulo Vypol(kV)';
cadIop_p_izq='modulo Iypop(kA)';
cadAngDir_p_izq='AngDirYp(º)';

cadVpol_p_der='modulo Vrpol(kV)';
cadIop_p_der='modulo Irpop(kV)';
cadAngDir_p_der='AngDirRp(º)';

%***************************************
% Para Salida 4  R-X (Tramo 3)
%***************************************
ZAB_m1=Xm_ZAB;
ZBC_m1=Xm_ZBC;
ZCA_m1=Xm_ZCA;
ZAG_m1=Xm_ZAG;
ZBG_m1=Xm_ZBG;
ZCG_m1=Xm_ZCG;

ZAB_m2=Rm_ZAB;
ZBC_m2=Rm_ZBC;
ZCA_m2=Rm_ZCA;
ZAG_m2=Rm_ZAG;
ZBG_m2=Rm_ZBG;
ZCG_m2=Rm_ZCG;

ZAB_p1=Yp_ZAB;
ZBC_p1=Yp_ZBC;
ZCA_p1=Yp_ZCA;
ZAG_p1=Yp_ZAG;
ZBG_p1=Yp_ZBG;
ZCG_p1=Yp_ZCG;

ZAB_p2=Rp_ZAB;
ZBC_p2=Rp_ZBC;
ZCA_p2=Rp_ZCA;
ZAG_p2=Rp_ZAG;
ZBG_p2=Rp_ZBG;
ZCG_p2=Rp_ZCG;

%***************************************

case{4}, %tramo4
   
%case{4}, % SX,RX
%Lma=L1+L4;
%Lpa=L4+L2+L3;
   
Ism0=Ie0(5);Ixm0=Ie0(7);Irp0=Ie0(9);Ixp0=Ie0(10);
Ism1=Ie1(5);Ixm1=Ie1(7);Irp1=Ie1(9);Ixp1=Ie1(10);
Ism2=Ie2(5);Ixm2=Ie2(7);Irp2=Ie2(9);Ixp2=Ie2(10);
Isy0=Ie0(6);
Iry0=-Ie0(8);
   
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2);
[Ixma,Ixmb,Ixmc]=abc(Ixm0,Ixm1,Ixm2);
[Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2);
[Ixpa,Ixpb,Ixpc]=abc(Ixp0,Ixp1,Ixp2);

switch(compens21)
case{1},
% COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
  PDistComp(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k,Km,Isy0);
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
   PDistComp(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k,Km,Ixp0);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
   PDistComp(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k,Km,Iry0);
[Xp_ZAB,Xp_ZBC,Xp_ZCA,Xp_ZAG,Xp_ZBG,Xp_ZCG]=...
   PDistComp(Vxa,Vxb,Vxc,Ixpa,Ixpb,Ixpc,Ixp0,k,Km,Ixm0);
   
case{2},
% NO COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
  PDist(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k);
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
   PDist(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
   PDist(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k);
[Xp_ZAB,Xp_ZBC,Xp_ZCA,Xp_ZAG,Xp_ZBG,Xp_ZCG]=...
   PDist(Vxa,Vxb,Vxc,Ixpa,Ixpb,Ixpc,Ixp0,k);
end


% MEDICION 21 EN LOS EXTREMOS DE LINEA
%



% MEDICION 67N EN LOS EXTREMOS DE LINEA
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0);
[Vxpol,Ixmop,AngDirXm]=PDire(Vx0,Ixm0);
[Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0);
[Vxpol,Ixpop,AngDirXp]=PDire(Vx0,Ixp0);

% Recuerdo: tramo4

%case{4}, % SX,RX
%Lma=L1+L4;
%Lpa=L4+L2+L3;

zLSm=abs(zl1u*(L1+m*L4));
zLXm=abs(zl1u*(1-m)*L4);
zLRp=abs(zl1u*(L3+L2+p*L4));
zLXp=abs(zl1u*(1-p)*L4);

[alcanceSm_ZABpu,alcanceSm_ZBCpu,alcanceSm_ZCApu,alcanceSm_ZAGpu,alcanceSm_ZBGpu,alcanceSm_ZCGpu]=...
   subalcance(Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG,zl1u,zLSm);
OrdenadaSm_ZABpu(1,i)=alcanceSm_ZABpu;
OrdenadaSm_ZBCpu(1,i)=alcanceSm_ZBCpu;
OrdenadaSm_ZCApu(1,i)=alcanceSm_ZCApu;
OrdenadaSm_ZAGpu(1,i)=alcanceSm_ZAGpu;
OrdenadaSm_ZBGpu(1,i)=alcanceSm_ZBGpu;
OrdenadaSm_ZCGpu(1,i)=alcanceSm_ZCGpu;

[alcanceXm_ZABpu,alcanceXm_ZBCpu,alcanceXm_ZCApu,alcanceXm_ZAGpu,alcanceXm_ZBGpu,alcanceXm_ZCGpu]=...
subalcance(Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG,zl1u,zLXm);
OrdenadaXm_ZABpu(1,i)=alcanceXm_ZABpu;
OrdenadaXm_ZBCpu(1,i)=alcanceXm_ZBCpu;
OrdenadaXm_ZCApu(1,i)=alcanceXm_ZCApu;
OrdenadaXm_ZAGpu(1,i)=alcanceXm_ZAGpu;
OrdenadaXm_ZBGpu(1,i)=alcanceXm_ZBGpu;

OrdenadaXm_ZCGpu(1,i)=alcanceXm_ZCGpu;

[alcanceRp_ZABpu,alcanceRp_ZBCpu,alcanceRp_ZCApu,alcanceRp_ZAGpu,alcanceRp_ZBGpu,alcanceRp_ZCGpu]=...
subalcance(Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG,zl1u,zLRp);
OrdenadaRp_ZABpu(1,i)=alcanceRp_ZABpu;
OrdenadaRp_ZBCpu(1,i)=alcanceRp_ZBCpu;
OrdenadaRp_ZCApu(1,i)=alcanceRp_ZCApu;
OrdenadaRp_ZAGpu(1,i)=alcanceRp_ZAGpu;
OrdenadaRp_ZBGpu(1,i)=alcanceRp_ZBGpu;
OrdenadaRp_ZCGpu(1,i)=alcanceRp_ZCGpu;

[alcanceXp_ZABpu,alcanceXp_ZBCpu,alcanceXp_ZCApu,alcanceXp_ZAGpu,alcanceXp_ZBGpu,alcanceXp_ZCGpu]=...
subalcance(Xp_ZAB,Xp_ZBC,Xp_ZCA,Xp_ZAG,Xp_ZBG,Xp_ZCG,zl1u,zLXp);
OrdenadaXp_ZABpu(1,i)=alcanceXp_ZABpu;
OrdenadaXp_ZBCpu(1,i)=alcanceXp_ZBCpu;
OrdenadaXp_ZCApu(1,i)=alcanceXp_ZCApu;
OrdenadaXp_ZAGpu(1,i)=alcanceXp_ZAGpu;
OrdenadaXp_ZBGpu(1,i)=alcanceXp_ZBGpu;
OrdenadaXp_ZCGpu(1,i)=alcanceXp_ZCGpu;

Abscisa(fila,i)=m;

[alcanceSm_ZABpua,alcanceSm_ZBCpua,alcanceSm_ZCApua,alcanceSm_ZAGpua,alcanceSm_ZBGpua,alcanceSm_ZCGpua]=...
subalcance(Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG,zl1u,zLma);
OrdenadaSm_ZABpua(1,i)=alcanceSm_ZABpua;
OrdenadaSm_ZBCpua(1,i)=alcanceSm_ZBCpua;
OrdenadaSm_ZCApua(1,i)=alcanceSm_ZCApua;
OrdenadaSm_ZAGpua(1,i)=alcanceSm_ZAGpua;
OrdenadaSm_ZBGpua(1,i)=alcanceSm_ZBGpua;
OrdenadaSm_ZCGpua(1,i)=alcanceSm_ZCGpua;

[alcanceXm_ZABpua,alcanceXm_ZBCpua,alcanceXm_ZCApua,alcanceXm_ZAGpua,alcanceXm_ZBGpua,alcanceXm_ZCGpua]=...
subalcance(Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG,zl1u,zLma);
OrdenadaXm_ZABpua(1,i)=alcanceXm_ZABpua;
OrdenadaXm_ZBCpua(1,i)=alcanceXm_ZBCpua;
OrdenadaXm_ZCApua(1,i)=alcanceXm_ZCApua;
OrdenadaXm_ZAGpua(1,i)=alcanceXm_ZAGpua;
OrdenadaXm_ZBGpua(1,i)=alcanceXm_ZBGpua;
OrdenadaXm_ZCGpua(1,i)=alcanceXm_ZCGpua;

[alcanceRp_ZABpua,alcanceRp_ZBCpua,alcanceRp_ZCApua,alcanceRp_ZAGpua,alcanceRp_ZBGpua,alcanceRp_ZCGpua]=...
subalcance(Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG,zl1u,zLpa);
OrdenadaRp_ZABpua(1,i)=alcanceRp_ZABpua;
OrdenadaRp_ZBCpua(1,i)=alcanceRp_ZBCpua;
OrdenadaRp_ZCApua(1,i)=alcanceRp_ZCApua;
OrdenadaRp_ZAGpua(1,i)=alcanceRp_ZAGpua;
OrdenadaRp_ZBGpua(1,i)=alcanceRp_ZBGpua;
OrdenadaRp_ZCGpua(1,i)=alcanceRp_ZCGpua;

[alcanceXp_ZABpua,alcanceXp_ZBCpua,alcanceXp_ZCApua,alcanceXp_ZAGpua,alcanceXp_ZBGpua,alcanceXp_ZCGpua]=...
subalcance(Xp_ZAB,Xp_ZBC,Xp_ZCA,Xp_ZAG,Xp_ZBG,Xp_ZCG,zl1u,zLpa);
OrdenadaXp_ZABpua(1,i)=alcanceXp_ZABpua;
OrdenadaXp_ZBCpua(1,i)=alcanceXp_ZBCpua;
OrdenadaXp_ZCApua(1,i)=alcanceXp_ZCApua;
OrdenadaXp_ZAGpua(1,i)=alcanceXp_ZAGpua;
OrdenadaXp_ZBGpua(1,i)=alcanceXp_ZBGpua;
OrdenadaXp_ZCGpua(1,i)=alcanceXp_ZCGpua;

%***************************************
%Para salida 67N
%***************************************
Vpol_m_izq=Vspol;
Iop_m_izq=Ismop;
AngDir_m_izq=AngDirSm;

Vpol_m_der=Vxpol;
Iop_m_der=Ixmop;
AngDir_m_der=AngDirXm;

cadVpol_m_izq='modulo Vspol(kV)';
cadIop_m_izq='modulo Ismop(kA)';
cadAngDir_m_izq='AngDirSm(º)';

cadVpol_m_der='modulo Vxpol(kV)';
cadIop_m_der='modulo Ixmop(kV)';
cadAngDir_m_der='AngDirXm(º)';

Vpol_p_izq=Vxpol;
Iop_p_izq=Ixpop;
AngDir_p_izq=AngDirXp;

Vpol_p_der=Vrpol;
Iop_p_der=Irpop;
AngDir_p_der=AngDirRp;

cadVpol_p_izq='modulo Vxpol(kV)';
cadIop_p_izq='modulo Ixpop(kA)';
cadAngDir_p_izq='AngDirXp(º)';

cadVpol_p_der='modulo Vrpol(kV)';
cadIop_p_der='modulo Irpop(kV)';
cadAngDir_p_der='AngDirRp(º)';

%***************************************
% Para Salida 4  R-X (Tramo 4)
%***************************************
ZAB_m1=Sm_ZAB;
ZBC_m1=Sm_ZBC;
ZCA_m1=Sm_ZCA;
ZAG_m1=Sm_ZAG;
ZBG_m1=Sm_ZBG;
ZCG_m1=Sm_ZCG;

ZAB_m2=Xm_ZAB;
ZBC_m2=Xm_ZBC;
ZCA_m2=Xm_ZCA;
ZAG_m2=Xm_ZAG;
ZBG_m2=Xm_ZBG;
ZCG_m2=Xm_ZCG;

ZAB_p1=Rp_ZAB;
ZBC_p1=Rp_ZBC;
ZCA_p1=Rp_ZCA;
ZAG_p1=Rp_ZAG;
ZBG_p1=Rp_ZBG;
ZCG_p1=Rp_ZCG;

ZAB_p2=Xp_ZAB;
ZBC_p2=Xp_ZBC;
ZCA_p2=Xp_ZCA;
ZAG_p2=Xp_ZAG;
ZBG_p2=Xp_ZBG;
ZCG_p2=Xp_ZCG;

%***************************************

case{5}, %tramo5
   
%case{5}, % SY,RY
%Lma=L1+L2+L5;
%Lpa=L3+L5;
   
Ism0=Ie0(5);Iym0=Ie0(7);Irp0=Ie0(9);Iyp0=Ie0(10);
Ism1=Ie1(5);Iym1=Ie1(7);Irp1=Ie1(9);Iyp1=Ie1(10);
Ism2=Ie2(5);Iym2=Ie2(7);Irp2=Ie2(9);Iyp2=Ie2(10);
Isx0=Ie0(6);
Irx0=-Ie0(8);
   
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2);
[Iyma,Iymb,Iymc]=abc(Iym0,Iym1,Iym2);
[Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2);
[Iypa,Iypb,Iypc]=abc(Iyp0,Iyp1,Iyp2);

switch(compens21)
case{1},
% COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
  PDistComp(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k,Km,Isx0);
[Ym_ZAB,Ym_ZBC,Ym_ZCA,Ym_ZAG,Ym_ZBG,Ym_ZCG]=...
   PDistComp(Vya,Vyb,Vyc,Iyma,Iymb,Iymc,Iym0,k,Km,Iyp0);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
   PDistComp(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k,Km,Irx0);
[Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG]=...
   PDistComp(Vya,Vyb,Vyc,Iypa,Iypb,Iypc,Iyp0,k,Km,Iym0);
   
case{2},
% NO COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
  PDist(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k);
[Ym_ZAB,Ym_ZBC,Ym_ZCA,Ym_ZAG,Ym_ZBG,Ym_ZCG]=...
   PDist(Vya,Vyb,Vyc,Iyma,Iymb,Iymc,Iym0,k);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
   PDist(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k);
[Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG]=...
   PDist(Vya,Vyb,Vyc,Iypa,Iypb,Iypc,Iyp0,k);
end

% MEDICION 67N EN LOS EXTREMOS DE LINEA
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0);
[Vypol,Iymop,AngDirYm]=PDire(Vy0,Iym0);

[Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0);

[Vypol,Iypop,AngDirYp]=PDire(Vy0,Iyp0);

% Recuerdo: tramo5

zLSm=abs(zl1u*(L1+L2+m*L5));
zLYm=abs(zl1u*(1-m)*L5);
zLRp=abs(zl1u*(L3+p*L5));
zLYp=abs(zl1u*(1-p)*L5);

[alcanceSm_ZABpu,alcanceSm_ZBCpu,alcanceSm_ZCApu,alcanceSm_ZAGpu,alcanceSm_ZBGpu,alcanceSm_ZCGpu]=...
subalcance(Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG,zl1u,zLSm);
OrdenadaSm_ZABpu(1,i)=alcanceSm_ZABpu;
OrdenadaSm_ZBCpu(1,i)=alcanceSm_ZBCpu;
OrdenadaSm_ZCApu(1,i)=alcanceSm_ZCApu;
OrdenadaSm_ZAGpu(1,i)=alcanceSm_ZAGpu;
OrdenadaSm_ZBGpu(1,i)=alcanceSm_ZBGpu;

OrdenadaSm_ZCGpu(1,i)=alcanceSm_ZCGpu;

[alcanceYm_ZABpu,alcanceYm_ZBCpu,alcanceYm_ZCApu,alcanceYm_ZAGpu,alcanceYm_ZBGpu,alcanceYm_ZCGpu]=...
   subalcance(Ym_ZAB,Ym_ZBC,Ym_ZCA,Ym_ZAG,Ym_ZBG,Ym_ZCG,zl1u,zLYm);
OrdenadaYm_ZABpu(1,i)=alcanceYm_ZABpu;
OrdenadaYm_ZBCpu(1,i)=alcanceYm_ZBCpu;
OrdenadaYm_ZCApu(1,i)=alcanceYm_ZCApu;
OrdenadaYm_ZAGpu(1,i)=alcanceYm_ZAGpu;
OrdenadaYm_ZBGpu(1,i)=alcanceYm_ZBGpu;
OrdenadaYm_ZCGpu(1,i)=alcanceYm_ZCGpu;

[alcanceRp_ZABpu,alcanceRp_ZBCpu,alcanceRp_ZCApu,alcanceRp_ZAGpu,alcanceRp_ZBGpu,alcanceRp_ZCGpu]=...
   subalcance(Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG,zl1u,zLRp);
OrdenadaRp_ZABpu(1,i)=alcanceRp_ZABpu;
OrdenadaRp_ZBCpu(1,i)=alcanceRp_ZBCpu;
OrdenadaRp_ZCApu(1,i)=alcanceRp_ZCApu;
OrdenadaRp_ZAGpu(1,i)=alcanceRp_ZAGpu;
OrdenadaRp_ZBGpu(1,i)=alcanceRp_ZBGpu;
OrdenadaRp_ZCGpu(1,i)=alcanceRp_ZCGpu;

[alcanceYp_ZABpu,alcanceYp_ZBCpu,alcanceYp_ZCApu,alcanceYp_ZAGpu,alcanceYp_ZBGpu,alcanceYp_ZCGpu]=...
   subalcance(Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG,zl1u,zLYp);
OrdenadaYp_ZABpu(1,i)=alcanceYp_ZABpu;
OrdenadaYp_ZBCpu(1,i)=alcanceYp_ZBCpu;
OrdenadaYp_ZCApu(1,i)=alcanceYp_ZCApu;
OrdenadaYp_ZAGpu(1,i)=alcanceYp_ZAGpu;
OrdenadaYp_ZBGpu(1,i)=alcanceYp_ZBGpu;
OrdenadaYp_ZCGpu(1,i)=alcanceYp_ZCGpu;

Abscisa(fila,i)=m;

[alcanceSm_ZABpua,alcanceSm_ZBCpua,alcanceSm_ZCApua,alcanceSm_ZAGpua,alcanceSm_ZBGpua,alcanceSm_ZCGpua]=...
   subalcance(Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG,zl1u,zLma);
OrdenadaSm_ZABpua(1,i)=alcanceSm_ZABpua;
OrdenadaSm_ZBCpua(1,i)=alcanceSm_ZBCpua;
OrdenadaSm_ZCApua(1,i)=alcanceSm_ZCApua;
OrdenadaSm_ZAGpua(1,i)=alcanceSm_ZAGpua;
OrdenadaSm_ZBGpua(1,i)=alcanceSm_ZBGpua;

OrdenadaSm_ZCGpua(1,i)=alcanceSm_ZCGpua;

[alcanceYm_ZABpua,alcanceYm_ZBCpua,alcanceYm_ZCApua,alcanceYm_ZAGpua,alcanceYm_ZBGpua,alcanceYm_ZCGpua]=...
   subalcance(Ym_ZAB,Ym_ZBC,Ym_ZCA,Ym_ZAG,Ym_ZBG,Ym_ZCG,zl1u,zLma);
OrdenadaYm_ZABpua(1,i)=alcanceYm_ZABpua;
OrdenadaYm_ZBCpua(1,i)=alcanceYm_ZBCpua;
OrdenadaYm_ZCApua(1,i)=alcanceYm_ZCApua;
OrdenadaYm_ZAGpua(1,i)=alcanceYm_ZAGpua;
OrdenadaYm_ZBGpua(1,i)=alcanceYm_ZBGpua;
OrdenadaYm_ZCGpua(1,i)=alcanceYm_ZCGpua;

[alcanceRp_ZABpua,alcanceRp_ZBCpua,alcanceRp_ZCApua,alcanceRp_ZAGpua,alcanceRp_ZBGpua,alcanceRp_ZCGpua]=...
   subalcance(Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG,zl1u,zLpa);
OrdenadaRp_ZABpua(1,i)=alcanceRp_ZABpua;
OrdenadaRp_ZBCpua(1,i)=alcanceRp_ZBCpua;
OrdenadaRp_ZCApua(1,i)=alcanceRp_ZCApua;
OrdenadaRp_ZAGpua(1,i)=alcanceRp_ZAGpua;
OrdenadaRp_ZBGpua(1,i)=alcanceRp_ZBGpua;
OrdenadaRp_ZCGpua(1,i)=alcanceRp_ZCGpua;

[alcanceYp_ZABpua,alcanceYp_ZBCpua,alcanceYp_ZCApua,alcanceYp_ZAGpua,alcanceYp_ZBGpua,alcanceYp_ZCGpua]=...
   subalcance(Yp_ZAB,Yp_ZBC,Yp_ZCA,Yp_ZAG,Yp_ZBG,Yp_ZCG,zl1u,zLpa);
OrdenadaYp_ZABpua(1,i)=alcanceYp_ZABpua;
OrdenadaYp_ZBCpua(1,i)=alcanceYp_ZBCpua;
OrdenadaYp_ZCApua(1,i)=alcanceYp_ZCApua;
OrdenadaYp_ZAGpua(1,i)=alcanceYp_ZAGpua;
OrdenadaYp_ZBGpua(1,i)=alcanceYp_ZBGpua;
OrdenadaYp_ZCGpua(1,i)=alcanceYp_ZCGpua;

%***************************************
%Para salida 67N
%***************************************
Vpol_m_izq=Vspol;
Iop_m_izq=Ismop;
AngDir_m_izq=AngDirSm;

Vpol_m_der=Vypol;
Iop_m_der=Iymop;
AngDir_m_der=AngDirYm;

cadVpol_m_izq='modulo Vspol(kV)';
cadIop_m_izq='modulo Ismop(kA)';
cadAngDir_m_izq='AngDirSm(º)';

cadVpol_m_der='modulo Vypol(kV)';
cadIop_m_der='modulo Iymop(kV)';
cadAngDir_m_der='AngDirYm(º)';

Vpol_p_izq=Vypol;
Iop_p_izq=Iypop;
AngDir_p_izq=AngDirYp;

Vpol_p_der=Vrpol;
Iop_p_der=Irpop;
AngDir_p_der=AngDirRp;

cadVpol_p_izq='modulo Vypol(kV)';
cadIop_p_izq='modulo Iypop(kA)';
cadAngDir_p_izq='AngDirYp(º)';

cadVpol_p_der='modulo Vrpol(kV)';
cadIop_p_der='modulo Irpop(kV)';
cadAngDir_p_der='AngDirRp(º)';

%***************************************
% Para Salida 4  R-X (Tramo 5)
%***************************************
ZAB_m1=Sm_ZAB;
ZBC_m1=Sm_ZBC;
ZCA_m1=Sm_ZCA;
ZAG_m1=Sm_ZAG;
ZBG_m1=Sm_ZBG;
ZCG_m1=Sm_ZCG;

ZAB_m2=Ym_ZAB;
ZBC_m2=Ym_ZBC;
ZCA_m2=Ym_ZCA;
ZAG_m2=Ym_ZAG;
ZBG_m2=Ym_ZBG;
ZCG_m2=Ym_ZCG;

ZAB_p1=Rp_ZAB;
ZBC_p1=Rp_ZBC;
ZCA_p1=Rp_ZCA;
ZAG_p1=Rp_ZAG;
ZBG_p1=Rp_ZBG;
ZCG_p1=Rp_ZCG;

ZAB_p2=Yp_ZAB;
ZBC_p2=Yp_ZBC;
ZCA_p2=Yp_ZCA;
ZAG_p2=Yp_ZAG;
ZBG_p2=Yp_ZBG;
ZCG_p2=Yp_ZCG;

%***************************************

end % end del switch(tramo)

%***************************************
%Para salida 3, 67N
%***************************************
vectorVpol_m_izq(1,i)=Vpol_m_izq;
vectorIop_m_izq(1,i)=Iop_m_izq;
vectorAngDir_m_izq(1,i)=AngDir_m_izq;

vectorVpol_m_der(1,i)=Vpol_m_der;
vectorIop_m_der(1,i)=Iop_m_der;
vectorAngDir_m_der(1,i)=AngDir_m_der;

vectorVpol_p_izq(1,i)=Vpol_p_izq;
vectorIop_p_izq(1,i)=Iop_p_izq;
vectorAngDir_p_izq(1,i)=AngDir_p_izq;

vectorVpol_p_der(1,i)=Vpol_p_der;
vectorIop_p_der(1,i)=Iop_p_der;
vectorAngDir_p_der(1,i)=AngDir_p_der;

%***************************************
% Para Salida 4  R-X
%***************************************
vectorZAB_m1(1,i)=ZAB_m1;
vectorZBC_m1(1,i)=ZBC_m1;
vectorZCA_m1(1,i)=ZCA_m1;
vectorZAG_m1(1,i)=ZAG_m1;
vectorZBG_m1(1,i)=ZBG_m1;
vectorZCG_m1(1,i)=ZCG_m1;

vectorZAB_m2(1,i)=ZAB_m2;
vectorZBC_m2(1,i)=ZBC_m2;
vectorZCA_m2(1,i)=ZCA_m2;
vectorZAG_m2(1,i)=ZAG_m2;
vectorZBG_m2(1,i)=ZBG_m2;
vectorZCG_m2(1,i)=ZCG_m2;

vectorZAB_p1(1,i)=ZAB_p1;
vectorZBC_p1(1,i)=ZBC_p1;
vectorZCA_p1(1,i)=ZCA_p1;
vectorZAG_p1(1,i)=ZAG_p1;
vectorZBG_p1(1,i)=ZBG_p1;
vectorZCG_p1(1,i)=ZCG_p1;

vectorZAB_p2(1,i)=ZAB_p2;
vectorZBC_p2(1,i)=ZBC_p2;
vectorZCA_p2(1,i)=ZCA_p2;
vectorZAG_p2(1,i)=ZAG_p2;
vectorZBG_p2(1,i)=ZBG_p2;
vectorZCG_p2(1,i)=ZCG_p2;

%***************************************

i=i+1;
m=m+0.01;
end  % end del while m


% ------------------------------------------------------
% SALIDA
% ------------------------------------------------------

variassalidas=1;
while variassalidas==1

switch(tramo)
   
case{1}, %  ******** salidas para falta en el tramo 1 ******** 
   
%case{1}, % SX,SY
%Lma=L1+L4;
%Lpa=L1+L2+L5;
fprintf('\n');
fprintf('Salida Grafica 1: Extremos m: Todas las Unidades (pu)\n');
fprintf('                : Extremos p: Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2: Extremos m: Todas las Unidades (pua)\n');
fprintf('                : Extremos p: Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3: Analisis Extremo m Izq, 67N\n');
fprintf('                : Analisis Extremo m Der, 67N\n');
fprintf('                : Analisis Extremo p Izq, 67N\n');
fprintf('                : Analisis Extremo p Der, 67N\n');
fprintf('Salida Grafica 4: Analisis Extremo m Izq, R-X\n');
fprintf('                : Analisis Extremo m Der, R-X\n');
fprintf('                : Analisis Extremo p Izq, R-X\n');
fprintf('                : Analisis Extremo p Der, R-X\n');
fprintf('\n');
   
salida=input('Salida: 1,2,3,4 :');

switch(salida)
case{1}, % salida 1, para el tramo 1
   
cadenaym1=sprintf('alcanceSmpu');
cadenaym2=sprintf('alcanceXmpu');
cadenayp1=sprintf('alcanceSppu');
cadenayp2=sprintf('alcanceYppu');

Ordenada1mAB=OrdenadaSm_ZABpu(1:1,1:(i-1));
Ordenada1mBC=OrdenadaSm_ZBCpu(1:1,1:(i-1));
Ordenada1mCA=OrdenadaSm_ZCApu(1:1,1:(i-1));
Ordenada1mAG=OrdenadaSm_ZAGpu(1:1,1:(i-1));
Ordenada1mBG=OrdenadaSm_ZBGpu(1:1,1:(i-1));
Ordenada1mCG=OrdenadaSm_ZCGpu(1:1,1:(i-1));

Ordenada4mAB=OrdenadaXm_ZABpu(1:1,1:(i-1));
Ordenada4mBC=OrdenadaXm_ZBCpu(1:1,1:(i-1));
Ordenada4mCA=OrdenadaXm_ZCApu(1:1,1:(i-1));
Ordenada4mAG=OrdenadaXm_ZAGpu(1:1,1:(i-1));
Ordenada4mBG=OrdenadaXm_ZBGpu(1:1,1:(i-1));
Ordenada4mCG=OrdenadaXm_ZCGpu(1:1,1:(i-1));

Ordenada7mAB=OrdenadaSp_ZABpu(1:1,1:(i-1));
Ordenada7mBC=OrdenadaSp_ZBCpu(1:1,1:(i-1));
Ordenada7mCA=OrdenadaSp_ZCApu(1:1,1:(i-1));
Ordenada7mAG=OrdenadaSp_ZAGpu(1:1,1:(i-1));
Ordenada7mBG=OrdenadaSp_ZBGpu(1:1,1:(i-1));
Ordenada7mCG=OrdenadaSp_ZCGpu(1:1,1:(i-1));

Ordenada10mAB=OrdenadaYp_ZABpu(1:1,1:(i-1));
Ordenada10mBC=OrdenadaYp_ZBCpu(1:1,1:(i-1));
Ordenada10mCA=OrdenadaYp_ZCApu(1:1,1:(i-1));
Ordenada10mAG=OrdenadaYp_ZAGpu(1:1,1:(i-1));
Ordenada10mBG=OrdenadaYp_ZBGpu(1:1,1:(i-1));
Ordenada10mCG=OrdenadaYp_ZCGpu(1:1,1:(i-1));


       case{2}, % salida 2, para el tramo 1
          
cadenaym1a=sprintf('alcanceSmpua');
cadenaym2a=sprintf('alcanceXmpua');
cadenayp1a=sprintf('alcanceSppua');
cadenayp2a=sprintf('alcanceYppua');
          
Ordenada1mABa=OrdenadaSm_ZABpua(1:1,1:(i-1));
Ordenada1mBCa=OrdenadaSm_ZBCpua(1:1,1:(i-1));
Ordenada1mCAa=OrdenadaSm_ZCApua(1:1,1:(i-1));
Ordenada1mAGa=OrdenadaSm_ZAGpua(1:1,1:(i-1));
Ordenada1mBGa=OrdenadaSm_ZBGpua(1:1,1:(i-1));
Ordenada1mCGa=OrdenadaSm_ZCGpua(1:1,1:(i-1));

Ordenada4mABa=OrdenadaXm_ZABpua(1:1,1:(i-1));
Ordenada4mBCa=OrdenadaXm_ZBCpua(1:1,1:(i-1));
Ordenada4mCAa=OrdenadaXm_ZCApua(1:1,1:(i-1));
Ordenada4mAGa=OrdenadaXm_ZAGpua(1:1,1:(i-1));
Ordenada4mBGa=OrdenadaXm_ZBGpua(1:1,1:(i-1));
Ordenada4mCGa=OrdenadaXm_ZCGpua(1:1,1:(i-1));

Ordenada7mABa=OrdenadaSp_ZABpua(1:1,1:(i-1));
Ordenada7mBCa=OrdenadaSp_ZBCpua(1:1,1:(i-1));
Ordenada7mCAa=OrdenadaSp_ZCApua(1:1,1:(i-1));
Ordenada7mAGa=OrdenadaSp_ZAGpua(1:1,1:(i-1));
Ordenada7mBGa=OrdenadaSp_ZBGpua(1:1,1:(i-1));
Ordenada7mCGa=OrdenadaSp_ZCGpua(1:1,1:(i-1));

Ordenada10mABa=OrdenadaYp_ZABpua(1:1,1:(i-1));
Ordenada10mBCa=OrdenadaYp_ZBCpua(1:1,1:(i-1));
Ordenada10mCAa=OrdenadaYp_ZCApua(1:1,1:(i-1));
Ordenada10mAGa=OrdenadaYp_ZAGpua(1:1,1:(i-1));
Ordenada10mBGa=OrdenadaYp_ZBGpua(1:1,1:(i-1));
Ordenada10mCGa=OrdenadaYp_ZCGpua(1:1,1:(i-1));

% Tramo 1
% SX:m (S=Extremo primero X=Extremo segundo) 
% SY:m (S=Extremo primero Y=Extremo segundo)

% amp: m=linea con falta m p=extremo primero  
% ams:                     s=extremo segundo  
% app: p=linea con falta p p=extremo primero  
% aps:                     s=extremo segundo 

amp=0; bmp=L1/Lma; ams=(L4+L1)/Lma; bms=-L1/Lma;
app=0; bpp=L1/Lpa; aps=(L5+L2+L1)/Lpa; bps=-L1/Lpa;

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam0(cont)=mm;
%Ordenadam0(cont)=mm*(L1/(L1+L4));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam00(cont)=mm;
%Ordenadam00(cont)=1-mm*(L1/(L1+L4));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap0(cont)=mm;
%Ordenadap0(cont)=mm*(L1/(L1+L2+L5));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap00(cont)=mm;
%Ordenadap00(cont)=1-mm*(L1/(L1+L2+L5));
%cont=cont+1;
%mm=mm+0.01;
%end

end % end del switch(salida), para falta en el tramo 1 

case{2}, %  ******** salidas para falta en el tramo 2 ******** 
   
%case{2}, % XR,SY
%Lma=L4+L2+L3;
%Lpa=L1+L2+L5;

fprintf('\n');
fprintf('Salida Grafica 1: Extremos m: Todas las Unidades (pu)\n');
fprintf('                : Extremos p: Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2: Extremos m: Todas las Unidades (pua)\n');
fprintf('                : Extremos p: Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3: Analisis Extremo m Izq, 67N\n');
fprintf('                : Analisis Extremo m Der, 67N\n');
fprintf('                : Analisis Extremo p Izq, 67N\n');
fprintf('                : Analisis Extremo p Der, 67N\n');
fprintf('Salida Grafica 4: Analisis Extremo m Izq, R-X\n');
fprintf('                : Analisis Extremo m Der, R-X\n');
fprintf('                : Analisis Extremo p Izq, R-X\n');
fprintf('                : Analisis Extremo p Der, R-X\n');
fprintf('\n');
   

salida=input('Salida: 1,2,3,4 :');

switch(salida)
case{1}, % salida 1, para el tramo 2
   
cadenaym1=sprintf('alcanceXmpu');
cadenaym2=sprintf('alcanceRmpu');
cadenayp1=sprintf('alcanceSppu');
cadenayp2=sprintf('alcanceYppu');
   
Ordenada1mAB=OrdenadaXm_ZABpu(1:1,1:(i-1));
Ordenada1mBC=OrdenadaXm_ZBCpu(1:1,1:(i-1));
Ordenada1mCA=OrdenadaXm_ZCApu(1:1,1:(i-1));
Ordenada1mAG=OrdenadaXm_ZAGpu(1:1,1:(i-1));
Ordenada1mBG=OrdenadaXm_ZBGpu(1:1,1:(i-1));
Ordenada1mCG=OrdenadaXm_ZCGpu(1:1,1:(i-1));
   
Ordenada4mAB=OrdenadaRm_ZABpu(1:1,1:(i-1));
Ordenada4mBC=OrdenadaRm_ZBCpu(1:1,1:(i-1));
Ordenada4mCA=OrdenadaRm_ZCApu(1:1,1:(i-1));
Ordenada4mAG=OrdenadaRm_ZAGpu(1:1,1:(i-1));
Ordenada4mBG=OrdenadaRm_ZBGpu(1:1,1:(i-1));
Ordenada4mCG=OrdenadaRm_ZCGpu(1:1,1:(i-1));

Ordenada7mAB=OrdenadaSp_ZABpu(1:1,1:(i-1));
Ordenada7mBC=OrdenadaSp_ZBCpu(1:1,1:(i-1));
Ordenada7mCA=OrdenadaSp_ZCApu(1:1,1:(i-1));
Ordenada7mAG=OrdenadaSp_ZAGpu(1:1,1:(i-1));
Ordenada7mBG=OrdenadaSp_ZBGpu(1:1,1:(i-1));
Ordenada7mCG=OrdenadaSp_ZCGpu(1:1,1:(i-1));

Ordenada10mAB=OrdenadaYp_ZABpu(1:1,1:(i-1));
Ordenada10mBC=OrdenadaYp_ZBCpu(1:1,1:(i-1));
Ordenada10mCA=OrdenadaYp_ZCApu(1:1,1:(i-1));
Ordenada10mAG=OrdenadaYp_ZAGpu(1:1,1:(i-1));
Ordenada10mBG=OrdenadaYp_ZBGpu(1:1,1:(i-1));
Ordenada10mCG=OrdenadaYp_ZCGpu(1:1,1:(i-1));

case{2}, % salida 2, para el tramo 2
   
cadenaym1a=sprintf('alcanceXmpua');
cadenaym2a=sprintf('alcanceRmpua');
cadenayp1a=sprintf('alcanceSppua');
cadenayp2a=sprintf('alcanceYppua');
   
Ordenada1mABa=OrdenadaXm_ZABpua(1:1,1:(i-1));
Ordenada1mBCa=OrdenadaXm_ZBCpua(1:1,1:(i-1));
Ordenada1mCAa=OrdenadaXm_ZCApua(1:1,1:(i-1));
Ordenada1mAGa=OrdenadaXm_ZAGpua(1:1,1:(i-1));
Ordenada1mBGa=OrdenadaXm_ZBGpua(1:1,1:(i-1));
Ordenada1mCGa=OrdenadaXm_ZCGpua(1:1,1:(i-1));
          
Ordenada4mABa=OrdenadaRm_ZABpua(1:1,1:(i-1));
Ordenada4mBCa=OrdenadaRm_ZBCpua(1:1,1:(i-1));
Ordenada4mCAa=OrdenadaRm_ZCApua(1:1,1:(i-1));
Ordenada4mAGa=OrdenadaRm_ZAGpua(1:1,1:(i-1));
Ordenada4mBGa=OrdenadaRm_ZBGpua(1:1,1:(i-1));
Ordenada4mCGa=OrdenadaRm_ZCGpua(1:1,1:(i-1));

Ordenada7mABa=OrdenadaSp_ZABpua(1:1,1:(i-1));
Ordenada7mBCa=OrdenadaSp_ZBCpua(1:1,1:(i-1));
Ordenada7mCAa=OrdenadaSp_ZCApua(1:1,1:(i-1));
Ordenada7mAGa=OrdenadaSp_ZAGpua(1:1,1:(i-1));
Ordenada7mBGa=OrdenadaSp_ZBGpua(1:1,1:(i-1));
Ordenada7mCGa=OrdenadaSp_ZCGpua(1:1,1:(i-1));

Ordenada10mABa=OrdenadaYp_ZABpua(1:1,1:(i-1));
Ordenada10mBCa=OrdenadaYp_ZBCpua(1:1,1:(i-1));
Ordenada10mCAa=OrdenadaYp_ZCApua(1:1,1:(i-1));
Ordenada10mAGa=OrdenadaYp_ZAGpua(1:1,1:(i-1));
Ordenada10mBGa=OrdenadaYp_ZBGpua(1:1,1:(i-1));
Ordenada10mCGa=OrdenadaYp_ZCGpua(1:1,1:(i-1));

% Tramo 2
% XR:m (X=Extremo primero R=Extremo segundo) 
% SY:p (S=Extremo primero Y=Extremo segundo)

% amp: m=linea con falta m p=extremo primero  
% ams:                     s=extremo segundo  
% app: p=linea con falta p p=extremo primero  
% aps:                     s=extremo segundo 

amp=L4/Lma; bmp=L2/Lma; ams=(L3+L2)/Lma; bms=-L2/Lma;
app=L1/Lpa; bpp=L2/Lpa; aps=(L5+L2)/Lpa; bps=-L2/Lpa;

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam0(cont)=mm;
%Ordenadam0(cont)=L4/Lma+mm*(L2/Lma);
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam00(cont)=mm;
%Ordenadam00(cont)=(L3+L2)/Lma-mm*(L2/Lma);
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap0(cont)=mm;
%Ordenadap0(cont)=L1/Lpa+mm*(L2/Lpa);
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap00(cont)=mm;
%Ordenadap00(cont)=(L5+L2)/Lpa-mm*(L2/Lpa);
%cont=cont+1;
%mm=mm+0.01;
%end

end % end del switch(salida), para el tramo 2 

case{3}, %  ******** salidas para falta en el tramo 3 ********
   
%case{3}, % XR,YR
%Lma=L4+L2+L3;
%Lpa=L3+L5;
   

fprintf('\n');
fprintf('Salida Grafica 1: Extremos m: Todas las Unidades (pu)\n');
fprintf('                : Extremos p: Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2: Extremos m: Todas las Unidades (pua)\n');
fprintf('                : Extremos p: Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3: Analisis Extremo m Izq, 67N\n');
fprintf('                : Analisis Extremo m Der, 67N\n');
fprintf('                : Analisis Extremo p Izq, 67N\n');
fprintf('                : Analisis Extremo p Der, 67N\n');
fprintf('Salida Grafica 4: Analisis Extremo m Izq, R-X\n');
fprintf('                : Analisis Extremo m Der, R-X\n');
fprintf('                : Analisis Extremo p Izq, R-X\n');
fprintf('                : Analisis Extremo p Der, R-X\n');
fprintf('\n');

salida=input('Salida: 1,2,3,4 :');

switch(salida)
         
case{1}, % salida 1, para el tramo 3
   
cadenaym1=sprintf('alcanceXmpu');
cadenaym2=sprintf('alcanceRmpu');
cadenayp1=sprintf('alcanceYppu');
cadenayp2=sprintf('alcanceRppu');

Ordenada1mAB=OrdenadaXm_ZABpu(1:1,1:(i-1));
Ordenada1mBC=OrdenadaXm_ZBCpu(1:1,1:(i-1));
Ordenada1mCA=OrdenadaXm_ZCApu(1:1,1:(i-1));
Ordenada1mAG=OrdenadaXm_ZAGpu(1:1,1:(i-1));
Ordenada1mBG=OrdenadaXm_ZBGpu(1:1,1:(i-1));
Ordenada1mCG=OrdenadaXm_ZCGpu(1:1,1:(i-1));
   
Ordenada4mAB=OrdenadaRm_ZABpu(1:1,1:(i-1));
Ordenada4mBC=OrdenadaRm_ZBCpu(1:1,1:(i-1));
Ordenada4mCA=OrdenadaRm_ZCApu(1:1,1:(i-1));
Ordenada4mAG=OrdenadaRm_ZAGpu(1:1,1:(i-1));
Ordenada4mBG=OrdenadaRm_ZBGpu(1:1,1:(i-1));
Ordenada4mCG=OrdenadaRm_ZCGpu(1:1,1:(i-1));

Ordenada7mAB=OrdenadaYp_ZABpu(1:1,1:(i-1));
Ordenada7mBC=OrdenadaYp_ZBCpu(1:1,1:(i-1));
Ordenada7mCA=OrdenadaYp_ZCApu(1:1,1:(i-1));
Ordenada7mAG=OrdenadaYp_ZAGpu(1:1,1:(i-1));
Ordenada7mBG=OrdenadaYp_ZBGpu(1:1,1:(i-1));
Ordenada7mCG=OrdenadaYp_ZCGpu(1:1,1:(i-1));

Ordenada10mAB=OrdenadaRp_ZABpu(1:1,1:(i-1));
Ordenada10mBC=OrdenadaRp_ZBCpu(1:1,1:(i-1));
Ordenada10mCA=OrdenadaRp_ZCApu(1:1,1:(i-1));
Ordenada10mAG=OrdenadaRp_ZAGpu(1:1,1:(i-1));
Ordenada10mBG=OrdenadaRp_ZBGpu(1:1,1:(i-1));
Ordenada10mCG=OrdenadaRp_ZCGpu(1:1,1:(i-1));


case{2}, % salida 2, para el tramo 3
   
   
cadenaym1a=sprintf('alcanceXmpua');
cadenaym2a=sprintf('alcanceRmpua');
cadenayp1a=sprintf('alcanceYppua');
cadenayp2a=sprintf('alcanceRppua');
   
Ordenada1mABa=OrdenadaXm_ZABpua(1:1,1:(i-1));
Ordenada1mBCa=OrdenadaXm_ZBCpua(1:1,1:(i-1));
Ordenada1mCAa=OrdenadaXm_ZCApua(1:1,1:(i-1));
Ordenada1mAGa=OrdenadaXm_ZAGpua(1:1,1:(i-1));
Ordenada1mBGa=OrdenadaXm_ZBGpua(1:1,1:(i-1));
Ordenada1mCGa=OrdenadaXm_ZCGpua(1:1,1:(i-1));
          
Ordenada4mABa=OrdenadaRm_ZABpua(1:1,1:(i-1));
Ordenada4mBCa=OrdenadaRm_ZBCpua(1:1,1:(i-1));
Ordenada4mCAa=OrdenadaRm_ZCApua(1:1,1:(i-1));
Ordenada4mAGa=OrdenadaRm_ZAGpua(1:1,1:(i-1));
Ordenada4mBGa=OrdenadaRm_ZBGpua(1:1,1:(i-1));
Ordenada4mCGa=OrdenadaRm_ZCGpua(1:1,1:(i-1));

Ordenada7mABa=OrdenadaYp_ZABpua(1:1,1:(i-1));
Ordenada7mBCa=OrdenadaYp_ZBCpua(1:1,1:(i-1));
Ordenada7mCAa=OrdenadaYp_ZCApua(1:1,1:(i-1));
Ordenada7mAGa=OrdenadaYp_ZAGpua(1:1,1:(i-1));
Ordenada7mBGa=OrdenadaYp_ZBGpua(1:1,1:(i-1));
Ordenada7mCGa=OrdenadaYp_ZCGpua(1:1,1:(i-1));

Ordenada10mABa=OrdenadaRp_ZABpua(1:1,1:(i-1));
Ordenada10mBCa=OrdenadaRp_ZBCpua(1:1,1:(i-1));
Ordenada10mCAa=OrdenadaRp_ZCApua(1:1,1:(i-1));
Ordenada10mAGa=OrdenadaRp_ZAGpua(1:1,1:(i-1));
Ordenada10mBGa=OrdenadaRp_ZBGpua(1:1,1:(i-1));
Ordenada10mCGa=OrdenadaRp_ZCGpua(1:1,1:(i-1));

% Tramo 3
% XR:m (X=Extremo primero R=Extremo segundo) 
% YR:p (Y=Extremo primero R=Extremo segundo)

% amp: m=linea con falta m p=extremo primero  
% ams:                     s=extremo segundo  
% app: p=linea con falta p p=extremo primero  
% aps:                     s=extremo segundo 

amp=(L4+L2)/Lma; bmp=L3/Lma; ams=L3/Lma; bms=-L3/Lma;
app=L5/Lpa; bpp=L3/Lpa; aps=L3/Lpa; bps=-L3/Lpa;

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam0(cont)=mm;
%Ordenadam0(cont)=(L4+L2)/(L4+L2+L3)+mm*(L3/(L4+L2+L3));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam00(cont)=mm;
%Ordenadam00(cont)=L3/(L4+L3+L2)-mm*(L3/(L4+L2+L3));
%cont=cont+1;
%mm=mm+0.01;
%end


%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap0(cont)=mm;
%Ordenadap0(cont)=L5/(L5+L3)+mm*(L3/(L5+L3));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap00(cont)=mm;
%Ordenadap00(cont)=L3/(L5+L3)-mm*(L3/(L5+L3));
%cont=cont+1;
%mm=mm+0.01;
%end

end % end del switch(salida), para el tramo 3 

case{4}, %  ******** salidas para falta en el tramo 4 ********
   
%case{4}, % SX,RX
%Lma=L1+L4;
%Lpa=L4+L2+L3;

fprintf('\n');
fprintf('Salida Grafica 1: Extremos m: Todas las Unidades (pu)\n');
fprintf('                : Extremos p: Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2: Extremos m: Todas las Unidades (pua)\n');
fprintf('                : Extremos p: Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3: Analisis Extremo m Izq, 67N\n');
fprintf('                : Analisis Extremo m Der, 67N\n');
fprintf('                : Analisis Extremo p Izq, 67N\n');
fprintf('                : Analisis Extremo p Der, 67N\n');
fprintf('Salida Grafica 4: Analisis Extremo m Izq, R-X\n');
fprintf('                : Analisis Extremo m Der, R-X\n');
fprintf('                : Analisis Extremo p Izq, R-X\n');
fprintf('                : Analisis Extremo p Der, R-X\n');
fprintf('\n');

salida=input('Salida: 1,2,3,4 :');

switch(salida)
case{1}, % salida 1, para el tramo 4
   
cadenaym1=sprintf('alcanceSmpu');
cadenaym2=sprintf('alcanceXmpu');
cadenayp1=sprintf('alcanceRppu');
cadenayp2=sprintf('alcanceXppu');
          
Ordenada1mAB=OrdenadaSm_ZABpu(1:1,1:(i-1));
Ordenada1mBC=OrdenadaSm_ZBCpu(1:1,1:(i-1));
Ordenada1mCA=OrdenadaSm_ZCApu(1:1,1:(i-1));
Ordenada1mAG=OrdenadaSm_ZAGpu(1:1,1:(i-1));
Ordenada1mBG=OrdenadaSm_ZBGpu(1:1,1:(i-1));
Ordenada1mCG=OrdenadaSm_ZCGpu(1:1,1:(i-1));

Ordenada4mAB=OrdenadaXm_ZABpu(1:1,1:(i-1));
Ordenada4mBC=OrdenadaXm_ZBCpu(1:1,1:(i-1));
Ordenada4mCA=OrdenadaXm_ZCApu(1:1,1:(i-1));
Ordenada4mAG=OrdenadaXm_ZAGpu(1:1,1:(i-1));
Ordenada4mBG=OrdenadaXm_ZBGpu(1:1,1:(i-1));
Ordenada4mCG=OrdenadaXm_ZCGpu(1:1,1:(i-1));

Ordenada7mAB=OrdenadaRp_ZABpu(1:1,1:(i-1));
Ordenada7mBC=OrdenadaRp_ZBCpu(1:1,1:(i-1));
Ordenada7mCA=OrdenadaRp_ZCApu(1:1,1:(i-1));
Ordenada7mAG=OrdenadaRp_ZAGpu(1:1,1:(i-1));
Ordenada7mBG=OrdenadaRp_ZBGpu(1:1,1:(i-1));
Ordenada7mCG=OrdenadaRp_ZCGpu(1:1,1:(i-1));

Ordenada10mAB=OrdenadaXp_ZABpu(1:1,1:(i-1));
Ordenada10mBC=OrdenadaXp_ZBCpu(1:1,1:(i-1));
Ordenada10mCA=OrdenadaXp_ZCApu(1:1,1:(i-1));
Ordenada10mAG=OrdenadaXp_ZAGpu(1:1,1:(i-1));
Ordenada10mBG=OrdenadaXp_ZBGpu(1:1,1:(i-1));
Ordenada10mCG=OrdenadaXp_ZCGpu(1:1,1:(i-1));

       case{2}, % salida 2, para el tramo 4
          
cadenaym1a=sprintf('alcanceSmpua');
cadenaym2a=sprintf('alcanceXmpua');
cadenayp1a=sprintf('alcanceRppua');
cadenayp2a=sprintf('alcanceXppua');
          
Ordenada1mABa=OrdenadaSm_ZABpua(1:1,1:(i-1));
Ordenada1mBCa=OrdenadaSm_ZBCpua(1:1,1:(i-1));
Ordenada1mCAa=OrdenadaSm_ZCApua(1:1,1:(i-1));
Ordenada1mAGa=OrdenadaSm_ZAGpua(1:1,1:(i-1));
Ordenada1mBGa=OrdenadaSm_ZBGpua(1:1,1:(i-1));
Ordenada1mCGa=OrdenadaSm_ZCGpua(1:1,1:(i-1));

Ordenada4mABa=OrdenadaXm_ZABpua(1:1,1:(i-1));
Ordenada4mBCa=OrdenadaXm_ZBCpua(1:1,1:(i-1));
Ordenada4mCAa=OrdenadaXm_ZCApua(1:1,1:(i-1));
Ordenada4mAGa=OrdenadaXm_ZAGpua(1:1,1:(i-1));
Ordenada4mBGa=OrdenadaXm_ZBGpua(1:1,1:(i-1));
Ordenada4mCGa=OrdenadaXm_ZCGpua(1:1,1:(i-1));

Ordenada7mABa=OrdenadaRp_ZABpua(1:1,1:(i-1));
Ordenada7mBCa=OrdenadaRp_ZBCpua(1:1,1:(i-1));
Ordenada7mCAa=OrdenadaRp_ZCApua(1:1,1:(i-1));
Ordenada7mAGa=OrdenadaRp_ZAGpua(1:1,1:(i-1));
Ordenada7mBGa=OrdenadaRp_ZBGpua(1:1,1:(i-1));
Ordenada7mCGa=OrdenadaRp_ZCGpua(1:1,1:(i-1));

Ordenada10mABa=OrdenadaXp_ZABpua(1:1,1:(i-1));
Ordenada10mBCa=OrdenadaXp_ZBCpua(1:1,1:(i-1));
Ordenada10mCAa=OrdenadaXp_ZCApua(1:1,1:(i-1));
Ordenada10mAGa=OrdenadaXp_ZAGpua(1:1,1:(i-1));
Ordenada10mBGa=OrdenadaXp_ZBGpua(1:1,1:(i-1));
Ordenada10mCGa=OrdenadaXp_ZCGpua(1:1,1:(i-1));

% Tramo 4
% SX:m (S=Extremo primero X=Extremo segundo) 
% RX:m (R=Extremo primero X=Extremo segundo)

% amp: m=linea con falta m p=extremo primero  
% ams:                     s=extremo segundo  
% app: p=linea con falta p p=extremo primero  
% aps:                     s=extremo segundo 

amp=L1/Lma; bmp=L4/Lma; ams=L4/Lma; bms=-L4/Lma;
app=(L3+L2)/Lpa; bpp=L4/Lpa; aps=L4/Lpa; bps=-L4/Lpa;

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam0(cont)=mm;
%Ordenadam0(cont)=L1/(L4+L1)+mm*(L4/(L4+L1));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam00(cont)=mm;
%Ordenadam00(cont)=L4/(L4+L1)-mm*(L4/(L4+L1));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap0(cont)=mm;
%Ordenadap0(cont)=(L3+L2)/Lpa+mm*(L4/Lpa);
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap00(cont)=mm;
%Ordenadap00(cont)=L4/Lpa-mm*(L4/Lpa);
%cont=cont+1;
%mm=mm+0.01;
%end

end % end del switch(salida), para el tramo 4 

case{5}, %  ******** salidas para falta en el tramo 5 ********
   
%case{5}, % SY,RY
%Lma=L1+L2+L5;
%Lpa=L3+L5;

fprintf('\n');
fprintf('Salida Grafica 1: Extremos m: Todas las Unidades (pu)\n');
fprintf('                : Extremos p: Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2: Extremos m: Todas las Unidades (pua)\n');
fprintf('                : Extremos p: Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3: Analisis Extremo m Izq, 67N\n');
fprintf('                : Analisis Extremo m Der, 67N\n');
fprintf('                : Analisis Extremo p Izq, 67N\n');
fprintf('                : Analisis Extremo p Der, 67N\n');
fprintf('Salida Grafica 4: Analisis Extremo m Izq, R-X\n');
fprintf('                : Analisis Extremo m Der, R-X\n');
fprintf('                : Analisis Extremo p Izq, R-X\n');
fprintf('                : Analisis Extremo p Der, R-X\n');
fprintf('\n');

salida=input('Salida: 1,2,3,4 :');

switch(salida)
         
case{1}, % salida 1, para el tramo 5

cadenaym1=sprintf('alcanceSmpu');
cadenaym2=sprintf('alcanceYmpu');
cadenayp1=sprintf('alcanceRppu');
cadenayp2=sprintf('alcanceYppu');

Ordenada1mAB=OrdenadaSm_ZABpu(1:1,1:(i-1));
Ordenada1mBC=OrdenadaSm_ZBCpu(1:1,1:(i-1));
Ordenada1mCA=OrdenadaSm_ZCApu(1:1,1:(i-1));
Ordenada1mAG=OrdenadaSm_ZAGpu(1:1,1:(i-1));
Ordenada1mBG=OrdenadaSm_ZBGpu(1:1,1:(i-1));
Ordenada1mCG=OrdenadaSm_ZCGpu(1:1,1:(i-1));
   
Ordenada4mAB=OrdenadaYm_ZABpu(1:1,1:(i-1));
Ordenada4mBC=OrdenadaYm_ZBCpu(1:1,1:(i-1));
Ordenada4mCA=OrdenadaYm_ZCApu(1:1,1:(i-1));
Ordenada4mAG=OrdenadaYm_ZAGpu(1:1,1:(i-1));
Ordenada4mBG=OrdenadaYm_ZBGpu(1:1,1:(i-1));
Ordenada4mCG=OrdenadaYm_ZCGpu(1:1,1:(i-1));

Ordenada7mAB=OrdenadaRp_ZABpu(1:1,1:(i-1));
Ordenada7mBC=OrdenadaRp_ZBCpu(1:1,1:(i-1));
Ordenada7mCA=OrdenadaRp_ZCApu(1:1,1:(i-1));
Ordenada7mAG=OrdenadaRp_ZAGpu(1:1,1:(i-1));
Ordenada7mBG=OrdenadaRp_ZBGpu(1:1,1:(i-1));
Ordenada7mCG=OrdenadaRp_ZCGpu(1:1,1:(i-1));

Ordenada10mAB=OrdenadaYp_ZABpu(1:1,1:(i-1));
Ordenada10mBC=OrdenadaYp_ZBCpu(1:1,1:(i-1));
Ordenada10mCA=OrdenadaYp_ZCApu(1:1,1:(i-1));
Ordenada10mAG=OrdenadaYp_ZAGpu(1:1,1:(i-1));
Ordenada10mBG=OrdenadaYp_ZBGpu(1:1,1:(i-1));
Ordenada10mCG=OrdenadaYp_ZCGpu(1:1,1:(i-1));


case{2}, % salida 2, para el tramo 5
   
   
cadenaym1a=sprintf('alcanceSmpua');
cadenaym2a=sprintf('alcanceYmpua');
cadenayp1a=sprintf('alcanceRppua');
cadenayp2a=sprintf('alcanceYppua');
   
Ordenada1mABa=OrdenadaSm_ZABpua(1:1,1:(i-1));
Ordenada1mBCa=OrdenadaSm_ZBCpua(1:1,1:(i-1));
Ordenada1mCAa=OrdenadaSm_ZCApua(1:1,1:(i-1));
Ordenada1mAGa=OrdenadaSm_ZAGpua(1:1,1:(i-1));
Ordenada1mBGa=OrdenadaSm_ZBGpua(1:1,1:(i-1));
Ordenada1mCGa=OrdenadaSm_ZCGpua(1:1,1:(i-1));
          
Ordenada4mABa=OrdenadaYm_ZABpua(1:1,1:(i-1));
Ordenada4mBCa=OrdenadaYm_ZBCpua(1:1,1:(i-1));
Ordenada4mCAa=OrdenadaYm_ZCApua(1:1,1:(i-1));
Ordenada4mAGa=OrdenadaYm_ZAGpua(1:1,1:(i-1));
Ordenada4mBGa=OrdenadaYm_ZBGpua(1:1,1:(i-1));
Ordenada4mCGa=OrdenadaYm_ZCGpua(1:1,1:(i-1));

Ordenada7mABa=OrdenadaRp_ZABpua(1:1,1:(i-1));
Ordenada7mBCa=OrdenadaRp_ZBCpua(1:1,1:(i-1));
Ordenada7mCAa=OrdenadaRp_ZCApua(1:1,1:(i-1));
Ordenada7mAGa=OrdenadaRp_ZAGpua(1:1,1:(i-1));
Ordenada7mBGa=OrdenadaRp_ZBGpua(1:1,1:(i-1));
Ordenada7mCGa=OrdenadaRp_ZCGpua(1:1,1:(i-1));

Ordenada10mABa=OrdenadaYp_ZABpua(1:1,1:(i-1));
Ordenada10mBCa=OrdenadaYp_ZBCpua(1:1,1:(i-1));
Ordenada10mCAa=OrdenadaYp_ZCApua(1:1,1:(i-1));
Ordenada10mAGa=OrdenadaYp_ZAGpua(1:1,1:(i-1));
Ordenada10mBGa=OrdenadaYp_ZBGpua(1:1,1:(i-1));
Ordenada10mCGa=OrdenadaYp_ZCGpua(1:1,1:(i-1));

% Tramo 5
% SY:m (S=Extremo primero Y=Extremo segundo) 
% RY:p (R=Extremo primero Y=Extremo segundo)

% amp: m=linea con falta m p=extremo primero  
% ams:                     s=extremo segundo  
% app: p=linea con falta p p=extremo primero  
% aps:                     s=extremo segundo 

amp=(L1+L2)/Lma; bmp=L5/Lma; ams=L5/Lma; bms=-L5/Lma;
app=L3/Lpa; bpp=L5/Lpa; aps=L5/Lpa; bps=-L5/Lpa;

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam0(cont)=mm;
%Ordenadam0(cont)=(L1+L2)/(L1+L2+L5)+mm*(L5/(L1+L2+L5));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisam00(cont)=mm;
%Ordenadam00(cont)=L5/(L1+L2+L5)-mm*(L5/(L1+L2+L5));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap0(cont)=mm;
%Ordenadap0(cont)=L3/(L5+L3)+mm*(L5/(L5+L3));
%cont=cont+1;
%mm=mm+0.01;
%end

%mm=0.02;
%cont=1;
%while mm<1;
%Abscisap00(cont)=mm;
%Ordenadap00(cont)=L5/(L3+L5)-mm*(L5/(L3+L5));
%cont=cont+1;
%mm=mm+0.01;
%end

end % end del switch(salida), para el tramo 5 
end % end del switch(tramo)

% --------------------------------------------------------

if compens21==1
   comp='21c';
elseif compens21==2
   comp='21nc';
end
% --------------------------------------------------------
%***************************************
% Titulo para las Salidas 1,2,4
%***************************************
cadenam1=sprintf('PG54ic,%s,%3.0fkV,kE=%-4.2f,T%-1.0f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm,L4=%3.0fkm,L5=%3.0fkm',...
   comp,Un,kE,tramo,L1,L2,L3,L4,L5);
cadenam2=sprintf('Sccs=%-5.0f,Sccx=%-5.0f \n Sccr=%-5.0f,Sccy=%-5.0f',...
   Sccs,Sccx,Sccr,Sccy);
cadenam3=sprintf('Kx0=%-6.2f,Kr0=%-6.2f,Ky0=%-6.2f,k0=%-6.2f',...
   Kx0,Kr0,Ky0,k0);
cadenam4=sprintf('m: %s, Rfm=%-3.0f',cfaltam,Rfm);
if faltam==7|faltam==8|faltam==9|faltam==11
cadenam5=sprintf('Rgm=%-3.0f',Rgm);
cadenam4=sprintf('%s,%s',cadenam4,cadenam5);
end
cadenam=sprintf('%s,%s,%s,%s',cadenam1,cadenam2,cadenam3,cadenam4);

cadenap=sprintf('p:%s,Rfp=%-3.0f',cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
cadenap1=sprintf('Rgp=%-3.0f',Rgp);
cadenap=sprintf('%s,%s',cadenap,cadenap1);
end

cadenamp=sprintf('%s, %s',cadenam,cadenap);

%***************************************
% Titulo para la Salida 3, gráficas 67N
%***************************************
dircadenam1=sprintf('PG54ic,67N,%3.0fkV,kE=%-4.2f,T%-1.0f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm,L4=%3.0fkm,L5=%3.0fkm',...
   Un,kE,tramo,L1,L2,L3,L4,L5);
dircadenam2=sprintf('Sccs=%-5.0f,Sccx=%-5.0f \n Sccr=%-5.0f,Sccy=%-5.0f',...
   Sccs,Sccx,Sccr,Sccy);
dircadenam3=sprintf('Kx0=%-6.2f,Kr0=%-6.2f,Ky0=%-6.2f,k0=%-6.2f',...
   Kx0,Kr0,Ky0,k0);
dircadenam4=sprintf('m: %s, Rfm=%-3.0f',cfaltam,Rfm);
if faltam==7|faltam==8|faltam==9|faltam==11
dircadenam5=sprintf('Rgm=%-3.0f',Rgm);
dircadenam4=sprintf('%s,%s',dircadenam4,dircadenam5);
end
dircadenam=sprintf('%s,%s,%s,%s',dircadenam1,dircadenam2,dircadenam3,dircadenam4);

dircadenap=sprintf('p:%s,Rfp=%-3.0f',cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
dircadenap1=sprintf('Rgp=%-3.0f',Rgp);
dircadenap=sprintf('%s,%s',dircadenap,dircadenap1);
end
dircadenamp=sprintf('%s, %s',dircadenam,dircadenap);
%*************************************

switch(salida)
   
case{1}, % Salida 1, pu
   
mm=0.02;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
cont=cont+1;
mm=mm+0.01;
end

extremos=input('Extremos: m=1, p=2: ');
switch(extremos)
   
case{1}, % Salida 1, Extremos m 
   
fprintf('%s =1,%s =2\n',em1,em2);
extr=input('Extremo:');
switch(extr)
case{1}, % salida 1, Extremo m1 
   
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,Ordenada1mAB,'r',Abscisa,Ordenada1mBC,'g',Abscisa,Ordenada1mCA,'b',...
   Abscisa,Ordenada1mAG,'m',Abscisa,Ordenada1mBG,'k'),grid on
hold off

% Colores
hold on
a=plot(Abscisa,Ordenada1mCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores

axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenaym1)
legend('Ref','AB','BC','CA','AG','BG','CG')

case{2}, % salida 2, Extremos m2
   
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,Ordenada4mAB,'r',Abscisa,Ordenada4mBC,'g',Abscisa,Ordenada4mCA,'b',...
      Abscisa,Ordenada4mAG,'m',Abscisa,Ordenada4mBG,'k'),grid on
hold off
  
% Colores
hold on
a=plot(Abscisa,Ordenada4mCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
   
axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenaym2)
legend('Ref','AB','BC','CA','AG','BG','CG')
end

case{2}, % Salida 1, Extremos p
fprintf('%s =1,%s =2\n',ep1,ep2);
extr=input('Extremo:');

switch(extr)
case{1}, % salida 1, Extremo p1
   
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,Ordenada7mAB,'r',Abscisa,Ordenada7mBC,'g',Abscisa,Ordenada7mCA,'b',...
   Abscisa,Ordenada7mAG,'m',Abscisa,Ordenada7mBG,'k'),grid on
hold off

% Colores
hold on
a=plot(Abscisa,Ordenada7mCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores

axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenayp1)
legend('Ref','AB','BC','CA','AG','BG','CG')

case{2}, % salida 1, Extremo p2
   
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
% Fin Colores

hold on
   plot(Abscisa,Ordenada10mAB,'r',Abscisa,Ordenada10mBC,'g',Abscisa,Ordenada10mCA,'b',...
      Abscisa,Ordenada10mAG,'m',Abscisa,Ordenada10mBG,'k'),grid on
hold off
   
% Colores
hold on
a=plot(Abscisa,Ordenada10mCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
   
axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenayp2)
legend('Ref','AB','BC','CA','AG','BG','CG')
end
end

case{2}, % Salida 2
   
AnalisisLinea=input('Extremos Falta m=1, Falta p=2: ');
switch(AnalisisLinea)
case{1}, % salida 2, Extremos m
   
mm=0.02;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenadacmp(cont)=amp+bmp*mm;  
Ordenadacms(cont)=ams+bms*mm;
cont=cont+1;
mm=mm+0.01;
end

fprintf('%s =1,%s =2\n',em1,em2);
extr=input('Extremo:');
switch(extr)
       case{1}, % salida 2, Extremo m1 
          %Abscisam0,Ordenadam0,'y',
          
% Colores
b=plot(Abscisa0,Ordenadacmp);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,Ordenada1mABa,'r',Abscisa,Ordenada1mBCa,'g',Abscisa,Ordenada1mCAa,'b',...
   Abscisa,Ordenada1mAGa,'m',Abscisa,Ordenada1mBGa,'k'),grid on
hold off

% Colores
hold on
a=plot(Abscisa,Ordenada1mCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores

axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenaym1a)
legend('45º,Cob','AB','BC','CA','AG','BG','CG','z1')

case{2}, % salida 2, Extremos m2
   %Abscisam00,Ordenadam00,'y',
   
% Colores
b=plot(Abscisa0,Ordenadacms);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,Ordenada4mABa,'r',Abscisa,Ordenada4mBCa,'g',Abscisa,Ordenada4mCAa,'b',...
   Abscisa,Ordenada4mAGa,'m',Abscisa,Ordenada4mBGa,'k'),grid on
hold off

% Colores
hold on
a=plot(Abscisa,Ordenada4mCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores


axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenaym2a)
legend('-45º,Cob','AB','BC','CA','AG','BG','CG','z1')
end

case{2}, % salida 2, Extremos p
   
mm=0.02;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenadacpp(cont)=app+bpp*mm;
Ordenadacps(cont)=aps+bps*mm;
cont=cont+1;
mm=mm+0.01;
end

fprintf('%s =1,%s =2\n',ep1,ep2);
extr=input('Extremo: ');
switch(extr)
       case{1}, % salida 2, Extremo p1
          %Abscisap0,Ordenadap0,'y',
          
% Colores
b=plot(Abscisa0,Ordenadacpp);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,Ordenada7mABa,'r',Abscisa,Ordenada7mBCa,'g',Abscisa,Ordenada7mCAa,'b',...
   Abscisa,Ordenada7mAGa,'m',Abscisa,Ordenada7mBGa,'k'),grid on
hold off

% Colores
hold on
a=plot(Abscisa,Ordenada7mCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores

axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenayp1a)
legend('45º,Cob','AB','BC','CA','AG','BG','CG','z1')

case{2}, % salida 2, Extremo p2
   %Abscisap00,Ordenadap00,'y',
   
% Colores
b=plot(Abscisa0,Ordenadacps);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,Ordenada10mABa,'r',Abscisa,Ordenada10mBCa,'g',Abscisa,Ordenada10mCAa,'b',...
   Abscisa,Ordenada10mAGa,'m',Abscisa,Ordenada10mBGa,'k'),grid on
hold off

% Colores
hold on
a=plot(Abscisa,Ordenada10mCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores

axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenayp2a)
legend('-45º,Cob','AB','BC','CA','AG','BG','CG','z1')
end % switch(extr 1,2)
end % switch(extremos m,p)

%*************************************************
case{3}  % Salida Proteccion direccional
   
   AnalisisLinea=input('Extremos Falta m=1, Falta p=2: ');
   switch(AnalisisLinea)
   case{1}   % Extremos falta m
      
      fprintf('%s =1,%s =2\n',em1,em2);
      extr=input('Extremo:');
      switch(extr)
        case{1}, % salida 3, Extremo m1 

   subplot(3,1,1)
   plot(Abscisa,vectorAngDir_m_izq,'b');
   axis([0 1 -180 180]);
   xlabel('m=p'),ylabel(cadAngDir_m_izq)
   title(dircadenamp)
   rayar(1);
   legend(cadAngDir_m_izq,'AngMaxPar')
   
   subplot(3,1,2)
   plot(Abscisa,vectorVpol_m_izq,'b');
   grid
   %title(dircadenamp)
   xlabel('m=p'),ylabel(cadVpol_m_izq)

   subplot(3,1,3)
   plot(Abscisa,vectorIop_m_izq,'b');
   grid
   xlabel('m=p'),ylabel(cadIop_m_izq)
   
case{2}  % salida 3, Extremo m2 
   
   subplot(3,1,1)
   plot(Abscisa,vectorAngDir_m_der,'b');
   axis([0 1 -180 180]);
   xlabel('m=p'),ylabel(cadAngDir_m_der)
   title(dircadenamp)
   rayar(1);
   legend(cadAngDir_m_der,'AngMaxPar')
   
   subplot(3,1,2)
   plot(Abscisa,vectorVpol_m_der,'b');
   grid
   %title(dircadenamp)
   xlabel('m=p'),ylabel(cadVpol_m_der)

   subplot(3,1,3)
   plot(Abscisa,vectorIop_m_der,'b');
   grid
   xlabel('m=p'),ylabel(cadIop_m_der)

end % end del switch extr

   case{2}   % Extremos falta p
      if tramo==1 | tramo==2 | tramo==3
      fprintf('%s =1,%s =2\n',ep1,ep2);
      extr=input('Extremo:');
      elseif tramo==4 | tramo==5
      fprintf('%s =1,%s =2\n',ep2,ep1);
      extr=input('Extremo:');
      end
      switch(extr)
        case{1}, % salida 3, Extremo p1 

   subplot(3,1,1)
   plot(Abscisa,vectorAngDir_p_izq,'b');
   axis([0 1 -180 180]);
   xlabel('m=p'),ylabel(cadAngDir_p_izq)
   title(dircadenamp)
   rayar(1);
   legend(cadAngDir_p_izq,'AngMaxPar')
   
   subplot(3,1,2)
   plot(Abscisa,vectorVpol_p_izq,'b');
   grid
   %title(dircadenamp)
   xlabel('m=p'),ylabel(cadVpol_p_izq)

   subplot(3,1,3)
   plot(Abscisa,vectorIop_p_izq,'b');
   grid
   xlabel('m=p'),ylabel(cadIop_p_izq)
   
case{2}  % salida 3, Extremo p2 
   
   subplot(3,1,1)
   plot(Abscisa,vectorAngDir_p_der,'b');
   axis([0 1 -180 180]);
   xlabel('m=p'),ylabel(cadAngDir_p_der)
   title(dircadenamp)
   rayar(1);
   legend(cadAngDir_p_der,'AngMaxPar')
   
   subplot(3,1,2)
   plot(Abscisa,vectorVpol_p_der,'b');
   grid
   %title(dircadenamp)
   xlabel('m=p'),ylabel(cadVpol_p_der)

   subplot(3,1,3)
   plot(Abscisa,vectorIop_p_der,'b');
   grid
   xlabel('m=p'),ylabel(cadIop_p_der)
   
   
end % end del switch extr
end % end del switch Analisis Linea

%*************************************************
case{4},  % *********  Salida 4  R-X  ************
%*************************************************
   
fprintf('Extremo: Analisis Extremo %s=1\n',m1);
fprintf('         Analisis Extremo %s=2\n',m2);
fprintf('         Analisis Extremo %s=3\n',p1);
fprintf('         Analisis Extremo %s=4\n',p2);
   
extr=input('Extremo: ');
   
         switch(extr)
         case{1}
extremo='m1';          
proteccion=m1; 
L=Lma;
pzona=0.8;
   hold on 
   plot19c(L,Lt,Latm,Lttm,Latp,Lttp,zl1u,pzona,proteccion,extremo,tramo,L1,L2,L3,L4,L5)
   hold off
   hold on
   % Colores
   a=plot(vectorZAB_m1);
   set(a,'Color',[1 0.60 0])
   % Fin Colores
   plot(vectorZBC_m1,'m')
   plot(vectorZCA_m1,'c')
   plot(vectorZAG_m1,'r')
   plot(vectorZBG_m1,'k')
   % Colores
   c=plot(vectorZCG_m1);
   set(c,'Color',[0.7216 0.4510 0.20])
   % Fin Colores
   hold off

   [subAbscisa,subvectorZAB_m1]=subconjunto21(Abscisa,vectorZAB_m1);
             subAbscisa(:,1)=[];
             subvectorZAB_m1(:,1)=[];
             hold on
             b=plot(subvectorZAB_m1,'+');
             set(b,'Color',[1 0.60 0])
             plot20(L,zl1u,subvectorZAB_m1,'naranja')
             hold off
   [subAbscisa,subvectorZBC_m1]=subconjunto21(Abscisa,vectorZBC_m1);
             subAbscisa(:,1)=[];
             subvectorZBC_m1(:,1)=[];
             hold on
             plot(subvectorZBC_m1,'m+');
             plot20(L,zl1u,subvectorZBC_m1,'m')
             hold off
    [subAbscisa,subvectorZCA_m1]=subconjunto21(Abscisa,vectorZCA_m1);
             subAbscisa(:,1)=[];
             subvectorZCA_m1(:,1)=[];
             hold on
             plot(subvectorZCA_m1,'c+');
             plot20(L,zl1u,subvectorZCA_m1,'c')
             hold off
   [subAbscisa,subvectorZAG_m1]=subconjunto21(Abscisa,vectorZAG_m1);
             subAbscisa(:,1)=[];
             subvectorZAG_m1(:,1)=[];
             hold on
             plot(subvectorZAG_m1,'r+');
             plot20(L,zl1u,subvectorZAG_m1,'r')
             hold off
     [subAbscisa,subvectorZBG_m1]=subconjunto21(Abscisa,vectorZBG_m1);
             subAbscisa(:,1)=[];
             subvectorZBG_m1(:,1)=[];
             hold on
             plot(subvectorZBG_m1,'k+');
             plot20(L,zl1u,subvectorZBG_m1,'k')
             hold off
   [subAbscisa,subvectorZCG_m1]=subconjunto21(Abscisa,vectorZCG_m1);
             subAbscisa(:,1)=[];
             subvectorZCG_m1(:,1)=[];
             hold on
             % Colores
             c=plot(subvectorZCG_m1,'+');
             set(c,'Color',[0.7216 0.4510 0.20])
             % Fin Colores
             plot20(L,zl1u,subvectorZCG_m1,'marron')
             hold off
             
plot21(L,zl1u,extremo,tramo,L1,L2,L3,L4,L5)
             
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
            
title(RXcadenamp)

legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG')

         case{2}
            L=Lma;
extremo='m2';          
proteccion=m2;
pzona=0.8;
   hold on 
   plot19c(L,Lt,Latm,Lttm,Latp,Lttp,zl1u,pzona,proteccion,extremo,tramo,L1,L2,L3,L4,L5)
   hold off
   hold on
   % Colores
   a=plot(vectorZAB_m2);
   set(a,'Color',[1 0.60 0])
   % Fin Colores
   plot(vectorZBC_m2,'m')
   plot(vectorZCA_m2,'c')
   plot(vectorZAG_m2,'r')
   plot(vectorZBG_m2,'k')
   % Colores
   c=plot(vectorZCG_m2);
   set(c,'Color',[0.7216 0.4510 0.20])
   % Fin Colores
   hold off
   
   [subAbscisa,subvectorZAB_m2]=subconjunto21(Abscisa,vectorZAB_m2);
             subAbscisa(:,11)=[];
             subvectorZAB_m2(:,11)=[];
             invsubvectorZAB_m2=subvectorZAB_m2(1:1,10:-1:1);
             hold on
             b=plot(invsubvectorZAB_m2,'+');
             set(b,'Color',[1 0.60 0])
             plot20(L,zl1u,invsubvectorZAB_m2,'naranja')
             hold off
   [subAbscisa,subvectorZBC_m2]=subconjunto21(Abscisa,vectorZBC_m2);
             subAbscisa(:,11)=[];
             subvectorZBC_m2(:,11)=[];
             invsubvectorZBC_m2=subvectorZBC_m2(1:1,10:-1:1);
             hold on
             plot(subvectorZBC_m2,'m+');
             plot20(L,zl1u,invsubvectorZBC_m2,'m')
             hold off
    [subAbscisa,subvectorZCA_m2]=subconjunto21(Abscisa,vectorZCA_m2);
             subAbscisa(:,11)=[];
             subvectorZCA_m2(:,11)=[];
             invsubvectorZCA_m2=subvectorZCA_m2(1:1,10:-1:1);
             hold on
             plot(subvectorZCA_m2,'c+');
             plot20(L,zl1u,invsubvectorZCA_m2,'c')
             hold off
    [subAbscisa,subvectorZAG_m2]=subconjunto21(Abscisa,vectorZAG_m2);
             subAbscisa(:,11)=[];
             subvectorZAG_m2(:,11)=[];
             invsubvectorZAG_m2=subvectorZAG_m2(1:1,10:-1:1);
             hold on
             plot(subvectorZAG_m2,'r+');
             plot20(L,zl1u,invsubvectorZAG_m2,'r')
             hold off
    [subAbscisa,subvectorZBG_m2]=subconjunto21(Abscisa,vectorZBG_m2);
             subAbscisa(:,11)=[];
             subvectorZBG_m2(:,11)=[];
             invsubvectorZBG_m2=subvectorZBG_m2(1:1,10:-1:1);
             hold on
             plot(subvectorZBG_m2,'k+');
             plot20(L,zl1u,invsubvectorZBG_m2,'k')
             hold off
    [subAbscisa,subvectorZCG_m2]=subconjunto21(Abscisa,vectorZCG_m2);
             subAbscisa(:,11)=[];
             subvectorZCG_m2(:,11)=[];
             invsubvectorZCG_m2=subvectorZCG_m2(1:1,10:-1:1);
             hold on
             % Colores
             c=plot(invsubvectorZCG_m2,'+');
             set(c,'Color',[0.7216 0.4510 0.20])
             % Fin Colores
             plot20(L,zl1u,invsubvectorZCG_m2,'marron')
             hold off
             
plot21(L,zl1u,extremo,tramo,L1,L2,L3,L4,L5)
             
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
            
title(RXcadenamp)

legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG')

case{3}
L=Lpa;   
extremo='p1';          
proteccion=p1;
pzona=0.8;
   hold on 
   plot19c(L,Lt,Latm,Lttm,Latp,Lttp,zl1u,pzona,proteccion,extremo,tramo,L1,L2,L3,L4,L5)
   hold off
   hold on
   % Colores
   a=plot(vectorZAB_p1);
   set(a,'Color',[1 0.60 0])
   % Fin Colores
   plot(vectorZBC_p1,'m')
   plot(vectorZCA_p1,'c')
   plot(vectorZAG_p1,'r')
   plot(vectorZBG_p1,'k')
   % Colores
   c=plot(vectorZCG_p1);
   set(c,'Color',[0.7216 0.4510 0.20])
   % Fin Colores
   hold off
   
   [subAbscisa,subvectorZAB_p1]=subconjunto21(Abscisa,vectorZAB_p1);
             subAbscisa(:,1)=[];
             subvectorZAB_p1(:,1)=[];
             hold on
             b=plot(subvectorZAB_p1,'+');
             set(b,'Color',[1 0.60 0])
             plot20(L,zl1u,subvectorZAB_p1,'naranja')
             hold off
   [subAbscisa,subvectorZBC_p1]=subconjunto21(Abscisa,vectorZBC_p1);
             subAbscisa(:,1)=[];
             subvectorZBC_p1(:,1)=[];
             hold on
             plot(subvectorZBC_p1,'m+');
             plot20(L,zl1u,subvectorZBC_p1,'m')
             hold off
    [subAbscisa,subvectorZCA_p1]=subconjunto21(Abscisa,vectorZCA_p1);
             subAbscisa(:,1)=[];
             subvectorZCA_p1(:,1)=[];
             hold on
             plot(subvectorZCA_p1,'c+');
             plot20(L,zl1u,subvectorZCA_p1,'c')
             hold off
   [subAbscisa,subvectorZAG_p1]=subconjunto21(Abscisa,vectorZAG_p1);
             subAbscisa(:,1)=[];
             subvectorZAG_p1(:,1)=[];
             hold on
             plot(subvectorZAG_p1,'r+');
             plot20(L,zl1u,subvectorZAG_p1,'r')
             hold off
     [subAbscisa,subvectorZBG_p1]=subconjunto21(Abscisa,vectorZBG_p1);
             subAbscisa(:,1)=[];
             subvectorZBG_p1(:,1)=[];
             hold on
             plot(subvectorZBG_p1,'k+');
             plot20(L,zl1u,subvectorZBG_p1,'k')
             hold off
   [subAbscisa,subvectorZCG_p1]=subconjunto21(Abscisa,vectorZCG_p1);
             subAbscisa(:,1)=[];
             subvectorZCG_p1(:,1)=[];
             hold on
             % Colores
             c=plot(subvectorZCG_p1,'+');
             set(c,'Color',[0.7216 0.4510 0.20])
             % Fin Colores
             plot20(L,zl1u,subvectorZCG_p1,'marron')
             hold off
             
plot21(L,zl1u,extremo,tramo,L1,L2,L3,L4,L5)
             
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
            
title(RXcadenamp)

legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG')

         case{4}
L=Lpa;
extremo='p2';          
proteccion=p2;
pzona=0.8;
   hold on 
   plot19c(L,Lt,Latm,Lttm,Latp,Lttp,zl1u,pzona,proteccion,extremo,tramo,L1,L2,L3,L4,L5)
   hold off
   hold on
   % Colores
   a=plot(vectorZAB_p2);
   set(a,'Color',[1 0.60 0])
   % Fin Colores
   plot(vectorZBC_p2,'m')
   plot(vectorZCA_p2,'c')
   plot(vectorZAG_p2,'r')
   plot(vectorZBG_p2,'k')
   % Colores
   c=plot(vectorZCG_p2);
   set(c,'Color',[0.7216 0.4510 0.20])
   % Fin Colores
   hold off
   
   [subAbscisa,subvectorZAB_p2]=subconjunto21(Abscisa,vectorZAB_p2);
             subAbscisa(:,11)=[];
             subvectorZAB_p2(:,11)=[];
             invsubvectorZAB_p2=subvectorZAB_p2(1:1,10:-1:1);
             hold on
             b=plot(invsubvectorZAB_p2,'+');
             set(b,'Color',[1 0.60 0])
             plot20(L,zl1u,invsubvectorZAB_p2,'naranja')
             hold off
   [subAbscisa,subvectorZBC_p2]=subconjunto21(Abscisa,vectorZBC_p2);
             subAbscisa(:,11)=[];
             subvectorZBC_p2(:,11)=[];
             invsubvectorZBC_p2=subvectorZBC_p2(1:1,10:-1:1);
             hold on
             plot(subvectorZBC_p2,'m+');
             plot20(L,zl1u,invsubvectorZBC_p2,'m')
             hold off
    [subAbscisa,subvectorZCA_p2]=subconjunto21(Abscisa,vectorZCA_p2);
             subAbscisa(:,11)=[];
             subvectorZCA_p2(:,11)=[];
             invsubvectorZCA_p2=subvectorZCA_p2(1:1,10:-1:1);
             hold on
             plot(subvectorZCA_p2,'c+');
             plot20(L,zl1u,invsubvectorZCA_p2,'c')
             hold off
    [subAbscisa,subvectorZAG_p2]=subconjunto21(Abscisa,vectorZAG_p2);
             subAbscisa(:,11)=[];
             subvectorZAG_p2(:,11)=[];
             invsubvectorZAG_p2=subvectorZAG_p2(1:1,10:-1:1);
             hold on
             plot(subvectorZAG_p2,'r+');
             plot20(L,zl1u,invsubvectorZAG_p2,'r')
             hold off
    [subAbscisa,subvectorZBG_p2]=subconjunto21(Abscisa,vectorZBG_p2);
             subAbscisa(:,11)=[];
             subvectorZBG_p2(:,11)=[];
             invsubvectorZBG_p2=subvectorZBG_p2(1:1,10:-1:1);
             hold on
             plot(subvectorZBG_p2,'k+');
             plot20(L,zl1u,invsubvectorZBG_p2,'k')
             hold off
    [subAbscisa,subvectorZCG_p2]=subconjunto21(Abscisa,vectorZCG_p2);
             subAbscisa(:,11)=[];
             subvectorZCG_p2(:,11)=[];
             invsubvectorZCG_p2=subvectorZCG_p2(1:1,10:-1:1);
             hold on
             % Colores
             c=plot(invsubvectorZCG_p2,'+');
             set(c,'Color',[0.7216 0.4510 0.20])
             % Fin Colores
             plot20(L,zl1u,invsubvectorZCG_p2,'marron')
             hold off
             
plot21(L,zl1u,extremo,tramo,L1,L2,L3,L4,L5)
             
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
            
title(RXcadenamp)

legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG')

end % end del switch extr de Salida 4

end % switch(salida 1,2,3,4)

variassalidas=input('Mas Salidas=1,Acabar=2 :');
end  % end del while variassalidas

variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :');
end  % end del while variasejecuciones
fprintf(1,'Pulse Enter\n');


<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="es">
	<head>
            <title>e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.</title>

                <meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>
<meta name="description" content="Buscador de proyectos ETSI" />
<meta name="author" content="Daniel Callejo"/>
<meta name="author" content="Federico L&oacute;pez"/>
<meta name="author" content="Juli&aacute;n P&eacute;rez"/>
<meta name="author" content="Jos&eacute; Mar&iacute;a Vidal Vidal"/>
<base href="https://biblus.us.es/bibing/proyectos/" /><link rel="shortcut icon" href="./favicon.ico" />
<link rel="stylesheet" href="./publico/css/estructura.css" type="text/css" />
<link rel="stylesheet" href="./publico/css/bootstrap.css" type="text/css" media="screen"/>
<link rel="stylesheet" href="./publico/css/enlaces-capas.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/pestanas.css" type="text/css"  /> 
<link rel="stylesheet" href="./publico/css/formularios.css" type="text/css" />
<link rel="stylesheet" href="./publico/css/lightbox.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/autocompleter.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/resultados.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/milista.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/explorador.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/ultimos.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/contactar.css" type="text/css"  />
<!-- <link rel="stylesheet" href="./css/print.css" type="text/css" media="print" /> -->
<script src="./publico/javascript/jquery-1.9.1.min.js" type="text/javascript"></script>
<script src="./publico/javascript/bootstrap.js" type="text/javascript"></script>

<script src="./publico/javascript/mootools-1.2.4-core.js" type="text/javascript"></script>
<script src="./publico/javascript/mootools-1-1.2.4.2-more.js" type="text/javascript"></script>
<script src="./publico/javascript/targetbl.js" type="text/javascript"></script>
<script src="./publico/javascript/mensajes.js" type="text/javascript"></script>
<script src="./publico/javascript/form.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Observer.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Autocompleter.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Autocompleter.Request.js" type="text/javascript"></script>
<script src="./publico/javascript/domready.js" type="text/javascript"></script>

<!-- Google tag (gtag.js) -->
<script async src="https://www.googletagmanager.com/gtag/js?id=G-2HP38G406P"></script>
<script>
  window.dataLayer = window.dataLayer || [];
  function gtag(){dataLayer.push(arguments);}
  gtag('js', new Date());

  gtag('config', 'G-2HP38G406P');
</script>



<script type="text/javascript">
function seleccionar_todo(){
   for (i=0;i<document.formulario_marcados.elements.length;i++)
      if(document.formulario_marcados.elements[i].type == "checkbox")
         document.formulario_marcados.elements[i].checked=1
}
function deseleccionar_todo(){
   for (i=0;i<document.formulario_marcados.elements.length;i++)
      if(document.formulario_marcados.elements[i].type == "checkbox")
         document.formulario_marcados.elements[i].checked=0
} 
</script>
	</head>
    <body style="min-width: 960px; background-color: #e2e2e2; height: 100%">	
		<a name="arriba"></a>
<div>
<div id="header">
    <table id="Tabla_01" height="105" border="0" cellpadding="0" cellspacing="0">
    <tr>
        <td>
            <img src="./publico/imagenes/cabecera/banner_new2_01.jpg" width="3" height="105" alt="">
        </td>
        <td>
            <img src="./publico/imagenes/cabecera/banner_new2_02.jpg" width="100" height="105" alt="">
        </td>
        <td width="100%" align="right"  style="min-width: 613px; background-repeat: no-repeat; background-position: top right; font-size: medium; font-weight: bold; color: white;">
            <table width="100%" height="105" border="0">
                <tr>
                    <td style=" vertical-align: top; width: 547px">
                        <!-- <a style="text-decoration: none" href="./"><img src="./publico/imagenes/cabecera/banner_new2_03.jpg"alt=""></a> -->
                        <a style="text-decoration: none" href="https://biblus.us.es/bibing/proyectos/"><img src="./publico/imagenes/cabecera/banner_new2_03.jpg"alt=""></a>
                     </td>
                     <td  background="./publico/imagenes/cabecera/banner_new2_04.jpg" style="background-repeat: no-repeat; background-color: #c6d4e8;"><br><br><br>
                         &nbsp;
                        </td>
                    </tr>
             </table>
        </td>
    </tr>
    <tr style="border-top: 1px solid #ffffff;">
        <td colspan="3" id="menu_sup" style="font-size: small; text-align: right; background-color: #9c1f2f; vertical-align: bottom">
                                     

                                <a href="https://biblus.us.es/bibing/proyectos/" class="enlace " style=" color: #ffffff">Inicio&nbsp;&nbsp;|</a>
                                <a href="https://biblus.us.es/bibing/proyectos/acerca/" style=" color: #ffffff" class="enlace ">&nbsp;&nbsp;Acerca de&nbsp;&nbsp;|</a>
                                
                                <a href="https://biblus.us.es/bibing/proyectos/contactar/" style=" color: #ffffff" class="enlace ">&nbsp;&nbsp;Contactar&nbsp;&nbsp;</a>
                         <!--       <a href="https://biblus.us.es/bibing/proyectos/ayuda/" style=" color: #ffffff" class="enlace  ">&nbsp;&nbsp;Ayuda&nbsp;&nbsp;|</a> -->
                                <!-- <form id="idioma" class=""  method="post" action="">
                                <div class="inline" style="color: white">
                                <input type="submit" class="tipo_letra " style=" color: #ffffff"name="idioma" value="Castellanossss" />&nbsp;&nbsp;
                                </div>	
                                </form> -->

                                </td>

    </tr>
    </table>
</div>
<br><br>
<div class="colmask rightmenu">
<div class="panel panel-default" style="margin:20px">
            <div class="panel-heading"><h4>AN�LISIS DE ACTUACIÓN DE LAS PROTECCIONES 21 Y 67N ANTE FALTAS SIMPLES, MULTIPLES E INTERCIRCUITO EN L�NEAS ELÉCTRICAS SOBRE LOS MISMOS APOYOS</h4><h5 style="color:grey"><span class="glyphicon glyphicon-user""></span> : <b>Ruiz Lozano, José Luis</b><br><span class="glyphicon glyphicon-book""></span>  :<b> Ingenier&iacute;a Industrial</b></h5></div><div class="panel-body"><div class="contenido_proyecto"><div class="titulo2_explorador">Contenido del proyecto:</div><div class="ruta"><a href="use/abreproy/4273">Directorio ra&iacute;z</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252F">z_Prog</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPGc%252F">PGc</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPGc%252FPG54ic2%252F">PG54ic2</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPGc%252FPG54ic2%252FFG54ic2.m%252F">FG54ic2.m</a></div><div class="listing"></div></div></div></div>
</div>
</div>

        
	</body>
 
 
<!-- Inicio Google Analytics -->
<!--
	<script type="text/javascript">
		var gaJsHost = (("https:" == document.location.protocol) ? "https://ssl." : "http://www.");
		document.write(unescape("%3Cscript src='" + gaJsHost + "google-analytics.com/ga.js' type='text/javascript'%3E%3C/script%3E"));
	</script>
	
	<script type="text/javascript">
		try {
		var pageTracker = _gat._getTracker("UA-7789620-1");
		pageTracker._trackPageview();
	} 
		catch(err) {}
	</script>
-->
<!-- Fin Google Analytics -->

</html>