function[z0]=subz0(zs0,zr0,zl0,zlM0,m,p); z0=zeros(6); z0(1,1)=zs0; z0(2,2)=zr0; z0(3,3)=m*zl0; z0(4,4)=p*zl0; z0(5,5)=(1-m)*zl0; z0(6,6)=(1-p)*zl0; z0(3,4)=m*zlM0; z0(5,6)=(1-p)*zlM0; z0(4,3)=z0(3,4); z0(6,5)=z0(5,6);