function[z1]=subz1T4(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,m,p) z1=zeros(10); z1(1,1)=zs1; z1(2,2)=zx1; z1(3,3)=zr1; z1(4,4)=zy1; z1(5,5)=zl1u*(L1+m*L4); z1(6,6)=zl1u*(L1+L2+L5); z1(7,7)=zl1u*(1-m)*L4; z1(8,8)=zl1u*(L5+L3); z1(9,9)=zl1u*(L3+L2+p*L4); z1(10,10)=zl1u*(1-p)*L4;