%[z1]=subz1(zs1,zr1,zt1,zl1,m); function[z1]=subz1(zs1,zr1,zt1,zl1,m); z1=zeros(6); z1(1,1)=zs1; % NS Elemento 1 z1(2,2)=zr1; % NR Elemento 2 z1(3,3)=zt1; % NT Elemento 3 z1(4,4)=m*zl1; % Sm Elemento 3 z1(5,5)=(1-m)*zl1; % Rm Elemento 4 z1(6,6)=zl1; % ST Elemento 6