function[z0]=subz0T(zs0,zr0,zsr0,zl0u,zlM0u,L,m,p)
z0=zeros(7);
z0(1,1)=zs0; % NS Elemento 1
z0(2,2)=zr0; % NR Elemento 2
z0(3,3)=zl0u*m*L; % Sm Elemento 3
z0(4,4)=zl0u*p*L; % Sp Elemento 4
z0(5,5)=zl0u*(1-m)*L; % Rm Elemento 5
z0(6,6)=zl0u*(1-p)*L; % Rp Elemento 6
z0(7,7)=zsr0; % zsr0 Elemento 7
z0(3,4)=zlM0u*m*L; % Mutua 3-4 Sm-Sp
z0(4,5)=-zlM0u*(p-m)*L; % Mutua 4-5 Sp-Rm
z0(5,6)=zlM0u*(1-p)*L; % Mutua 6-7 Rm-Rp
z0(4,3)=z0(3,4);
z0(5,4)=z0(4,5);
z0(6,5)=z0(5,6);
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
