function[z2]=subz2(zs2,zr2,zt2,zl2u,L1,L2,L3,a,b,c,d) z2=zeros(11); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zt2; z2(4,4)=zl2u*a*L1; z2(5,5)=zl2u*c*L1; z2(6,6)=zl2u*(1-b)*L2; z2(7,7)=zl2u*(1-d)*L3; z2(8,8)=zl2u*(1-a)*L1; z2(9,9)=zl2u*b*L2; z2(10,10)=zl2u*(1-c)*L1; z2(11,11)=zl2u*d*L3;