% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG54ni\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('Compensacion 21: Si,No\n');
fprintf('Falta Doble. Para cada Tramo: m=p=0...1\n');
fprintf('Eleccion: kE,L1,L2,L3,L4,L5,Sccs,Sccx,Sccr,Sccy\n');
fprintf('          Kx0,Kr0,Ky0\n');
fprintf('*** Parametro k0=zs0/zr0=0,1-1-10 ***\n');
fprintf('Para cada Tramo: Barra m=1, Barra p=2\n');
fprintf('Salida Grafica 1: Extr.m(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Extr.m(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
fprintf('\n');
% ---------------------------------------------------------------
%PG54ni
datos=input('datos: standard=1,input=2: ');
if datos==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;
% ---------------------------------------------------------------
else
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
% ---------------------------------------------------------------
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
% ---------------------------------------------------------------
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;   
Km=zlM0u/zl1u;
fprintf('\n');
% ---------------------------------------------------------------
pfaltam=1;
pfaltap=2;
% ---------------------------------------------------------------
fila=1;
k0=0.1;
while k0<11
zs0=k0*zr0;
i=1;
m=0.01;

while m<1
p=m;

switch(tramo)
case{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},
[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},
[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},
[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},
[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
   
% 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);

% PROTECCION DE DISTANCIA (21)
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);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
  PDistComp(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k,Km,Ism0);
[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=...
   PDistComp(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k,Km,Ixr0);
[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 21 
%
% 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

if faltam==1
   cfaltam='ag';
   elementom='AG';
   alcanceSm=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u));
   alcanceXm=abs(Xm_ZAG)/cos(angle(Xm_ZAG)-angle(zl1u));

elseif faltam==2
   cfaltam='bg';
   elementom='BG';
   alcanceSm=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u));
   alcanceXm=abs(Xm_ZBG)/cos(angle(Xm_ZBG)-angle(zl1u));
   
elseif faltam==3
   cfaltam='cg';
   elementom='CG';
   alcanceSm=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u));
   alcanceXm=abs(Xm_ZCG)/cos(angle(Xm_ZCG)-angle(zl1u));

elseif faltam==4
   cfaltam='ab';
   elementom='AB';
   alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
   alcanceXm=abs(Xm_ZAB)/cos(angle(Xm_ZAB)-angle(zl1u));

elseif faltam==5
   cfaltam='bc';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));

elseif faltam==6
   cfaltam='ca';
   elementom='CA';
   alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
   alcanceXm=abs(Xm_ZCA)/cos(angle(Xm_ZCA)-angle(zl1u));

elseif faltam==7
   cfaltam='abg';
   elementom='AB';
   alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
   alcanceXm=abs(Xm_ZAB)/cos(angle(Xm_ZAB)-angle(zl1u));

elseif faltam==8
   cfaltam='bcg';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));

elseif faltam==9
   cfaltam='cag';
   elementom='CA';
   alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
   alcanceXm=abs(Xm_ZCA)/cos(angle(Xm_ZCA)-angle(zl1u));

elseif faltam==10
   cfaltam='abc';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));

elseif faltam==11
   cfaltam='abcg';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));

end

% Recuerdo: tramo1

if faltap==1
   cfaltap='ag';
   elementop='AG';
   alcanceSp=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u));
   alcanceYp=abs(Yp_ZAG)/cos(angle(Yp_ZAG)-angle(zl1u));

elseif faltap==2
   cfaltap='bg';
   elementop='BG';
   alcanceSp=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u));
   alcanceYp=abs(Yp_ZBG)/cos(angle(Yp_ZBG)-angle(zl1u));
   
elseif faltap==3
   cfaltap='cg';
   elementop='CG';
   alcanceSp=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u));
   alcanceYp=abs(Yp_ZCG)/cos(angle(Yp_ZCG)-angle(zl1u));
   
elseif faltap==4
   cfaltap='ab';
   elementop='AB';
   alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
   alcanceYp=abs(Yp_ZAB)/cos(angle(Yp_ZAB)-angle(zl1u));
   
elseif faltap==5
   cfaltap='bc';
   elementop='BC';
   alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   
elseif faltap==6
   cfaltap='ca';
   elementop='CA';
   alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
   alcanceYp=abs(Yp_ZCA)/cos(angle(Yp_ZCA)-angle(zl1u));
   
elseif faltap==7
   cfaltap='abg';
   elementop='AB';
alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcanceYp=abs(Yp_ZAB)/cos(angle(Yp_ZAB)-angle(zl1u));

elseif faltap==8
   cfaltap='bcg';
   elementop='BC';
   alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   
elseif faltap==9
   cfaltap='cag';
   elementop='CA';
   alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
   alcanceYp=abs(Yp_ZCA)/cos(angle(Yp_ZCA)-angle(zl1u));
   
elseif faltap==10
   cfaltap='abc';
   elementop='BC';
   alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   
elseif faltap==11
   cfaltap='abcg';
   elementop='BC';
   alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   
end

alcanceSmpu=alcanceSm/abs(zl1u*m*L1);
alcanceXmpu=alcanceXm/abs(zl1u*(L4+(1-m)*L1));
alcanceSppu=alcanceSp/abs(zl1u*p*L1);
alcanceYppu=alcanceYp/abs(zl1u*(L5+L2+(1-p)*L1));

LimSm=(0.8*(L1+L4))/(m*L1);
LimXm=(0.8*(L1+L4))/(L4+(1-m)*L1);
LimSp=(0.8*(L1+L2+L5))/(m*L1);
LimYp=(0.8*(L1+L2+L5))/(L5+L2+(1-m)*L1);

vectorLimSm(1,i)=LimSm;
vectorLimXm(1,i)=LimXm;
vectorLimSp(1,i)=LimSp;
vectorLimYp(1,i)=LimYp;

Abscisa(fila,i)=m;
OrdenadaSm(fila,i)=alcanceSmpu;
OrdenadaXm(fila,i)=alcanceXmpu;
OrdenadaSp(fila,i)=alcanceSppu;
OrdenadaYp(fila,i)=alcanceYppu;

cadenaym1=sprintf('alcanceSmpu');
cadenaym2=sprintf('alcanceXmpu');
cadenayp1=sprintf('alcanceSppu');
cadenayp2=sprintf('alcanceYppu');

cadenaym1a=sprintf('alcanceSmpua');
cadenaym2a=sprintf('alcanceXmpua');
cadenayp1a=sprintf('alcanceSppua');
cadenayp2a=sprintf('alcanceYppua');

alcanceSmpua=alcanceSm/abs(zl1u*(L1+L4));
alcanceXmpua=alcanceXm/abs(zl1u*(L4+L1));
alcanceSppua=alcanceSp/abs(zl1u*(L1+L2+L5));
alcanceYppua=alcanceYp/abs(zl1u*(L5+L2+L1));

