function [z2]=subz2(zs2,zSAm2,zSBp2,zpTS,zsTS) z2=zeros(4); z2(1,1)=zs2; z2(2,2)=zSAm2; z2(3,3)=zSBp2; z2(4,4)=zpTS+zsTS;