We develop a theory for the dynamics of the density matrix describing a multimode polariton condensate. In such a condensate several single-particle orbitals become highly occupied, due to stimulated scattering from reservoirs ...