% PG52n function[z1]=subz1(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r) z1=zeros(24); z1(1,1)=zs1; z1(2,2)=zx1; z1(3,3)=zr1; z1(4,4)=zy1; z1(5,5)=zl1u*a*L1; z1(6,6)=zl1u*r*L1; z1(7,7)=zl1u*(1-b)*L4; z1(8,8)=zl1u*(1-c)*L4; z1(9,9)=zl1u*(1-e)*L3; z1(10,10)=zl1u*(1-f)*L3; z1(11,11)=zl1u*(1-h)*L5; z1(12,12)=zl1u*(1-g)*L5; z1(13,13)=zl1u*(1-a)*L1; z1(14,14)=zl1u*b*L4; z1(15,15)=zl1u*c*L4; z1(16,16)=zl1u*d*L2; z1(17,17)=zl1u*(1-d)*L2; z1(18,18)=zl1u*e*L3; z1(19,19)=zl1u*f*L3; z1(20,20)=zl1u*g*L5; z1(21,21)=zl1u*h*L5; z1(22,22)=zl1u*(1-q)*L2; z1(23,23)=zl1u*q*L2; z1(24,24)=zl1u*(1-r)*L1;