function[z2]=subz2(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r) z2=zeros(24); z2(1,1)=zs2; z2(2,2)=zx2; z2(3,3)=zr2; z2(4,4)=zy2; z2(5,5)=zl2u*a*L1; z2(6,6)=zl2u*r*L1; z2(7,7)=zl2u*(1-b)*L4; z2(8,8)=zl2u*(1-c)*L4; z2(9,9)=zl2u*(1-e)*L3; z2(10,10)=zl2u*(1-f)*L3; z2(11,11)=zl2u*(1-h)*L5; z2(12,12)=zl2u*(1-g)*L5; z2(13,13)=zl2u*(1-a)*L1; z2(14,14)=zl2u*b*L4; z2(15,15)=zl2u*c*L4; z2(16,16)=zl2u*d*L2; z2(17,17)=zl2u*(1-d)*L2; z2(18,18)=zl2u*e*L3; z2(19,19)=zl2u*f*L3; z2(20,20)=zl2u*g*L5; z2(21,21)=zl2u*h*L5; z2(22,22)=zl2u*(1-q)*L2; z2(23,23)=zl2u*q*L2; z2(24,24)=zl2u*(1-r)*L1;