OrdenadaSma(fila,i)=alcanceSmpua;
OrdenadaXma(fila,i)=alcanceXmpua;
OrdenadaSpa(fila,i)=alcanceSppua;
OrdenadaYpa(fila,i)=alcanceYppua;

case{2},% tramo2
   
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);

% PROTECCION DE DISTANCIA (21)
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: tramo2

if faltam==1
   cfaltam='ag';
   elementom='AG';
   alcanceXm=abs(Xm_ZAG)/cos(angle(Xm_ZAG)-angle(zl1u));
   alcanceRm=abs(Rm_ZAG)/cos(angle(Rm_ZAG)-angle(zl1u));
   
elseif faltam==2
   cfaltam='bg';
   elementom='BG';
   alcanceXm=abs(Xm_ZBG)/cos(angle(Xm_ZBG)-angle(zl1u));
   alcanceRm=abs(Rm_ZBG)/cos(angle(Rm_ZBG)-angle(zl1u));
   
elseif faltam==3
   cfaltam='cg';
   elementom='CG';
   alcanceXm=abs(Xm_ZCG)/cos(angle(Xm_ZCG)-angle(zl1u));
   alcanceRm=abs(Rm_ZCG)/cos(angle(Rm_ZCG)-angle(zl1u));


elseif faltam==4
   cfaltam='ab';
   elementom='AB';
   alcanceXm=abs(Xm_ZAB)/cos(angle(Xm_ZAB)-angle(zl1u));
   alcanceRm=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));

elseif faltam==5
   cfaltam='bc';
   elementom='BC';
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));
   alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));

elseif faltam==6
   cfaltam='ca';
   elementom='CA';
   alcanceXm=abs(Xm_ZCA)/cos(angle(Xm_ZCA)-angle(zl1u));
   alcanceRm=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));

elseif faltam==7
   cfaltam='abg';
   elementom='AB';
   alcanceXm=abs(Xm_ZAB)/cos(angle(Xm_ZAB)-angle(zl1u));
   alcanceRm=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));

elseif faltam==8

   cfaltam='bcg';
   elementom='BC';
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));
   alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
     
elseif faltam==9
   cfaltam='cag';
   elementom='CA';
   alcanceXm=abs(Xm_ZCA)/cos(angle(Xm_ZCA)-angle(zl1u));
   alcanceRm=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
   
elseif faltam==10
   cfaltam='abc';
   elementom='BC';
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));
   alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
   
elseif faltam==11
   cfaltam='abcg';
   elementom='BC';
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));
   alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
   
end



   
% Recuerdo: tramo2
   
if faltap==1
   cfaltap='ag';
   elementop='AG';
   alcanceSp=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u));
   alcanceYp=abs(Yp_ZAG)/cos(angle(Yp_ZAG)-angle(zl1u));

elseif faltap==2
   cfaltap='bg';
   elementop='BG';
   alcanceSp=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u));
   alcanceYp=abs(Yp_ZBG)/cos(angle(Yp_ZBG)-angle(zl1u));
   
elseif faltap==3
   cfaltap='cg';
   elementop='CG';
   alcanceSp=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u));
   alcanceYp=abs(Yp_ZCG)/cos(angle(Yp_ZCG)-angle(zl1u));
   
elseif faltap==4
   cfaltap='ab';
   elementop='AB';
   alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
   alcanceYp=abs(Yp_ZAB)/cos(angle(Yp_ZAB)-angle(zl1u));
   
elseif faltap==5
   cfaltap='bc';
   elementop='BC';
   alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   
elseif faltap==6
   cfaltap='ca';
   elementop='CA';
   alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
   alcanceYp=abs(Yp_ZCA)/cos(angle(Yp_ZCA)-angle(zl1u));
   
elseif faltap==7
   cfaltap='abg';

   elementop='AB';
alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcanceYp=abs(Yp_ZAB)/cos(angle(Yp_ZAB)-angle(zl1u));

elseif faltap==8
   cfaltap='bcg';
   elementop='BC';
   alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   
elseif faltap==9
   cfaltap='cag';
   elementop='CA';
   alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
   alcanceYp=abs(Yp_ZCA)/cos(angle(Yp_ZCA)-angle(zl1u));
   
elseif faltap==10
   cfaltap='abc';
   elementop='BC';
   alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   
elseif faltap==11
   cfaltap='abcg';
   elementop='BC';
   alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   
end

alcanceXmpu=alcanceXm/abs(zl1u*(L4+m*L2));
alcanceRmpu=alcanceRm/abs(zl1u*(L3+(1-m)*L2));
alcanceSppu=alcanceSp/abs(zl1u*(L1+p*L2));
alcanceYppu=alcanceYp/abs(zl1u*(L5+(1-p)*L2));

LimXm=(0.8*(L4+L2+L3))/(L4+m*L2);
LimRm=(0.8*(L4+L2+L3))/(L3+(1-m)*L2);
LimSp=(0.8*(L1+L2+L5))/(L1+m*L2);
LimYp=(0.8*(L1+L2+L5))/(L5+(1-m)*L2);

vectorLimXm(1,i)=LimXm;
vectorLimRm(1,i)=LimRm;
vectorLimSp(1,i)=LimSp;
vectorLimYp(1,i)=LimYp;

Abscisa(fila,i)=m;
OrdenadaXm(fila,i)=alcanceXmpu;
OrdenadaRm(fila,i)=alcanceRmpu;
OrdenadaSp(fila,i)=alcanceSppu;
OrdenadaYp(fila,i)=alcanceYppu;

cadenaym1=sprintf('alcanceXmpu');
cadenaym2=sprintf('alcanceRmpu');
cadenayp1=sprintf('alcanceSppu');
cadenayp2=sprintf('alcanceYppu');

cadenaym1a=sprintf('alcanceXmpua');
cadenaym2a=sprintf('alcanceRmpua');
cadenayp1a=sprintf('alcanceSppua');
cadenayp2a=sprintf('alcanceYppua');

alcanceXmpua=alcanceXm/abs(zl1u*(L4+L2+L3));
alcanceRmpua=alcanceRm/abs(zl1u*(L3+L2+L4));
alcanceSppua=alcanceSp/abs(zl1u*(L1+L2+L5));
alcanceYppua=alcanceYp/abs(zl1u*(L5+L2+L1));

OrdenadaXma(fila,i)=alcanceXmpua;
OrdenadaRma(fila,i)=alcanceRmpua;
OrdenadaSpa(fila,i)=alcanceSppua;
OrdenadaYpa(fila,i)=alcanceYppua;



