% [z1]=subz1(zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1,zpTS,zsTS,zpTR,zsTR); function [z1]=subz1(zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1,zpTS,zsTS,zpTR,zsTR); z1=zeros(8); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=zSAm1; z1(4,4)=zRAm1; z1(5,5)=zpTS+zsTS; z1(6,6)=zpTR+zsTR; z1(7,7)=zSBp1; z1(8,8)=zRBp1;