function[Im0,Im1,Im2,Vm0,Vm1,Vm2]=mag2m(Z0,Z1,Z2,E,faltam,Rfm,Rgm,pfaltam);
a=cos(2*pi/3)+sin(2*pi/3)*i;
Rgm=0;
Mm=zeros(6);
%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 :');
switch(faltam)
case{1},
%Rfm=input('Resistencia de Falta en m. Rf(Ohm) = ');
Mm(1,1)=1;Mm(1,2)=-1;
Mm(2,1)=1;Mm(2,3)=-1;
Mm(3,1)=-3*Rfm;Mm(3,4)=1;Mm(3,5)=1;Mm(3,6)=1;
case{2},
%Rfm=input('Resistencia de Falta en m. Rf(Ohm): ');
Mm(1,1)=1;Mm(1,2)=-a^2;
Mm(2,1)=1;Mm(2,3)=-a;
Mm(3,1)=-Rfm;Mm(3,2)=-a^2*Rfm;Mm(3,3)=-a*Rfm;Mm(3,4)=1;Mm(3,5)=a^2;Mm(3,6)=a;
case{3},
%Rfm=input('Resistencia de Falta en m. Rf(Ohm): ');
Mm(1,1)=1;Mm(1,2)=-a;
Mm(2,1)=1;Mm(2,3)=-a^2;
Mm(3,1)=-Rfm;Mm(3,2)=-a*Rfm;Mm(3,3)=-a^2*Rfm;Mm(3,4)=1;Mm(3,5)=a;Mm(3,6)=a^2;
case{4},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');
Mm(1,1)=1;
Mm(2,2)=1;Mm(2,3)=a;
Mm(3,2)=-Rfm;Mm(3,3)=a*Rfm;Mm(3,5)=1;Mm(3,6)=-a;
case{5},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');
Mm(1,1)=1;
Mm(2,2)=1;Mm(2,3)=1;
Mm(3,2)=-Rfm;Mm(3,3)=Rfm;Mm(3,5)=1;Mm(3,6)=-1;
case{6},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');
Mm(1,1)=1;
Mm(2,2)=1;Mm(2,3)=a^2;
Mm(3,2)=-Rfm;Mm(3,3)=a^2*Rfm;Mm(3,5)=1;Mm(3,6)=-a^2;
case{7},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');
%Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): ');
Mm(1,1)=1;Mm(1,2)=a;Mm(1,3)=a^2;
Mm(2,2)=-a^2*Rfm;Mm(2,3)=Rfm;Mm(2,5)=a^2;Mm(2,6)=-1;
Mm(3,1)=-(Rfm+3*Rgm);Mm(3,2)=-Rfm;Mm(3,3)=-Rfm;
Mm(3,4)=1;Mm(3,5)=1;Mm(3,6)=1;
case{8},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');
%Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): ');
Mm(1,1)=1;Mm(1,2)=1;Mm(1,3)=1;
Mm(2,2)=-Rfm;Mm(2,3)=Rfm;Mm(2,5)=1;Mm(2,6)=-1;
Mm(3,1)=-(Rfm+3*Rgm);Mm(3,2)=-a^2*Rfm;Mm(3,3)=-a*Rfm;
Mm(3,4)=1;Mm(3,5)=a^2;Mm(3,6)=a;
case{9},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');
%Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): ');
Mm(1,1)=1;Mm(1,2)=a^2;Mm(1,3)=a;
Mm(2,2)=-Rfm;Mm(2,3)=a^2*Rfm;Mm(2,5)=1;Mm(2,6)=-a^2;
Mm(3,1)=-(Rfm+3*Rgm);Mm(3,2)=-a*Rfm;Mm(3,3)=-a^2*Rfm;
Mm(3,4)=1;Mm(3,5)=a;Mm(3,6)=a^2;
case{10},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm) = ');
Mm(1,1)=1;
Mm(2,2)=Rfm*a^2;Mm(2,3)=-Rfm;Mm(2,5)=-a^2;Mm(2,6)=1;
Mm(3,2)=-Rfm;Mm(3,3)=Rfm*a^2;Mm(3,5)=1;Mm(3,6)=-a^2;
case{11},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm) = ');
%Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): ');
Mm(1,1)=-(Rfm+3*Rgm);Mm(1,4)=1;
Mm(2,2)=-Rfm;Mm(2,3)=-Rfm;Mm(2,5)=1;Mm(2,6)=1;
Mm(3,2)=-a*Rfm;Mm(3,3)=-Rfm;Mm(3,5)=a;Mm(3,6)=1;
end
Mm(4,1)=Z0(pfaltam,pfaltam);Mm(4,4)=1;
Mm(5,2)=Z1(pfaltam,pfaltam);Mm(5,5)=1;
Mm(6,3)=Z2(pfaltam,pfaltam);Mm(6,6)=1;
P=[0;0;0;0;E;0];
N=Mm\P;
Im0=N(1);Im1=N(2);Im2=N(3);
Vm0=N(4);Vm1=N(5);Vm2=N(6);
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
