% PC28Mb function [z0]=subz0(zs0,zSAm0,zSBp0,zMSm0,zpTS,zsTS,ztTS) z0=zeros(4); z0(1,1)=zs0; z0(2,2)=zSAm0; z0(2,3)=zMSm0; z0(3,3)=zSBp0; z0(3,2)=z0(2,3); z0(4,4)=ztTS+zsTS;