%[z0]=subz0(zs0,zr0,zl0,zlM0,m); function[z0]=subz0(zs0,zr0,zl0,zlM0,m) z0=zeros(5); z0(1,1)=zs0; % NS Elemento 1 z0(2,2)=zr0; % NR Elemento 2 z0(3,3)=m*zl0; % Sm Elemento 3 z0(4,4)=zl0; % SR Elemento 4 z0(5,5)=(1-m)*zl0; % Rm Elemento 5 z0(3,4)=m*zlM0; % Mutua 3-4 Sm-SR z0(4,5)=-(1-m)*zlM0; % Mutua 4-5 SR-Rm z0(4,3)=z0(3,4); z0(5,4)=z0(4,5);