function[z2]=subz2(zs2,zr2,zt2,zl2,m,p) z2=zeros(7); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zt2; z2(4,4)=zl2*m; z2(5,5)=zl2*p; z2(6,6)=zl2*(1-m); z2(7,7)=zl2*(1-p);