case{3}, % tramo3
   
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);

% PROTECCION DE DISTANCIA (21)
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

if faltam==1
   cfaltam='ag';
   elementom='AG';
   alcanceXm=abs(Xm_ZAG)/cos(angle(Xm_ZAG)-angle(zl1u));
   alcanceRm=abs(Rm_ZAG)/cos(angle(Rm_ZAG)-angle(zl1u));
   
elseif faltam==2
   cfaltam='bg';
   elementom='BG';
   alcanceXm=abs(Xm_ZBG)/cos(angle(Xm_ZBG)-angle(zl1u));
   alcanceRm=abs(Rm_ZBG)/cos(angle(Rm_ZBG)-angle(zl1u));
   
elseif faltam==3
   cfaltam='cg';
   elementom='CG';
   alcanceXm=abs(Xm_ZCG)/cos(angle(Xm_ZCG)-angle(zl1u));
   alcanceRm=abs(Rm_ZCG)/cos(angle(Rm_ZCG)-angle(zl1u));


elseif faltam==4
   cfaltam='ab';
   elementom='AB';
   alcanceXm=abs(Xm_ZAB)/cos(angle(Xm_ZAB)-angle(zl1u));
   alcanceRm=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));

elseif faltam==5
   cfaltam='bc';
   elementom='BC';
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));
   alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));

elseif faltam==6
   cfaltam='ca';
   elementom='CA';
   alcanceXm=abs(Xm_ZCA)/cos(angle(Xm_ZCA)-angle(zl1u));
   alcanceRm=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));

elseif faltam==7
   cfaltam='abg';
   elementom='AB';
   alcanceXm=abs(Xm_ZAB)/cos(angle(Xm_ZAB)-angle(zl1u));
   alcanceRm=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));

elseif faltam==8

   cfaltam='bcg';
   elementom='BC';
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));
   alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
     
elseif faltam==9
   cfaltam='cag';
   elementom='CA';
   alcanceXm=abs(Xm_ZCA)/cos(angle(Xm_ZCA)-angle(zl1u));
   alcanceRm=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
   
elseif faltam==10
   cfaltam='abc';
   elementom='BC';
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));
   alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
   
elseif faltam==11
   cfaltam='abcg';
   elementom='BC';
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));
   alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
   
end

% Recuerdo: tramo3
   
if faltap==1
   cfaltap='ag';
   elementop='AG';
   alcanceYp=abs(Yp_ZAG)/cos(angle(Yp_ZAG)-angle(zl1u));
   alcanceRp=abs(Rp_ZAG)/cos(angle(Rp_ZAG)-angle(zl1u));

elseif faltap==2
   cfaltap='bg';
   elementop='BG';
   alcanceYp=abs(Yp_ZBG)/cos(angle(Yp_ZBG)-angle(zl1u));
   alcanceRp=abs(Rp_ZBG)/cos(angle(Rp_ZBG)-angle(zl1u));
   
elseif faltap==3
   cfaltap='cg';

   elementop='CG';
   alcanceYp=abs(Yp_ZCG)/cos(angle(Yp_ZCG)-angle(zl1u));
   alcanceRp=abs(Rp_ZCG)/cos(angle(Rp_ZCG)-angle(zl1u));
   
elseif faltap==4
   cfaltap='ab';
   elementop='AB';

   alcanceYp=abs(Yp_ZAB)/cos(angle(Yp_ZAB)-angle(zl1u));
   alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
   
elseif faltap==5
   cfaltap='bc';
   elementop='BC';
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==6
   cfaltap='ca';
   elementop='CA';
   alcanceYp=abs(Yp_ZCA)/cos(angle(Yp_ZCA)-angle(zl1u));
   alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
   
elseif faltap==7
   cfaltap='abg';
   elementop='AB';
alcanceYp=abs(Yp_ZAB)/cos(angle(Yp_ZAB)-angle(zl1u));
alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));

elseif faltap==8
   cfaltap='bcg';
   elementop='BC';
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));

   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==9
   cfaltap='cag';
   elementop='CA';
   alcanceYp=abs(Yp_ZCA)/cos(angle(Yp_ZCA)-angle(zl1u));
   alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
   
elseif faltap==10
   cfaltap='abc';
   elementop='BC';
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==11
   cfaltap='abcg';
   elementop='BC';
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
end

alcanceXmpu=alcanceXm/abs(zl1u*(L4+L2+m*L3));
alcanceRmpu=alcanceRm/abs(zl1u*(1-m)*L3);
alcanceYppu=alcanceYp/abs(zl1u*(L5+p*L3));
alcanceRppu=alcanceRp/abs(zl1u*(1-p)*L3);

LimXm=(0.8*(L4+L2+L3))/(L4+L2+m*L3);
LimRm=(0.8*(L4+L2+L3))/((1-m)*L3);
LimRp=(0.8*(L3+L5))/((1-m)*L3);
LimYp=(0.8*(L3+L5))/(L5+m*L3);

vectorLimXm(1,i)=LimXm;
vectorLimRm(1,i)=LimRm;
vectorLimRp(1,i)=LimRp;
vectorLimYp(1,i)=LimYp;

Abscisa(fila,i)=m;
OrdenadaXm(fila,i)=alcanceXmpu;
OrdenadaRm(fila,i)=alcanceRmpu;
OrdenadaYp(fila,i)=alcanceYppu;
OrdenadaRp(fila,i)=alcanceRppu;

cadenaym1=sprintf('alcanceXmpu');
cadenaym2=sprintf('alcanceRmpu');
cadenayp1=sprintf('alcanceYppu');
cadenayp2=sprintf('alcanceRppu');

cadenaym1a=sprintf('alcanceXmpua');
cadenaym2a=sprintf('alcanceRmpua');
cadenayp1a=sprintf('alcanceYppua');
cadenayp2a=sprintf('alcanceRppua');

alcanceXmpua=alcanceXm/abs(zl1u*(L4+L2+L3));
alcanceRmpua=alcanceRm/abs(zl1u*(L4+L2+L3));
alcanceYppua=alcanceYp/abs(zl1u*(L5+L3));
alcanceRppua=alcanceRp/abs(zl1u*(L5+L3));

OrdenadaXma(fila,i)=alcanceXmpua;
OrdenadaRma(fila,i)=alcanceRmpua;
OrdenadaYpa(fila,i)=alcanceYppua;
OrdenadaRpa(fila,i)=alcanceRppua;

case{4}, %tramo4
   
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);

