% PC62T %[z0]=subz0(zs0,zr0,zl0,zlM0,m); function [z0,P0,Q0]=subz0(zs0,zr0,zl0,zlM0,m) zz0=zeros(5); zz0(1,1)=zs0; % NS Elemento 1 zz0(2,2)=zr0; % NR Elemento 2 zz0(3,3)=m*zl0; % Sm Elemento 3 zz0(4,4)=(1-m)*zl0; % Rm Elemento 4 zz0(5,5)=zl0; % Lazo Elemento 5 zz0(3,5)=-m*zlM0; % Mutua 3-5 Sm-Lazo zz0(4,5)=(1-m)*zlM0; % Mutua 4-5 Rm-Lazo zz0(5,3)=zz0(3,5); zz0(5,4)=zz0(4,5); M0=zz0(1:4,1:4); N0=zz0(1:4,5:5); P0=zz0(5:5,1:4); Q0=zz0(5:5,5:5); z0=M0-N0*inv(Q0)*P0;