function[z1]=subz1_2sr(zs1,zsr1,zr1,zl1,m,p); z1=zeros(8); z1(1,1)=zs1; % NS Elemento 1 z1(2,2)=zr1; % NR Elemento 2 z1(3,3)=zl1*m; % Sm Elemento 3 z1(4,4)=zl1*p; % Sp Elemento 4 z1(5,5)=zl1*(1-m); % Rm Elemento 5 z1(6,6)=zl1*(1-p); % Rp Elemento 6 z1(7,7)=-z1(3,3); % Sm Elemento 7 z1(8,8)=zsr1;