  A Monte Carlo simulation program is presented which can be used to determine
the small-$x$ evolution of a heavy onium using Mueller's colour dipole
formulation, giving the full distribution of dipoles in rapidity and impact
parameter. Routines are also provided which calculate onium-onium scattering
amplitudes between individual pairs of onium configurations, making it possible
to establish the contribution of multiple pomeron exchange terms to onium-onium
scattering (the unitarisation corrections).
