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