% PROTECCION DE DISTANCIA (21)
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);
[Xp_ZAB,Xp_ZBC,Xp_ZCA,Xp_ZAG,Xp_ZBG,Xp_ZCG]=...
   PDistComp(Vxa,Vxb,Vxc,Ixpa,Ixpb,Ixpc,Ixp0,k,Km,Ixm0);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
  PDistComp(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k,Km,Iry0);
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);
[Xp_ZAB,Xp_ZBC,Xp_ZCA,Xp_ZAG,Xp_ZBG,Xp_ZCG]=...
   PDist(Vxa,Vxb,Vxc,Ixpa,Ixpb,Ixpc,Ixp0,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
[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

if faltam==1
   cfaltam='ag';
   elementom='AG';
   alcanceSm=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u));
   alcanceXm=abs(Xm_ZAG)/cos(angle(Xm_ZAG)-angle(zl1u));

elseif faltam==2
   cfaltam='bg';
   elementom='BG';
   alcanceSm=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u));
   alcanceXm=abs(Xm_ZBG)/cos(angle(Xm_ZBG)-angle(zl1u));
   
elseif faltam==3
   cfaltam='cg';
   elementom='CG';
   alcanceSm=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u));
   alcanceXm=abs(Xm_ZCG)/cos(angle(Xm_ZCG)-angle(zl1u));

elseif faltam==4
   cfaltam='ab';
   elementom='AB';
   alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
   alcanceXm=abs(Xm_ZAB)/cos(angle(Xm_ZAB)-angle(zl1u));

elseif faltam==5
   cfaltam='bc';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));

elseif faltam==6
   cfaltam='ca';
   elementom='CA';
   alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
   alcanceXm=abs(Xm_ZCA)/cos(angle(Xm_ZCA)-angle(zl1u));

elseif faltam==7
   cfaltam='abg';
   elementom='AB';
   alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
   alcanceXm=abs(Xm_ZAB)/cos(angle(Xm_ZAB)-angle(zl1u));

elseif faltam==8
   cfaltam='bcg';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));

elseif faltam==9
   cfaltam='cag';
   elementom='CA';
   alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
   alcanceXm=abs(Xm_ZCA)/cos(angle(Xm_ZCA)-angle(zl1u));

elseif faltam==10
   cfaltam='abc';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));

elseif faltam==11
   cfaltam='abcg';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceXm=abs(Xm_ZBC)/cos(angle(Xm_ZBC)-angle(zl1u));

end

% Recuerdo: tramo4
   
if faltap==1
   cfaltap='ag';
   elementop='AG';
   alcanceXp=abs(Xp_ZAG)/cos(angle(Xp_ZAG)-angle(zl1u));
   alcanceRp=abs(Rp_ZAG)/cos(angle(Rp_ZAG)-angle(zl1u));

elseif faltap==2
   cfaltap='bg';
   elementop='BG';
   alcanceXp=abs(Xp_ZBG)/cos(angle(Xp_ZBG)-angle(zl1u));
   alcanceRp=abs(Rp_ZBG)/cos(angle(Rp_ZBG)-angle(zl1u));
   
elseif faltap==3
   cfaltap='cg';
   elementop='CG';
   alcanceXp=abs(Xp_ZCG)/cos(angle(Xp_ZCG)-angle(zl1u));
   alcanceRp=abs(Rp_ZCG)/cos(angle(Rp_ZCG)-angle(zl1u));
   
elseif faltap==4
   cfaltap='ab';
   elementop='AB';
   alcanceXp=abs(Xp_ZAB)/cos(angle(Xp_ZAB)-angle(zl1u));
   alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
   
elseif faltap==5
   cfaltap='bc';
   elementop='BC';
   alcanceXp=abs(Xp_ZBC)/cos(angle(Xp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==6
   cfaltap='ca';
   elementop='CA';
   alcanceXp=abs(Xp_ZCA)/cos(angle(Xp_ZCA)-angle(zl1u));
   alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
   
elseif faltap==7
   cfaltap='abg';
   elementop='AB';
alcanceXp=abs(Xp_ZAB)/cos(angle(Xp_ZAB)-angle(zl1u));
alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));

elseif faltap==8
   cfaltap='bcg';
   elementop='BC';
   alcanceXp=abs(Xp_ZBC)/cos(angle(Xp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==9
   cfaltap='cag';
   elementop='CA';
   alcanceXp=abs(Xp_ZCA)/cos(angle(Xp_ZCA)-angle(zl1u));
   alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
   
elseif faltap==10
   cfaltap='abc';
   elementop='BC';
   alcanceXp=abs(Xp_ZBC)/cos(angle(Xp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==11
   cfaltap='abcg';
   elementop='BC';
   alcanceXp=abs(Xp_ZBC)/cos(angle(Xp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
end

alcanceSmpu=alcanceSm/abs(zl1u*(L1+m*L4));
alcanceXmpu=alcanceXm/abs(zl1u*(1-m)*L4);
alcanceXppu=alcanceXp/abs(zl1u*(1-p)*L4);
alcanceRppu=alcanceRp/abs(zl1u*(L3+L2+p*L4));

LimSm=(0.8*(L1+L4))/(L1+m*L4);
LimXm=(0.8*(L1+L4))/((1-m)*L4);
LimXp=(0.8*(L4+L2+L3))/((1-m)*L4);
LimRp=(0.8*(L4+L2+L3))/(L3+L2+m*L4);

vectorLimSm(1,i)=LimSm;
vectorLimXm(1,i)=LimXm;
vectorLimXp(1,i)=LimXp;
vectorLimRp(1,i)=LimRp;

Abscisa(fila,i)=m;
OrdenadaSm(fila,i)=alcanceSmpu;
OrdenadaXm(fila,i)=alcanceXmpu;
OrdenadaXp(fila,i)=alcanceXppu;
OrdenadaRp(fila,i)=alcanceRppu;

cadenaym1=sprintf('alcanceSmpu');
cadenaym2=sprintf('alcanceXmpu');
cadenayp1=sprintf('alcanceXppu');
cadenayp2=sprintf('alcanceRppu');

cadenaym1a=sprintf('alcanceSmpua');
cadenaym2a=sprintf('alcanceXmpua');
cadenayp1a=sprintf('alcanceXppua');
cadenayp2a=sprintf('alcanceRppua');

alcanceSmpua=alcanceSm/abs(zl1u*(L1+L4));
alcanceXmpua=alcanceXm/abs(zl1u*(L1+L4));
alcanceXppua=alcanceXp/abs(zl1u*(L4+L2+L3));
alcanceRppua=alcanceRp/abs(zl1u*(L4+L2+L3));

OrdenadaSma(fila,i)=alcanceSmpua;
OrdenadaXma(fila,i)=alcanceXmpua;
OrdenadaXpa(fila,i)=alcanceXppua;
OrdenadaRpa(fila,i)=alcanceRppua;

case{5}, %tramo5
   
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);

% PROTECCION DE DISTANCIA (21)
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

if faltam==1
   cfaltam='ag';
   elementom='AG';
   alcanceSm=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u));
   alcanceYm=abs(Ym_ZAG)/cos(angle(Ym_ZAG)-angle(zl1u));

elseif faltam==2
   cfaltam='bg';
   elementom='BG';
   alcanceSm=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u));
   alcanceYm=abs(Ym_ZBG)/cos(angle(Ym_ZBG)-angle(zl1u));
   
elseif faltam==3
   cfaltam='cg';
   elementom='CG';
   alcanceSm=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u));
   alcanceYm=abs(Ym_ZCG)/cos(angle(Ym_ZCG)-angle(zl1u));

