function[z1]=subz1(zs1,zr1,zl1u,L,m,p); z1=zeros(6); z1(1,1)=zs1; % NS Elemento 1 z1(2,2)=zr1; % NR Elemento 2 z1(3,3)=zl1u*m*L; % Sm Elemento 3 z1(4,4)=zl1u*p*L; % Sp Elemento 4 z1(5,5)=zl1u*(1-m)*L; % Rm Elemento 5 z1(6,6)=zl1u*(1-p)*L; % Rp Elemento 6