function[z0]=subz0(zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0,zpTS,zsTS,ztTS,zpTR,zsTR,ztTR) z0=zeros(12); z0(1,1)=zs0; z0(2,2)=zr0; z0(3,3)=zSAm0; z0(3,7)=zMSm0; z0(4,4)=zRAm0; z0(4,8)=zMRm0; z0(5,5)=ztTS; z0(6,6)=ztTR; z0(7,7)=zSBp0; z0(7,3)=z0(3,7); z0(8,8)=zRBp0; z0(8,4)=z0(4,8); z0(9,9)=zpTS; z0(10,10)=zsTS; z0(11,11)=zpTR; z0(12,12)=zsTR;