In a quantum system damping and noise enter through the coupling with an external large bath. The usual treatments computes the density matrix of the system by expanding in the coupling strength to the bath and averaging over the bath. The result, however, is known to be unacceptable as it leads to negative probabilities. It is here proposed to study correlation functions rather than the density matrix itself. Equations for them are derived and compared with the result obtained by the usual method.