elseif faltam==4
   cfaltam='ab';
   elementom='AB';
   alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
   alcanceYm=abs(Ym_ZAB)/cos(angle(Ym_ZAB)-angle(zl1u));

elseif faltam==5
   cfaltam='bc';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceYm=abs(Ym_ZBC)/cos(angle(Ym_ZBC)-angle(zl1u));

elseif faltam==6
   cfaltam='ca';
   elementom='CA';
   alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
   alcanceYm=abs(Ym_ZCA)/cos(angle(Ym_ZCA)-angle(zl1u));

elseif faltam==7
   cfaltam='abg';
   elementom='AB';
   alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
   alcanceYm=abs(Ym_ZAB)/cos(angle(Ym_ZAB)-angle(zl1u));

elseif faltam==8
   cfaltam='bcg';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceYm=abs(Ym_ZBC)/cos(angle(Ym_ZBC)-angle(zl1u));

elseif faltam==9
   cfaltam='cag';
   elementom='CA';
   alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
   alcanceYm=abs(Ym_ZCA)/cos(angle(Ym_ZCA)-angle(zl1u));

elseif faltam==10
   cfaltam='abc';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceYm=abs(Ym_ZBC)/cos(angle(Ym_ZBC)-angle(zl1u));

elseif faltam==11
   cfaltam='abcg';
   elementom='BC';
   alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
   alcanceYm=abs(Ym_ZBC)/cos(angle(Ym_ZBC)-angle(zl1u));

end

% Recuerdo: tramo5
   
if faltap==1
   cfaltap='ag';
   elementop='AG';
   alcanceYp=abs(Yp_ZAG)/cos(angle(Yp_ZAG)-angle(zl1u));
   alcanceRp=abs(Rp_ZAG)/cos(angle(Rp_ZAG)-angle(zl1u));

elseif faltap==2
   cfaltap='bg';
   elementop='BG';
   alcanceYp=abs(Yp_ZBG)/cos(angle(Yp_ZBG)-angle(zl1u));
   alcanceRp=abs(Rp_ZBG)/cos(angle(Rp_ZBG)-angle(zl1u));
   
elseif faltap==3
   cfaltap='cg';
   elementop='CG';
   alcanceYp=abs(Yp_ZCG)/cos(angle(Yp_ZCG)-angle(zl1u));
   alcanceRp=abs(Rp_ZCG)/cos(angle(Rp_ZCG)-angle(zl1u));
   
elseif faltap==4
   cfaltap='ab';
   elementop='AB';
   alcanceYp=abs(Yp_ZAB)/cos(angle(Yp_ZAB)-angle(zl1u));
   alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
   
elseif faltap==5
   cfaltap='bc';
   elementop='BC';
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==6
   cfaltap='ca';

   elementop='CA';
   alcanceYp=abs(Yp_ZCA)/cos(angle(Yp_ZCA)-angle(zl1u));
   alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
   
elseif faltap==7
   cfaltap='abg';
   elementop='AB';
alcanceYp=abs(Yp_ZAB)/cos(angle(Yp_ZAB)-angle(zl1u));
alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));

elseif faltap==8
   cfaltap='bcg';
   elementop='BC';
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==9
   cfaltap='cag';
   elementop='CA';
   alcanceYp=abs(Yp_ZCA)/cos(angle(Yp_ZCA)-angle(zl1u));
   alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
   
elseif faltap==10
   cfaltap='abc';
   elementop='BC';
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
elseif faltap==11
   cfaltap='abcg';
   elementop='BC';
   alcanceYp=abs(Yp_ZBC)/cos(angle(Yp_ZBC)-angle(zl1u));
   alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
   
end

alcanceSmpu=alcanceSm/abs(zl1u*(L1+L2+m*L5));
alcanceYmpu=alcanceYm/abs(zl1u*(1-m)*L5);
alcanceYppu=alcanceYp/abs(zl1u*(1-p)*L5);
alcanceRppu=alcanceRp/abs(zl1u*(L3+p*L5));

LimSm=(0.8*(L1+L2+L5))/(L1+L2+m*L5);
LimYm=(0.8*(L1+L2+L5))/((1-m)*L5);
LimYp=(0.8*(L5+L3))/((1-m)*L5);
LimRp=(0.8*(L5+L3))/(L3+m*L5);

vectorLimSm(1,i)=LimSm;
vectorLimYm(1,i)=LimYm;
vectorLimYp(1,i)=LimYp;
vectorLimRp(1,i)=LimRp;

Abscisa(fila,i)=m;
OrdenadaSm(fila,i)=alcanceSmpu;
OrdenadaYm(fila,i)=alcanceYmpu;
OrdenadaYp(fila,i)=alcanceYppu;
OrdenadaRp(fila,i)=alcanceRppu;

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


cadenaym1a=sprintf('alcanceSmpua');
cadenaym2a=sprintf('alcanceYmpua');
cadenayp1a=sprintf('alcanceYppua');
cadenayp2a=sprintf('alcanceRppua');


alcanceSmpua=alcanceSm/abs(zl1u*(L1+L2+L5));
alcanceYmpua=alcanceYm/abs(zl1u*(L1+L2+L5));
alcanceYppua=alcanceYp/abs(zl1u*(L3+L5));
alcanceRppua=alcanceRp/abs(zl1u*(L3+L5));

OrdenadaSma(fila,i)=alcanceSmpua;
OrdenadaYma(fila,i)=alcanceYmpua;
OrdenadaYpa(fila,i)=alcanceYppua;
OrdenadaRpa(fila,i)=alcanceRppua;


end % end del switch(tramo)

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

fila=fila+1;
k0=k0*10;
end  % end del while k0

if compens21==1
   comp='21c';
elseif compens21==2
   comp='21nc';
end

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

Abscisa1m=Abscisa(1:1,1:(i-1));
Abscisa2m=Abscisa(2:2,1:(i-1));
Abscisa3m=Abscisa(3:3,1:(i-1));

