  A method is presented for numerical implementation of the extended TDHF
theory in which two-body correlations beyond the mean-field approximation are
incorporated in the form of a quantal collision term. The method is tested in a
model problem in which the exact solution can be obtained numerically. Whereas
the usual TDHF fails to reproduce the long time evolution, a very good
agreement is found between the extended TDHF and the exact solution.
