clc
variasejecuciones=1;
while variasejecuciones==1
clear
% ---------------------------------------------------------------
salida=input('Pantalla=1, Archivo=2 :');
if salida==2
cadena=input('Poner NombreArchivo.txt:','s');
cadena1=strcat('Salida\',cadena);
end
if salida==1
fid=1;
elseif salida==2
fid=fopen(cadena1,'w+');
end
% ---------------------------------------------------------------
fprintf(fid,'PC28Mb\n');
fprintf(fid,'FC28Mb.m\n');
fprintf(fid,'Falta Multicircuito, Lineas a Distintas Tensiones (m=p)\n');
fprintf(fid,'Barra m=3, Barra p=4\n');
% ---------------------------------------------------------------
%PC28Mb
datos=input('datos: standard=1,input=2: ');
if datos==1
UA=220;
LA=66;
LB=LA;
UB=66;
kE=1;
Compens21=2;
Sccs=5000;
Sp=120; Ss=120; St=40; Up=220; Us=66; Ut=33;
tccps=12; tccpt=5; tccst=2;
else
% ---------------------------------------------------------------
UA=input('UA(Fase-Fase) kV: ');
LA=input('Longitud Linea LA km: ');
UB=input('UB(Fase-Fase) kV: ');
kE=input('kE=E/Enominal: ');
Compens21=input('Compensacion 21: Si=1, No=2 :');
Sccs=input('Sccs(MVA): ');
Sp=input('Sp(MVA): ');
Ss=input('Ss(MVA): ');
St=input('St(MVA): ');
tccps=input('tccps(%): ');
tccpt=input('tccpt(%): ');
tccst=input('tccst(%): ');
end
Ks0=input('Ks0=zs0/zs1: ');
m=input('m=p pu : ');
LB=LA;
p=m;
fprintf(fid,'\n');
% ---------------------------------------------------------------
if UA==400
zlA1u=0.0346+0.3208i; zlA2u=zlA1u; zlA0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i;
elseif UA==220
zlA1u=0.0829+0.4135i; zlA2u=zlA1u; zlA0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i;
elseif UA==132
zlA1u=0.1342+0.4042i; zlA2u=zl1u; zlA0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i;
elseif UA==66
zlA1u=0.2156+0.3952i; zlA2u=zl1u; zlA0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i;
end
if UB==400
zlB1u=0.0346+0.3208i; zlB2u=zlB1u; zlB0u=0.1828+1.1482i;
elseif UB==220
zlB1u=0.0829+0.4135i; zlB2u=zlB1u; zlB0u=0.2311+1.3085i;
elseif UB==132
zlB1u=0.1342+0.4042i; zlB2u=zlB1u; zlB0u=0.2824+1.3718i;
elseif UB==66
zlB1u=0.2156+0.3952i; zlB2u=zlB1u; zlB0u=0.3636+1.4278i;
end
% ---------------------------------------------------------------
% CALCULOS
% -----------------------------------------------------------------------
EAv=UA*1000/sqrt(3)+0i;
EA=kE;
EBv=UB*1000/sqrt(3)+0i;
EB=kE;
E=kE;
kA=(zlA0u-zlA1u)/(3*zlA1u);
kB=(zlB0u-zlB1u)/(3*zlB1u);
SAm=zlA1u*m*LA; SARA=zlA1u*LA; % Ohm 220 kV
SBp=zlB1u*p*LB; SBRB=zlB1u*LB; % Ohm 66 kV
Sbase=100; % MVA Trifasicos
UbaseA=UA; % kV Fase-Fase
UbaseB=UB; % kV Fase-Fase
VbaseA=UbaseA/sqrt(3); % kV Fase-Neutro
VbaseB=UbaseB/sqrt(3); % kV Fase-Neutro
zbaseA=UA^2/Sbase; % Ohm
zbaseB=UB^2/Sbase; % Ohm
IbaseA=(UA/sqrt(3))/zbaseA; % kA
IbaseB=(UB/sqrt(3))/zbaseB; % kA
VbaseA=VbaseA*1000; % V
VbaseB=VbaseB*1000; % V
IbaseA=IbaseA*1000; % A
IbaseB=IbaseB*1000; % A
KmA=zlM0u/zlA1u;
KmA=KmA*(IbaseB/IbaseA);
KmB=zlM0u/zlB1u;
KmB=KmB*(IbaseA/IbaseB);
% REDES DE SECUENCIA 0,1,2 pu
[zs1,zSAm1,zSBp1]=pusec1(UA,Sccs,zlA1u,m,LA,zlB1u,p,LB,zbaseA,zbaseB);%pu
[zs2,zSAm2,zSBp2]=pusec2(UA,Sccs,zlA2u,m,LA,zlB2u,p,LB,zbaseA,zbaseB);%pu
[zs0,zSAm0,zSBp0,zMSm0]=...
pusec0(Ks0,zs1,zlA0u,m,LA,zlB0u,p,LB,zlM0u,zbaseA,zbaseB,UbaseA,UbaseB,Sbase);%pu
zpTS=(1/2)*((tccps/100)*UA^2/Sp+(tccpt/100)*UA^2/St-(tccst/100)*UA^2/St);%Ohm 220
zpTS=zpTS/zbaseA;%pu
zsTS=(1/2)*((tccps/100)*UA^2/Sp-(tccpt/100)*UA^2/St+(tccst/100)*UA^2/St);%Ohm 220
zsTS=zsTS/zbaseA;%pu
ztTS=(1/2)*(-(tccps/100)*UA^2/Sp+(tccpt/100)*UA^2/St+(tccst/100)*UA^2/St);%Ohm 220
ztTS=ztTS/zbaseA;%pu
% CALCULO DE LA REDUCCION DE REDES 0,1,2 A LOS PUNTOS m,p DE FALTA
[Minc0]=subMinc0(1);
[Minc1]=subMinc1(1);
[Minc2]=subMinc2(1);
[Minca1]=subMinca1(1);
[z0]=subz0(zs0,zSAm0,zSBp0,zMSm0,zpTS,zsTS,ztTS); %pu
[z1]=subz1(zs1,zSAm1,zSBp1,zpTS,zsTS); %pu
[z2]=subz2(zs2,zSAm2,zSBp2,zpTS,zsTS); %pu
% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc0,z0);%pu
[Y1]=subY1(Minc1,z1);%pu
[Y2]=subY2(Minc2,z2);%pu
% ZBUS(0,1,2)
Z0=inv(Y0);%pu
Z1=inv(Y1);%pu
Z2=inv(Y2);%pu
% 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,falta,Rf,Rg]=...
maMDmp(Z0,Z1,Z2,E,VbaseA,VbaseB,IbaseA,IbaseB,fid); %pu
% IBUS(0,1,2)
I0=[0;0;-Im0;-Ip0];%pu
I1=[0;0;-Im1;-Ip1];%pu
I2=[0;0;-Im2;-Ip2];%pu
% VBUS(0,1,2)
V0=Z0*I0;%pu
V1=Z1*I1+E; V1a=[V1;E];%pu
V2=Z2*I2;%pu
% ASIGNACIONES V,I(0,1,2)
[Ie0]=subIe0(Minc0,V0,z0);%pu
[Ie1]=subIe1(Minca1,V1a,z1);%pu
[Ie2]=subIe2(Minc2,V2,z2);%pu
VSA0=V0(1); VSA1=V1(1); VSA2=V2(1);%pu
VSB0=V0(2); VSB1=V1(2); VSB2=V2(2);%pu
Is0=Ie0(1); Is1=Ie1(1); Is2=Ie2(1); %pu
ISAm0=Ie0(2); ISAm1=Ie1(2); ISAm2=Ie2(2); %pu
ISBp0=Ie0(3); ISBp1=Ie1(3); ISBp2=Ie2(3); %pu
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS PUNTOS m,p DE FALTA
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2); %pu
[Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2); %pu
[Ima,Imb,Imc]=abc(Im0,Im1,Im2); %pu
[Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2); %pu
[Isa,Isb,Isc]=abc(Is0,Is1,Is2); %pu
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS EXTREMOS DE LINEA
[VSAa,VSAb,VSAc]=abc(VSA0,VSA1,VSA2); %pu
[VSBa,VSBb,VSBc]=abc(VSB0,VSB1,VSB2); %pu
[ISAma,ISAmb,ISAmc]=abc(ISAm0,ISAm1,ISAm2); %pu
[ISBpa,ISBpb,ISBpc]=abc(ISBp0,ISBp1,ISBp2); %pu
% PROTECCION DE DISTANCIA (21)
switch(Compens21)
case{1},
% COMPENSADA
[SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG]=...
PDistComp(VSAa,VSAb,VSAc,ISAma,ISAmb,ISAmc,ISAm0,kA,KmA,ISBp0); %pu
[SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG]=...
PDistComp(VSBa,VSBb,VSBc,ISBpa,ISBpb,ISBpc,ISBp0,kB,KmB,ISAm0); %pu
case{2},
% NO COMPENSADA
[SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG]=...
PDist(VSAa,VSAb,VSAc,ISAma,ISAmb,ISAmc,ISAm0,kA);
[SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG]=...
PDist(VSBa,VSBb,VSBc,ISBpa,ISBpb,ISBpc,ISBp0,kB);
end
%
% CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION
% DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA
%
[VSApol,ISAmop,AngDirSAm]=PDire(VSA0,ISAm0); %pu
[VSBpol,ISBpop,AngDirSBp]=PDire(VSB0,ISBp0); %pu
%
if falta==1 cfalta='Aa';elseif falta==2 cfalta='Bc';elseif falta==3 cfalta='Cb';
elseif falta==4 cfalta='Aag';elseif falta==5 cfalta='Bcg';elseif falta==6 cfalta='Cbg';
elseif falta==7 cfalta='Aab';elseif falta==8 cfalta='Bac';elseif falta==9 cfalta='Cbc';
elseif falta==10 cfalta='Aabg';elseif falta==11 cfalta='Bacg';elseif falta==12 cfalta='Cbcg';
elseif falta==13 cfalta='Aabc';elseif falta==14 cfalta='Aabcg';elseif falta==15 cfalta='ABCabc';
elseif falta==16 cfalta='ABCabcg';
end
% Conversion I(0,1,2)pu a I(0,1,2)Amp
[IIm0,IIm1,IIm2]=conv3(Im0,Im1,Im2,IbaseA);
[IIp0,IIp1,IIp2]=conv3(Ip0,Ip1,Ip2,IbaseB);
[IIs0,IIs1,IIs2]=conv3(Is0,Is1,Is2,IbaseA);
[IISAm0,IISAm1,IISAm2]=...
conv3(ISAm0,ISAm1,ISAm2,IbaseA);
[IISBp0,IISBp1,IISBp2]=...
conv3(ISBp0,ISBp1,ISBp2,IbaseB);
% Conversion V(0,1,2)pu a V(0,1,2)Volt
[VVm0,VVm1,VVm2,VVSA0,VVSA1,VVSA2]=...
conv6(Vm0,Vm1,Vm2,VSA0,VSA1,VSA2,VbaseA);
[VVp0,VVp1,VVp2,VVSB0,VVSB1,VVSB2]=...
conv6(Vp0,Vp1,Vp2,VSB0,VSB1,VSB2,VbaseB);
% Conversion I(a,b,c)pu a I(a,b,c)Amp
[IIma,IImb,IImc]=conv3(Ima,Imb,Imc,IbaseA);
[IIpa,IIpb,IIpc]=conv3(Ipa,Ipb,Ipc,IbaseB);
[IIsa,IIsb,IIsc]=conv3(Isa,Isb,Isc,IbaseA);
[IISAma,IISAmb,IISAmc]=...
conv3(ISAma,ISAmb,ISAmc,IbaseA);
[IISBpa,IISBpb,IISBpc]=...
conv3(ISBpa,ISBpb,ISBpc,IbaseB);
% Conversion V(a,b,c)pu a V(a,b,c)Volt
[VVma,VVmb,VVmc]=conv3(Vma,Vmb,Vmc,VbaseA);
[VVpa,VVpb,VVpc]=conv3(Vpa,Vpb,Vpc,VbaseB);
[VVSAa,VVSAb,VVSAc]=...
conv3(VSAa,VSAb,VSAc,VbaseA);
[VVSBa,VVSBb,VVSBc]=...
conv3(VSBa,VSBb,VSBc,VbaseB);
% 21: Conversion Z(pu) a Z(Ohm)
[SAm_ZZAB,SAm_ZZBC,SAm_ZZCA,SAm_ZZAG,SAm_ZZBG,SAm_ZZCG]=...
conv6(SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG,zbaseA);
[SBp_ZZAB,SBp_ZZBC,SBp_ZZCA,SBp_ZZAG,SBp_ZZBG,SBp_ZZCG]=...
conv6(SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG,zbaseB);
% 67N: Conversion pu a Volt
VVSApol=VSApol*VbaseA;
VVSBpol=VSBpol*VbaseB;
% 67N: Conversion pu a Amp
IISAmop=ISAmop*IbaseA;
IISBpop=ISBpop*IbaseB;
% -----------------------------------------------------------------------
% SALIDAS
% -----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'**********************************************************\n');
fprintf(fid,'UA =%-4.0f(kV)\n',UA);
fprintf(fid,'UB =%-4.0f(kV)\n',UB);
fprintf(fid,'LA=%-4.0f(km)\n',LA);
fprintf(fid,'LB=%-4.0f(km)\n',LB);
fprintf(fid,'kE=E/Enominal =%-4.0f\n',kE);
fprintf(fid,'Sccs=%-5.0f\n',Sccs);
fprintf(fid,'Ks0=%-5.2f\n',Ks0);
fprintf(fid,'Compensacion 21, Si(1), No(2) =%-1.0f\n',Compens21);
fprintf(fid,'\n');
fprintf(fid,'m =%-3.2f\n',m);
fprintf(fid,'Falta:%s\n',cfalta);
fprintf(fid,'Rf:%-10.2f\n',Rf);
if (falta==4)|(falta==5)|(falta==6)|(falta==10)|(falta==11)|(falta==12)|(falta==14)|(falta==16)
fprintf(fid,'Rg:%-10.2f\n',Rg);
else
fprintf(fid,'\n');
end
fprintf(fid,'**********************************************************\n');
variassalidas=1;
while variassalidas==1
fprintf(fid,'\n');
output=input('Salida: pu(cart)=1, pu(polar)=2, Ohm-Volt-Amp(cart)=3, Ohm-Volt-Amp(polar)=4 :');
if salida==2
fprintf(fid,'Salida: pu(cart)=1, pu(polar)=2, Ohm-Volt-Amp(cart)=3, Ohm-Volt-Amp(polar)=4 : %1.0f',output);
fprintf(fid,'\n');
end
switch(output)
% OPCIONES DE SALIDA
% --------------------------------
% 21 67N V(0,1,2) V(a,b,c)
% I(0,1,2) I(a,b,c)
% 1: pu(cartesiana) *** ***
% 2: pu(polar)
% 3: OVA(cartesiana) *** *** *** ***
% 4: OVA(polar) *** *** *** ***
% **********************************************************************
case{1}, % pu(cartesiana): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
% ----------------- -----------------
% **********************************************************************
% -----------------------------------------------------------------------
% C1.3 SALIDA V(0,1,2),I(0,1,2) pu cartesiana
% -----------------------------------------------------------------------
% Vm,Vp,VSA,VRA,VSB,VRB
% Im,Ip,Is,Ir,ISAm,IRAm,ISBp,IRBp,Izsr
% -----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(0,1,2),I(0,1,2) pu cartesiana\n');
fprintf(fid,'\n');
fprintf(fid,'Vm0 :%8.4f+i(%-8.4f)pu',real(Vm0),imag(Vm0));
fprintf(fid,'\tVm1 :%8.4f+i(%-8.4f)pu',real(Vm1),imag(Vm1));
fprintf(fid,'\tVm2 :%8.4f+i(%-8.4f)pu',real(Vm2),imag(Vm2));
fprintf(fid,'\n');
fprintf(fid,'Vp0 :%8.4f+i(%-8.4f)pu',real(Vp0),imag(Vp0));
fprintf(fid,'\tVp1 :%8.4f+i(%-8.4f)pu',real(Vp1),imag(Vp1));
fprintf(fid,'\tVp2 :%8.4f+i(%-8.4f)pu',real(Vp2),imag(Vp2));
fprintf(fid,'\n');
fprintf(fid,'VSA0 :%8.4f+i(%-8.4f)pu',real(VSA0),imag(VSA0));
fprintf(fid,'\tVSA1 :%8.4f+i(%-8.4f)pu',real(VSA1),imag(VSA1));
fprintf(fid,'\tVSA2 :%8.4f+i(%-8.4f)pu',real(VSA2),imag(VSA2));
fprintf(fid,'\n');
fprintf(fid,'VSB0 :%8.4f+i(%-8.4f)pu',real(VSB0),imag(VSB0));
fprintf(fid,'\tVSB1 :%8.4f+i(%-8.4f)pu',real(VSB1),imag(VSB1));
fprintf(fid,'\tVSB2 :%8.4f+i(%-8.4f)pu',real(VSB2),imag(VSB2));
fprintf(fid,'\n');
fprintf(fid,'Im0 :%8.2f+i(%-8.2f)',real(Im0),imag(Im0));
fprintf(fid,'\tIm1 :%8.2f+i(%-8.2f)',real(Im1),imag(Im1));
fprintf(fid,'\tIm2 :%8.2f+i(%-8.2f)',real(Im2),imag(Im2));
fprintf(fid,'\n');
fprintf(fid,'Ip0 :%8.2f+i(%-8.2f)',real(Ip0),imag(Ip0));
fprintf(fid,'\tIp1 :%8.2f+i(%-8.2f)',real(Ip1),imag(Ip1));
fprintf(fid,'\tIp2 :%8.2f+i(%-8.2f)',real(Ip2),imag(Ip2));
fprintf(fid,'\n');
fprintf(fid,'Is0 :%8.4f+i(%-8.4f)pu',real(Is0),imag(Is0));
fprintf(fid,'\tIs1 :%8.4f+i(%-8.4f)pu',real(Is1),imag(Is1));
fprintf(fid,'\tIs2 :%8.4f+i(%-8.4f)pu',real(Is2),imag(Is2));
fprintf(fid,'\n');
fprintf(fid,'ISAm0:%8.4f+i(%-8.4f)',real(ISAm0),imag(ISAm0));
fprintf(fid,'\tISAm1:%8.4f+i(%-8.4f)',real(ISAm1),imag(ISAm1));
fprintf(fid,'\tISAm2:%8.4f+i(%-8.4f)',real(ISAm2),imag(ISAm2));
fprintf(fid,'\n');
fprintf(fid,'ISBp0:%8.4f+i(%-8.4f)',real(ISBp0),imag(ISBp0));
fprintf(fid,'\tISBp1:%8.4f+i(%-8.4f)',real(ISBp1),imag(ISBp1));
fprintf(fid,'\tISBp2:%8.4f+i(%-8.4f)',real(ISBp2),imag(ISBp2));
fprintf(fid,'\n');
% ------------------------------------------------------------------
% C1.4 SALIDA V(a,b,c),I(a,b,c) pu cartesiana
% ------------------------------------------------------------------
% Vm,Vp,VSA,VRA,VSB,VRB
% Im,Ip,Is,Ir,ISAm,IRAm,ISBp,IRBp,Izsr
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(a,b,c),I(a,b,c) pu cartesiana\n\n');
fprintf(fid,'\n');
fprintf(fid,'Vma :%8.4f+i(%-8.4f)pu',real(Vma),imag(Vma));
fprintf(fid,'\tVmb :%8.4f+i(%-8.4f)pu',real(Vmb),imag(Vmb));
fprintf(fid,'\tVmc :%8.4f+i(%-8.4f)pu',real(Vmc),imag(Vmc));
fprintf(fid,'\n');
fprintf(fid,'Vpa :%8.4f+i(%-8.4f)pu',real(Vpa),imag(Vpa));
fprintf(fid,'\tVpb :%8.4f+i(%-8.4f)pu',real(Vpb),imag(Vpb));
fprintf(fid,'\tVpc :%8.4f+i(%-8.4f)pu',real(Vpc),imag(Vpc));
fprintf(fid,'\n');
fprintf(fid,'VSAa :%8.4f+i(%-8.4f)pu',real(VSAa),imag(VSAa));
fprintf(fid,'\tVSAb :%8.4f+i(%-8.4f)pu',real(VSAb),imag(VSAb));
fprintf(fid,'\tVSAc :%8.4f+i(%-8.4f)pu',real(VSAc),imag(VSAc));
fprintf(fid,'\n');
fprintf(fid,'VSBa :%8.4f+i(%-8.4f)pu',real(VSBa),imag(VSBa));
fprintf(fid,'\tVSBb :%8.4f+i(%-8.4f)pu',real(VSBb),imag(VSBb));
fprintf(fid,'\tVSBc :%8.4f+i(%-8.4f)pu',real(VSBc),imag(VSBc));
fprintf(fid,'\n');
fprintf(fid,'Ima :%8.4f+i(%-8.4f)pu',real(Ima),imag(Ima));
fprintf(fid,'\tImb :%8.4f+i(%-8.4f)pu',real(Imb),imag(Imb));
fprintf(fid,'\tImc :%8.4f+i(%-8.4f)pu',real(Imc),imag(Imc));
fprintf(fid,'\n');
fprintf(fid,'Ipa :%8.4f+i(%-8.4f)pu',real(Ipa),imag(Ipa));
fprintf(fid,'\tIpb :%8.4f+i(%-8.4f)pu',real(Ipb),imag(Ipb));
fprintf(fid,'\tIpc :%8.4f+i(%-8.4f)pu',real(Ipc),imag(Ipc));
fprintf(fid,'\n');
fprintf(fid,'Isa :%8.4f+i(%-8.4f)pu',real(Isa),imag(Isa));
fprintf(fid,'\tIsb :%8.4f+i(%-8.4f)pu',real(Isb),imag(Isb));
fprintf(fid,'\tIsc :%8.4f+i(%-8.4f)pu',real(Isc),imag(Isc));
fprintf(fid,'\n');
fprintf(fid,'ISAma:%8.4f+i(%-8.4f)',real(ISAma),imag(ISAma));
fprintf(fid,'\tISAmb:%8.4f+i(%-8.4f)',real(ISAmb),imag(ISAmb));
fprintf(fid,'\tISAmc:%8.4f+i(%-8.4f)',real(ISAmc),imag(ISAmc));
fprintf(fid,'\n');
fprintf(fid,'ISBpa:%8.4f+i(%-8.4f)',real(ISBpa),imag(ISBpa));
fprintf(fid,'\tISBpb:%8.4f+i(%-8.4f)',real(ISBpb),imag(ISBpb));
fprintf(fid,'\tISBpc:%8.4f+i(%-8.4f)',real(ISBpc),imag(ISBpc));
fprintf(fid,'\n');
% **********************************************************************
case{2}, % pu(polar): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
%
% **********************************************************************
% **********************************************************************
case{3}, % OVA(cartesiana): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
% -- --- ----------------- -----------------
% **********************************************************************
% ------------------------------------------------------------------
% C3.1 SALIDA 21 Ohm cartesiana
% ------------------------------------------------------------------
% SAm,RAm,SBp,RBp
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SAm =%-5.2f+i(%-5.2f)\n',real(SAm),imag(SAm));
fprintf(fid,'SARA =%-5.2f+i(%-5.2f)\n',real(SARA),imag(SARA));
fprintf(fid,'\n');
fprintf('p =%-3.2f\n',p);
fprintf('SBp = %5.2f+i(%-5.2f)\n',real(SBp),imag(SBp));
fprintf('SBRB = %5.2f+i(%-5.2f)\n',real(SBRB),imag(SBRB));
fprintf(fid,'\n');
[SAm_ZZABre,SAm_ZZABim]=zmed(SAm_ZZAB,1);
fprintf(fid,'SAm_ZAB=%8.2f+i(%-8.2f)(Ohm)\n',SAm_ZZABre,SAm_ZZABim);
[SAm_ZZBCre,SAm_ZZBCim]=zmed(SAm_ZZBC,1);
fprintf(fid,'SAm_ZBC=%8.2f+i(%-8.2f)(Ohm)(Ohm)\n',SAm_ZZBCre,SAm_ZZBCim);
[SAm_ZZCAre,SAm_ZZCAim]=zmed(SAm_ZZCA,1);
fprintf(fid,'SAm_ZCA=%8.2f+i(%-8.2f)(Ohm)\n',SAm_ZZCAre,SAm_ZZCAim);
[SAm_ZZAGre,SAm_ZZAGim]=zmed(SAm_ZZAG,1);
fprintf(fid,'SAm_ZAG=%8.2f+i(%-8.2f)(Ohm)\n',SAm_ZZAGre,SAm_ZZAGim);
[SAm_ZZBGre,SAm_ZZBGim]=zmed(SAm_ZZBG,1);
fprintf(fid,'SAm_ZBG=%8.2f+i(%-8.2f)(Ohm)\n',SAm_ZZBGre,SAm_ZZBGim);
[SAm_ZZCGre,SAm_ZZCGim]=zmed(SAm_ZZCG,1);
fprintf(fid,'SAm_ZCG=%8.2f+i(%-8.2f)(Ohm)\n',SAm_ZZCGre,SAm_ZZCGim);
fprintf('\n\n');
[SBp_ZZABre,SBp_ZZABim]=zmed(SBp_ZZAB,1);
fprintf(fid,'SBp_ZAB=%8.2f+i(%-8.2f)(Ohm)\n',SBp_ZZABre,SBp_ZZABim);
[SBp_ZZBCre,SBp_ZZBCim]=zmed(SBp_ZZBC,1);
fprintf(fid,'SBp_ZBC=%8.2f+i(%-8.2f)(Ohm)\n',SBp_ZZBCre,SBp_ZZBCim);
[SBp_ZZCAre,SBp_ZZCAim]=zmed(SBp_ZZCA,1);
fprintf(fid,'SBp_ZCA=%8.2f+i(%-8.2f)(Ohm)\n',SBp_ZZCAre,SBp_ZZCAim);
[SBp_ZZAGre,SBp_ZZAGim]=zmed(SBp_ZZAG,1);
fprintf(fid,'SBp_ZAG=%8.2f+i(%-8.2f)(Ohm)\n',SBp_ZZAGre,SBp_ZZAGim);
[SBp_ZZBGre,SBp_ZZBGim]=zmed(SBp_ZZBG,1);
fprintf(fid,'SBp_ZBG=%8.2f+i(%-8.2f)(Ohm)\n',SBp_ZZBGre,SBp_ZZBGim);
[SBp_ZZCGre,SBp_ZZCGim]=zmed(SBp_ZZCG,1);
fprintf(fid,'SBp_ZCG=%8.2f+i(%-8.2f)(Ohm)\n',SBp_ZZCGre,SBp_ZZCGim);
fprintf(fid,'\n');
% ------------------------------------------------------------------
% C3.2 SALIDA 67N OVA cartesiana
% ------------------------------------------------------------------
% SAm,RAm,SBp,RBp
% ------------------------------------------------------------------
% ------------------------------------------------------------------
% C3.3 SALIDA V(0,1,2),I(0,1,2) VA cartesiana
% ------------------------------------------------------------------
% Vm,Vp,VSA,VRA,VSB,VRB
% Im,Ip,Is,Ir,ISAm,IRAm,ISBp,IRBp,Izsr
% ------------------------------------------------------------------
fprintf(fid,'SALIDA V(0,1,2),I(0,1,2) OVA cartesiana\n');
fprintf(fid,'\n');
fprintf(fid,'VVm0:%8.0f+i(%-8.0f)',real(VVm0),imag(VVm0));
fprintf(fid,'\tVVm1:%8.0f+i(%-8.0f)',real(VVm1),imag(VVm1));
fprintf(fid,'\tVVm2:%8.0f+i(%-8.0f)\n',real(VVm2),imag(VVm2));
fprintf(fid,'VVp0:%8.0f+i(%-8.0f)',real(VVp0),imag(VVp0));
fprintf(fid,'\tVVp1:%8.0f+i(%-8.0f)',real(VVp1),imag(VVp1));
fprintf(fid,'\tVVp2:%8.0f+i(%-8.0f)',real(VVp2),imag(VVp2));
fprintf(fid,'\n');
fprintf(fid,'VVSA0:%8.0f+i(%-8.0f)',real(VVSA0),imag(VVSA0));
fprintf(fid,'\tVVSA1:%8.0f+i(%-8.0f)',real(VVSA1),imag(VVSA1));
fprintf(fid,'\tVVSA2:%8.0f+i(%-8.0f)\n',real(VVSA2),imag(VVSA2));
fprintf(fid,'VVSB0:%8.0f+i(%-8.0f)',real(VVSB0),imag(VVSB0));
fprintf(fid,'\tVVSB1:%8.0f+i(%-8.0f)',real(VVSB1),imag(VVSB1));
fprintf(fid,'\tVVSB2:%8.0f+i(%-8.0f)\n',real(VVSB2),imag(VVSB2));
fprintf(fid,'IIm0 :%8.0f+i(%-8.0f)',real(IIm0),imag(IIm0));
fprintf(fid,'\tIIm1 :%8.0f+i(%-8.0f)',real(IIm1),imag(IIm1));
fprintf(fid,'\tIIm2 :%8.0f+i(%-8.0f)\n',real(IIm2),imag(IIm2));
fprintf(fid,'IIp0 :%8.0f+i(%-8.0f)',real(IIp0),imag(IIp0));
fprintf(fid,'\tIIp1 :%8.0f+i(%-8.0f)',real(IIp1),imag(IIp1));
fprintf(fid,'\tIIp2 :%8.0f+i(%-8.0f)',real(IIp2),imag(IIp2));
fprintf(fid,'\n');
fprintf(fid,'IISAm0:%8.0f+i(%-8.0f)',real(IISAm0),imag(IISAm0));
fprintf(fid,'\tIISAm1:%8.0f+i(%-8.0f)',real(IISAm1),imag(IISAm1));
fprintf(fid,'\tIISAm2:%8.0f+i(%-8.0f)\n',real(IISAm2),imag(IISAm2));
fprintf(fid,'IISBp0:%8.0f+i(%-8.0f)',real(IISBp0),imag(IISBp0));
fprintf(fid,'\tIISBp1:%8.0f+i(%-8.0f)',real(IISBp1),imag(IISBp1));
fprintf(fid,'\tIISBp2:%8.0f+i(%-8.0f)\n',real(IISBp2),imag(IISBp2));
% ------------------------------------------------------------------
% C3.4 SALIDA V(a,b,c),I(a,b,c) VA cartesiana
% ------------------------------------------------------------------
% Vm,Vp,VSA,VRA,VSB,VRB
% Im,Ip,Is,Ir,ISAm,IRAm,ISBp,IRBp,Izsr
% ------------------------------------------------------------------
fprintf(fid,'SALIDA V(a,b,c),I(a,b,c) OVA cartesiana\n\n');
fprintf(fid,'VVma:%8.0f+i(%-8.0f)',real(VVma),imag(VVma));
fprintf(fid,'\tVVmb:%8.0f+i(%-8.0f)',real(VVmb),imag(VVmb));
fprintf(fid,'\tVVmc:%8.0f+i(%-8.0f)\n',real(VVmc),imag(VVmc));
fprintf(fid,'VVpa:%8.0f+i(%-8.0f)',real(VVpa),imag(VVpa));
fprintf(fid,'\tVVpb:%8.0f+i(%-8.0f)',real(VVpb),imag(VVpb));
fprintf(fid,'\tVVpc:%8.0f+i(%-8.0f)',real(VVpc),imag(VVpc));
fprintf(fid,'\n');
fprintf(fid,'VVSAa:%8.0f+i(%-8.0f)',real(VVSAa),imag(VVSAa));
fprintf(fid,'\tVVSAb:%8.0f+i(%-8.0f)',real(VVSAb),imag(VVSAb));
fprintf(fid,'\tVVSAc:%8.0f+i(%-8.0f)\n',real(VVSAc),imag(VVSAc));
fprintf(fid,'VVSBa:%8.0f+i(%-8.0f)',real(VVSBa),imag(VVSBa));
fprintf(fid,'\tVVSBb:%8.0f+i(%-8.0f)',real(VVSBb),imag(VVSBb));
fprintf(fid,'\tVVSBc:%8.0f+i(%-8.0f)\n',real(VVSBc),imag(VVSBc));
fprintf(fid,'IIma :%8.0f+i(%-8.0f)',real(IIma),imag(IIma));
fprintf(fid,'\tIImb :%8.0f+i(%-8.0f)',real(IImb),imag(IImb));
fprintf(fid,'\tIImc :%8.0f+i(%-8.0f)\n',real(IImc),imag(IImc));
fprintf(fid,'IIpa :%8.0f+i(%-8.0f)',real(IIpa),imag(IIpa));
fprintf(fid,'\tIIpb :%8.0f+i(%-8.0f)',real(IIpb),imag(IIpb));
fprintf(fid,'\tIIpc :%8.0f+i(%-8.0f)',real(IIpc),imag(IIpc));
fprintf(fid,'\n');
fprintf(fid,'IISAma:%8.0f+i(%-8.0f)',real(IISAma),imag(IISAma));
fprintf(fid,'\tIISAmb:%8.0f+i(%-8.0f)',real(IISAmb),imag(IISAmb));
fprintf(fid,'\tIISAmc:%8.0f+i(%-8.0f)\n',real(IISAmc),imag(IISAmc));
fprintf(fid,'IISBpa:%8.0f+i(%-8.0f)',real(IISBpa),imag(IISBpa));
fprintf(fid,'\tIISBpb:%8.0f+i(%-8.0f)',real(IISBpb),imag(IISBpb));
fprintf(fid,'\tIISBpc:%8.0f+i(%-8.0f)\n',real(IISBpc),imag(IISBpc));
% **********************************************************************
case{4}, % OVA(polar): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
%
% **********************************************************************
% ------------------------------------------------------------------
% C4.1 SALIDA 21, Ohm , polar
% ------------------------------------------------------------------
% SAm,RAm,SBp,RBp
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA 21, Ohm , polar \n');
fprintf(fid,'\n');
fprintf(fid,'Nota: Componentes de Z >=10000 Ohm se representan por Inf\n');
fprintf(fid,'Nota: Componentes de Z <=-10000 Ohm se representan por -Inf\n');
fprintf(fid,'\n');
[mod,arg]=mod_arg(SAm);fprintf(fid,'SAm=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
[mod,arg]=mod_arg(SARA);fprintf(fid,'SARA=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
[mod,arg]=mod_arg(SBp);fprintf(fid,'SBp=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
[mod,arg]=mod_arg(SBRB);fprintf(fid,'SBRB=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_argz(SAm_ZZAB);fprintf(fid,'SAm_ZAB=%7.2f(%-5.1fº)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(SAm_ZZBC);fprintf(fid,'SAm_ZBC=%7.2f(%-5.1fº)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(SAm_ZZCA);fprintf(fid,'SAm_ZCA=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(SAm_ZZAG);fprintf(fid,'SAm_ZAG=%7.2f(%-5.1fº)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(SAm_ZZBG);fprintf(fid,'SAm_ZBG=%7.2f(%-5.1fº)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(SAm_ZZCG);fprintf(fid,'SAm_ZCG=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
fprintf('\n\n');
[mod,arg]=mod_argz(SBp_ZZAB);fprintf(fid,'SBp_ZAB=%7.2f(%-5.1fº)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(SBp_ZZBC);fprintf(fid,'SBp_ZBC=%7.2f(%-5.1fº)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(SBp_ZZCA);fprintf(fid,'SBp_ZCA=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(SBp_ZZAG);fprintf(fid,'SBp_ZAG=%7.2f(%-5.1fº)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(SBp_ZZBG);fprintf(fid,'SBp_ZBG=%7.2f(%-5.1fº)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(SBp_ZZCG);fprintf(fid,'SBp_ZCG=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
fprintf('\n\n');
% -----------------------------------------------------------------------
% C4.2 SALIDA 67N Volt,Amp polar
% -----------------------------------------------------------------------
% SAm,RAm,SBp,RBp
% -----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA 67N, Volt, Amp, polar\n');
fprintf(fid,'\n');
fprintf(fid,'SAm_67N:\tVSApol=%-8.0f(V)\tISAmop=%-8.0f(A)\tAngDirSAm=%-5.1fº\n',VVSApol,IISAmop,AngDirSAm);
fprintf(fid,'SBp_67N:\tVSBpol=%-8.0f(V)\tISBpop=%-8.0f(A)\tAngDirSBp=%-5.1fº\n',VVSBpol,IISBpop,AngDirSBp);
% ------------------------------------------------------------------
% C4.3 SALIDA V(0,1,2), I(0,1,2) Volt,Amp polar
% ------------------------------------------------------------------
% Vm,Vp,VSA,VRA,VSB,VRB
% Im,Ip,Is,Ir,ISAm,IRAm,ISBp,IRBp,Izsr
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(0,1,2), I(0,1,2) pu polar\n');
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVm0);fprintf(fid,'Vm0 =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVm1);fprintf(fid,'Vm1 =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVm2);fprintf(fid,'Vm2 =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVp0);fprintf(fid,'Vp0 =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVp1);fprintf(fid,'Vp1 =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVp2);fprintf(fid,'Vp2 =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVSA0);fprintf(fid,'VSA0 =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVSA1);fprintf(fid,'VSA1 =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVSA2);fprintf(fid,'VSA2 =%6.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVSB0);fprintf(fid,'VSB0 =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVSB1);fprintf(fid,'VSB1 =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVSB2);fprintf(fid,'VSB2 =%6.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
[mod,arg]=mod_arg(IIm0);fprintf(fid,'Im0 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIm1);fprintf(fid,'Im1 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIm2);fprintf(fid,'Im2 =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIp0);fprintf(fid,'Ip0 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIp1);fprintf(fid,'Ip1 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIp2);fprintf(fid,'Ip2 =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIs0);fprintf(fid,'Is0 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIs1);fprintf(fid,'Is1 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIs2);fprintf(fid,'Is2 =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IISAm0);fprintf(fid,'ISAm0=%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IISAm1);fprintf(fid,'ISAm1=%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IISAm2);fprintf(fid,'ISAm2=%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IISBp0);fprintf(fid,'ISBp0=%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IISBp1);fprintf(fid,'ISBp1=%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IISBp2);fprintf(fid,'ISBp2=%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
% ------------------------------------------------------------------
% C4.4 SALIDA V(a,b,c), I(a,b,c) Volt,Amp polar
% ------------------------------------------------------------------
% Vm,Vp,VSA,VRA,VSB,VRB
% Im,Ip,Is,Ir,ISAm,IRAm,ISBp,IRBp,Izsr
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(a,b,c), I(a,b,c) Volt,Amp Polar\n');
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVma);fprintf(fid,'Vma =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVmb);fprintf(fid,'Vmb =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVmc);fprintf(fid,'Vmc =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVpa);fprintf(fid,'Vpa =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVpb);fprintf(fid,'Vpb =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVpc);fprintf(fid,'Vpc =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVSAa);fprintf(fid,'VSAa =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVSAb);fprintf(fid,'VSAb =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVSAc);fprintf(fid,'VSAc =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVSBa);fprintf(fid,'VSBa =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVSBb);fprintf(fid,'VSBb =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVSBc);fprintf(fid,'VSBc =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
[mod,arg]=mod_arg(IIma);fprintf(fid,'Ima =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IImb);fprintf(fid,'Imb =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IImc);fprintf(fid,'Imc =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIpa);fprintf(fid,'Ipa =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIpb);fprintf(fid,'Ipb =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIpc);fprintf(fid,'Ipc =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIsa);fprintf(fid,'Isa =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsb);fprintf(fid,'Isb =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsc);fprintf(fid,'Isc =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IISAma);fprintf(fid,'ISAma=%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IISAmb);fprintf(fid,'ISAmb=%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IISAmc);fprintf(fid,'ISAmc=%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IISBpa);fprintf(fid,'ISBpa=%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IISBpb);fprintf(fid,'ISBpb=%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(IISBpc);fprintf(fid,'ISBpc=%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
end % end del switch(output)
fprintf(fid,'\n');
variassalidas=input('Mas Salidas=1,Acabar=2 :');
end % end del while variassalidas
if salida==2
fclose(fid);
end
fprintf(1,'\n');
variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :');
end % end del while variasejecuciones
fprintf(1,'Pulse Enter\n');
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