Abscisa4m=Abscisa(1:1,1:(i-1));
Abscisa5m=Abscisa(2:2,1:(i-1));
Abscisa6m=Abscisa(3:3,1:(i-1));

Abscisa7m=Abscisa(1:1,1:(i-1));
Abscisa8m=Abscisa(2:2,1:(i-1));
Abscisa9m=Abscisa(3:3,1:(i-1));

Abscisa10m=Abscisa(1:1,1:(i-1));
Abscisa11m=Abscisa(2:2,1:(i-1));
Abscisa12m=Abscisa(3:3,1:(i-1));

Abscisa1ma=Abscisa(1:1,1:(i-1));
Abscisa2ma=Abscisa(2:2,1:(i-1));
Abscisa3ma=Abscisa(3:3,1:(i-1));

Abscisa4ma=Abscisa(1:1,1:(i-1));
Abscisa5ma=Abscisa(2:2,1:(i-1));
Abscisa6ma=Abscisa(3:3,1:(i-1));

Abscisa7ma=Abscisa(1:1,1:(i-1));
Abscisa8ma=Abscisa(2:2,1:(i-1));
Abscisa9ma=Abscisa(3:3,1:(i-1));

Abscisa10ma=Abscisa(1:1,1:(i-1));
Abscisa11ma=Abscisa(2:2,1:(i-1));
Abscisa12ma=Abscisa(3:3,1:(i-1));

variassalidas=1;
while variassalidas==1

switch(tramo)
case{1}, % tramo 1
   
