  We present a rigorous framework to obtain evolution equations for the
momentum distribution and higher order correlation functions in weakly
interacting systems based on the Projection Operator Technique. These equations
can be numerically solved in an efficient way. We compare the solution of the
equations with known results for 1D models and find an excellent agreement.
