A general formulation of space-energy-time iterative synthesis in neutron diffusion theory is developed. A description of the iterative synthesis procedure is presented, followed by applications to the steady-state and time-dependent diffusion equation. It is shown how spatially and/or temporally discontinuous shape functions can be included. Finally, application is made to the discontinuous energy (multigroup) diffusion theory model.