function[z0]=subz0_2sr(zs0,zsr0,zr0,zl0,zlM0,m,p)
z0=zeros(8);
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*p; % Sp Elemento 4
z0(5,5)=zl0*(1-m); % Rm Elemento 5
z0(6,6)=zl0*(1-p); % Rp Elemento 6
z0(7,7)=-z0(3,3); % Sm Elemento 3
z0(8,8)=zsr0;
if m<=p
z0(3,4)=zlM0*m; % Mutua 3-4 Sm-Sp
z0(4,5)=-zlM0*(p-m); % Mutua 4-5 Sp-Rm
z0(5,6)=zlM0*(1-p); % Mutua 5-6 Rm-Rp
z0(7,4)=z0(3,4); % Mutua 7-4 Sm-Sp
z0(4,3)=z0(3,4);
z0(5,4)=z0(4,5);
z0(6,5)=z0(5,6);
z0(4,7)=z0(7,4);
else
z0(3,4)=zlM0*p; % Mutua 3-4 Sm-Sp
z0(3,6)=-zlM0*(m-p); % Mutua 3-6 Sm-Rp
z0(5,6)=zlM0*(1-m); % Mutua 5-6 Rm-Rp
z0(7,4)=z0(3,4); % Mutua 7-4 Sm-Sp
z0(7,6)=z0(3,6); % Mutua 7-6 Sm-Rp
z0(4,3)=z0(3,4);
z0(6,3)=z0(3,6);
z0(6,5)=z0(5,6);
z0(4,7)=z0(7,4);
z0(6,7)=z0(7,6);
end
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
