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