% [z2]=subz2(zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2,zpTS,zsTS,zpTR,zsTR); function [z2]=subz2(zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2,zpTS,zsTS,zpTR,zsTR) z2=zeros(8); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zSAm2; z2(4,4)=zRAm2; z2(5,5)=zpTS+zsTS; z2(6,6)=zpTR+zsTR; z2(7,7)=zSBp2; z2(8,8)=zRBp2;