function[z1]=subz1(zs1,zr1,zl1,m,p) z1=zeros(6); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=m*zl1; z1(4,4)=p*zl1; z1(5,5)=(1-m)*zl1; z1(6,6)=(1-p)*zl1;