A dynamic positron emission tomography problem to reconstruct temporally changing mass distributions is considered. We derive the reconstruction method in detail, analyze the associated optimization problem, and provide information on the numerical implementation.