fprintf('Salida Grafica 1: Extr.m(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Extr.m(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
   
   salida=input('Salida: 1,2 :');
   
  
switch(salida)
       case{1}, % salida 1, para el tramo 1

Ordenada1m=OrdenadaSm(1:1,1:(i-1));
Ordenada2m=OrdenadaSm(2:2,1:(i-1));
Ordenada3m=OrdenadaSm(3:3,1:(i-1));

Ordenada4m=OrdenadaXm(1:1,1:(i-1));
Ordenada5m=OrdenadaXm(2:2,1:(i-1));
Ordenada6m=OrdenadaXm(3:3,1:(i-1));

Ordenada7m=OrdenadaSp(1:1,1:(i-1));
Ordenada8m=OrdenadaSp(2:2,1:(i-1));
Ordenada9m=OrdenadaSp(3:3,1:(i-1));

Ordenada10m=OrdenadaYp(1:1,1:(i-1));
Ordenada11m=OrdenadaYp(2:2,1:(i-1));
Ordenada12m=OrdenadaYp(3:3,1:(i-1));

Ordenadami=vectorLimSm(1:1,1:(i-1));
Ordenadamd=vectorLimXm(1:1,1:(i-1));
Ordenadapi=vectorLimSp(1:1,1:(i-1));
Ordenadapd=vectorLimYp(1:1,1:(i-1));


       case{2}, % salida 2, para el tramo 1

Ordenada1ma=OrdenadaSma(1:1,1:(i-1));
Ordenada2ma=OrdenadaSma(2:2,1:(i-1));
Ordenada3ma=OrdenadaSma(3:3,1:(i-1));

Ordenada4ma=OrdenadaXma(1:1,1:(i-1));
Ordenada5ma=OrdenadaXma(2:2,1:(i-1));
Ordenada6ma=OrdenadaXma(3:3,1:(i-1));

Ordenada7ma=OrdenadaSpa(1:1,1:(i-1));
Ordenada8ma=OrdenadaSpa(2:2,1:(i-1));
Ordenada9ma=OrdenadaSpa(3:3,1:(i-1));

Ordenada10ma=OrdenadaYpa(1:1,1:(i-1));
Ordenada11ma=OrdenadaYpa(2:2,1:(i-1));
Ordenada12ma=OrdenadaYpa(3:3,1:(i-1));

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

mm=0.01;
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.01;
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.01;
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 el tramo 1 

case{2}, % tramo 2
   
fprintf('Salida Grafica 1: Extr.m(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Extr.m(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
   
   salida=input('Salida: 1,2 :');

switch(salida)
       case{1}, % salida 1, para el tramo 2
   
Ordenada1m=OrdenadaXm(1:1,1:(i-1));
Ordenada2m=OrdenadaXm(2:2,1:(i-1));
Ordenada3m=OrdenadaXm(3:3,1:(i-1));

Ordenada4m=OrdenadaRm(1:1,1:(i-1));
Ordenada5m=OrdenadaRm(2:2,1:(i-1));
Ordenada6m=OrdenadaRm(3:3,1:(i-1));

Ordenada7m=OrdenadaSp(1:1,1:(i-1));
Ordenada8m=OrdenadaSp(2:2,1:(i-1));
Ordenada9m=OrdenadaSp(3:3,1:(i-1));

Ordenada10m=OrdenadaYp(1:1,1:(i-1));
Ordenada11m=OrdenadaYp(2:2,1:(i-1));
Ordenada12m=OrdenadaYp(3:3,1:(i-1));

Ordenadami=vectorLimXm(1:1,1:(i-1));
Ordenadamd=vectorLimRm(1:1,1:(i-1));
Ordenadapi=vectorLimSp(1:1,1:(i-1));
Ordenadapd=vectorLimYp(1:1,1:(i-1));

       case{2}, % salida 2, para el tramo 2

Ordenada1ma=OrdenadaXma(1:1,1:(i-1));
Ordenada2ma=OrdenadaXma(2:2,1:(i-1));
Ordenada3ma=OrdenadaXma(3:3,1:(i-1));

Ordenada4ma=OrdenadaRma(1:1,1:(i-1));
Ordenada5ma=OrdenadaRma(2:2,1:(i-1));
Ordenada6ma=OrdenadaRma(3:3,1:(i-1));

Ordenada7ma=OrdenadaSpa(1:1,1:(i-1));
Ordenada8ma=OrdenadaSpa(2:2,1:(i-1));
Ordenada9ma=OrdenadaSpa(3:3,1:(i-1));

Ordenada10ma=OrdenadaYpa(1:1,1:(i-1));
Ordenada11ma=OrdenadaYpa(2:2,1:(i-1));
Ordenada12ma=OrdenadaYpa(3:3,1:(i-1));

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

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

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

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

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

case{3}, %tramo 3
   
fprintf('Salida Grafica 1: Extr.m(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Extr.m(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
   
      salida=input('Salida: 1,2 :');

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

Ordenada1m=OrdenadaXm(1:1,1:(i-1));
Ordenada2m=OrdenadaXm(2:2,1:(i-1));
Ordenada3m=OrdenadaXm(3:3,1:(i-1));

Ordenada4m=OrdenadaRm(1:1,1:(i-1));
Ordenada5m=OrdenadaRm(2:2,1:(i-1));
Ordenada6m=OrdenadaRm(3:3,1:(i-1));

Ordenada7m=OrdenadaYp(1:1,1:(i-1));
Ordenada8m=OrdenadaYp(2:2,1:(i-1));
Ordenada9m=OrdenadaYp(3:3,1:(i-1));

Ordenada10m=OrdenadaRp(1:1,1:(i-1));
Ordenada11m=OrdenadaRp(2:2,1:(i-1));
Ordenada12m=OrdenadaRp(3:3,1:(i-1));

Ordenadami=vectorLimXm(1:1,1:(i-1));
Ordenadamd=vectorLimRm(1:1,1:(i-1));
Ordenadapi=vectorLimYp(1:1,1:(i-1));
Ordenadapd=vectorLimRp(1:1,1:(i-1));

case{2}, % salida 2, para el tramo 3
   
Ordenada1ma=OrdenadaXma(1:1,1:(i-1));
Ordenada2ma=OrdenadaXma(2:2,1:(i-1));
Ordenada3ma=OrdenadaXma(3:3,1:(i-1));

Ordenada4ma=OrdenadaRma(1:1,1:(i-1));
Ordenada5ma=OrdenadaRma(2:2,1:(i-1));
Ordenada6ma=OrdenadaRma(3:3,1:(i-1));

Ordenada7ma=OrdenadaYpa(1:1,1:(i-1));
Ordenada8ma=OrdenadaYpa(2:2,1:(i-1));
Ordenada9ma=OrdenadaYpa(3:3,1:(i-1));

Ordenada10ma=OrdenadaRpa(1:1,1:(i-1));
Ordenada11ma=OrdenadaRpa(2:2,1:(i-1));
Ordenada12ma=OrdenadaRpa(3:3,1:(i-1));

mm=0.01;
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.01;
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.01;
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.01;
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}, % tramo 4
   
fprintf('Salida Grafica 1: Extr.m(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Extr.m(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
   
      salida=input('Salida: 1,2 :');

switch(salida)
         
case{1}, % salida 1, para el tramo 4
   
Ordenada1m=OrdenadaSm(1:1,1:(i-1));
Ordenada2m=OrdenadaSm(2:2,1:(i-1));
Ordenada3m=OrdenadaSm(3:3,1:(i-1));

Ordenada4m=OrdenadaXm(1:1,1:(i-1));
Ordenada5m=OrdenadaXm(2:2,1:(i-1));
Ordenada6m=OrdenadaXm(3:3,1:(i-1));

Ordenada7m=OrdenadaXp(1:1,1:(i-1));

Ordenada8m=OrdenadaXp(2:2,1:(i-1));
Ordenada9m=OrdenadaXp(3:3,1:(i-1));

Ordenada10m=OrdenadaRp(1:1,1:(i-1));
Ordenada11m=OrdenadaRp(2:2,1:(i-1));
Ordenada12m=OrdenadaRp(3:3,1:(i-1));

Ordenadami=vectorLimSm(1:1,1:(i-1));
Ordenadamd=vectorLimXm(1:1,1:(i-1));
Ordenadapi=vectorLimXp(1:1,1:(i-1));
Ordenadapd=vectorLimRp(1:1,1:(i-1));

case{2}, % salida 2, para el tramo 4
   
Ordenada1ma=OrdenadaSma(1:1,1:(i-1));
Ordenada2ma=OrdenadaSma(2:2,1:(i-1));
Ordenada3ma=OrdenadaSma(3:3,1:(i-1));

Ordenada4ma=OrdenadaXma(1:1,1:(i-1));
Ordenada5ma=OrdenadaXma(2:2,1:(i-1));
Ordenada6ma=OrdenadaXma(3:3,1:(i-1));

Ordenada7ma=OrdenadaXpa(1:1,1:(i-1));
Ordenada8ma=OrdenadaXpa(2:2,1:(i-1));
Ordenada9ma=OrdenadaXpa(3:3,1:(i-1));

Ordenada10ma=OrdenadaRpa(1:1,1:(i-1));
Ordenada11ma=OrdenadaRpa(2:2,1:(i-1));
Ordenada12ma=OrdenadaRpa(3:3,1:(i-1));

mm=0.01;
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.01;
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.01;
cont=1;
while mm<1;
Abscisap0(cont)=mm;
Ordenadap0(cont)=L4/(L4+L2+L3)-mm*(L4/(L4+L2+L3));
cont=cont+1;
mm=mm+0.01;
end

mm=0.01;
cont=1;
while mm<1;
Abscisap00(cont)=mm;
Ordenadap00(cont)=(L3+L2)/(L3+L2+L4)+mm*(L4/(L3+L2+L4));
cont=cont+1;
mm=mm+0.01;
end

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

case{5},
   
fprintf('Salida Grafica 1: Extr.m(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Extr.m(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
fprintf('                : Extr.p(2Subpl), Unidad Pref.(pua),k0=0.1-1-10\n');
   
    salida=input('Salida: 1,2 :');


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

Ordenada1m=OrdenadaSm(1:1,1:(i-1));
Ordenada2m=OrdenadaSm(2:2,1:(i-1));
Ordenada3m=OrdenadaSm(3:3,1:(i-1));

Ordenada4m=OrdenadaYm(1:1,1:(i-1));
Ordenada5m=OrdenadaYm(2:2,1:(i-1));
Ordenada6m=OrdenadaYm(3:3,1:(i-1));

Ordenada7m=OrdenadaYp(1:1,1:(i-1));
Ordenada8m=OrdenadaYp(2:2,1:(i-1));
Ordenada9m=OrdenadaYp(3:3,1:(i-1));

Ordenada10m=OrdenadaRp(1:1,1:(i-1));
Ordenada11m=OrdenadaRp(2:2,1:(i-1));
Ordenada12m=OrdenadaRp(3:3,1:(i-1));

Ordenadami=vectorLimSm(1:1,1:(i-1));
Ordenadamd=vectorLimYm(1:1,1:(i-1));
Ordenadapi=vectorLimYp(1:1,1:(i-1));
Ordenadapd=vectorLimRp(1:1,1:(i-1));

case{2}, % salida 2, para el tramo 5
   
Ordenada1ma=OrdenadaSma(1:1,1:(i-1));
Ordenada2ma=OrdenadaSma(2:2,1:(i-1));
Ordenada3ma=OrdenadaSma(3:3,1:(i-1));


Ordenada4ma=OrdenadaYma(1:1,1:(i-1));
Ordenada5ma=OrdenadaYma(2:2,1:(i-1));
Ordenada6ma=OrdenadaYma(3:3,1:(i-1));

Ordenada7ma=OrdenadaYpa(1:1,1:(i-1));
Ordenada8ma=OrdenadaYpa(2:2,1:(i-1));
Ordenada9ma=OrdenadaYpa(3:3,1:(i-1));

Ordenada10ma=OrdenadaRpa(1:1,1:(i-1));
Ordenada11ma=OrdenadaRpa(2:2,1:(i-1));
Ordenada12ma=OrdenadaRpa(3:3,1:(i-1));

mm=0.01;
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.01;
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.01;
cont=1;
while mm<1;
Abscisap0(cont)=mm;
Ordenadap0(cont)=L5/(L5+L3)-mm*(L5/(L5+L3));
cont=cont+1;
mm=mm+0.01;
end

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

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

end % end del switch(tramo)
% -----------------------------------------------------------------------
cadenam1=sprintf('PG54ni,%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\nSccr=%-5.0f,Sccy=%-5.0f',...
   Sccs,Sccx,Sccr,Sccy);
cadenam3=sprintf('Kx0=%-6.2f,Kr0=%-6.2f,Ky0=%-6.2f',...
   Kx0,Kr0,Ky0);
cadenam4=sprintf('m: %s, %s, Rfm=%-3.0f',cfaltam,elementom,Rfm);
if faltam==7|faltam==8|faltam==9|faltam==11
cadenam5=sprintf('Rgm=%-3.0f',Rgm);
cadena4=sprintf('%s,%s',cadenam4,cadenam5);
end
cadenam=sprintf('%s,%s,%s,%s',cadenam1,cadenam2,cadenam3,cadenam4);

cadenap=sprintf('p: %s, %s, Rfp=%-3.0f',cfaltap,elementop,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);
% -----------------------------------------------------------------------

switch(salida)
   
case{1}, % Salida 1
   
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
cont=cont+1;
mm=mm+0.01;
end
   
AnalisisLinea=input('Extremos m=1, p=2: ');
switch(AnalisisLinea)
   
case{1}, % Salida 1, Extremos m  
   
ejemasy1=input('Subplot 1: Eje +y: ');
ejemenosy1=input('Subplot 1: Eje -y (con signo): ');

ejemasy2=input('Subplot 2: Eje +y: ');
ejemenosy2=input('Subplot 2: Eje -y (con signo): ');
   
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa1m,Ordenadami,'m',Abscisa1m,Ordenada1m,'r',Abscisa2m,Ordenada2m,'g',...
   Abscisa3m,Ordenada3m,'b'),grid on
axis([0 1 ejemenosy1 ejemasy1]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenaym1)
%legend('Ref','Lim','k0=0.1','k0=1','k0=10')

subplot(2,1,2) 
plot(Abscisa0,Ordenada0,'c',Abscisa4m,Ordenadamd,'m',Abscisa4m,Ordenada4m,'r',Abscisa5m,Ordenada5m,'g',...
   Abscisa6m,Ordenada6m,'b'),grid on
axis([0 1 ejemenosy2 ejemasy2]);
%title(cadenamp)
xlabel('m=p'),ylabel(cadenaym2)
legend('Ref','Lim','k0=0.1','k0=1','k0=10')

case{2}, % Salida 1, Extremos p
   
ejemasy1=input('Subplot 1: Eje +y: ');
ejemenosy1=input('Subplot 1: Eje -y (con signo): ');

ejemasy2=input('Subplot 2: Eje +y: ');
ejemenosy2=input('Subplot 2: Eje -y (con signo): ');

subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa7m,Ordenadapi,'m',Abscisa7m,Ordenada7m,'r',Abscisa8m,Ordenada8m,'g',...
   Abscisa9m,Ordenada9m,'b'),grid on
axis([0 1 ejemenosy1 ejemasy1]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenayp1)
%legend('Ref','Lim','k0=0.1','k0=1','k0=10')

subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa10m,Ordenadapd,'m',Abscisa10m,Ordenada10m,'r',Abscisa11m,Ordenada11m,'g',...
   Abscisa12m,Ordenada12m,'b'),grid on
axis([0 1 ejemenosy2 ejemasy2]);
%title(cadenamp)
xlabel('m=p'),ylabel(cadenayp2)
legend('Ref','Lim','k0=0.1','k0=1','k0=10')

end

case{2}, % Salida 2
   
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenadaz1(cont)=0.8;
cont=cont+1;
mm=mm+0.01;
end
   
   
AnalisisLinea=input('Extremos m=1, p=2: ');
switch(AnalisisLinea)
   
case{1}, % salida 2, Extremos m

 subplot(2,1,1)
 plot(Abscisa0,Ordenadaz1,'m',Abscisam0,Ordenadam0,'c',Abscisa1ma,Ordenada1ma,'r',...
    Abscisa2ma,Ordenada2ma,'g',Abscisa3ma,Ordenada3ma,'b'),grid on
%axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenaym1a)
%legend('z1','Ref,Cob','k0=0.1','k0=1','k0=10')

subplot(2,1,2)
plot(Abscisa0,Ordenadaz1,'m',Abscisam00,Ordenadam00,'c',Abscisa4ma,Ordenada4ma,'r',...
   Abscisa5ma,Ordenada5ma,'g',Abscisa6ma,Ordenada6ma,'b'),grid on
%axis([0 1 0 2]);
%title(cadenamp)
xlabel('m=p'),ylabel(cadenaym2a)
legend('z1','Ref,Cob','k0=0.1','k0=1','k0=10')

case{2}, % salida 2, Extremos p

subplot(2,1,1)
plot(Abscisa0,Ordenadaz1,'m',Abscisap0,Ordenadap0,'c',Abscisa7ma,Ordenada7ma,'r',...
   Abscisa8ma,Ordenada8ma,'g',Abscisa9ma,Ordenada9ma,'b'),grid on
%axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel(cadenayp1a)
%legend('z1','Ref,Cob','k0=0.1','k0=1','k0=10')

subplot(2,1,2)
plot(Abscisa0,Ordenadaz1,'m',Abscisap00,Ordenadap00,'c',Abscisa10ma,Ordenada10ma,'r',...
   Abscisa11ma,Ordenada11ma,'g',Abscisa12ma,Ordenada12ma,'b'),grid on
%axis([0 1 0 2]);
%title(cadenamp)
xlabel('m=p'),ylabel(cadenayp2a)
legend('z1','Ref,Cob','k0=0.1','k0=1','k0=10')
end 
end

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%252FPGn%252F">PGn</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPGn%252FPG54ni%252F">PG54ni</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPGn%252FPG54ni%252FFG54ni.m%252F">FG54ni.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>