  The new approach outlined in Paper I (Spurzem \& Giersz 1996) to follow the
individual formation and evolution of binaries in an evolving, equal point-mass
star cluster is extended for the self-consistent treatment of relaxation and
close three- and four-body encounters for many binaries (typically a few
percent of the initial number of stars in the cluster). The distribution of
single stars is treated as a conducting gas sphere with a standard anisotropic
gaseous model. A Monte Carlo technique is used to model the motion of binaries,
their formation and subsequent hardening by close encounters, and their
relaxation (dynamical friction) with single stars and other binaries. The
results are a further approach towards a realistic model of globular clusters
with primordial binaries without using special hardware. We present, as our
main result, the self-consistent evolution of a cluster consisting of 300.000
equal point-mass stars, plus 30.000 equal mass binaries over several hundred
half-mass relaxation times, well into the phase where most of the binaries have
been dissolved and evacuated from the core. In a self-consistent model it is
the first time that such a realistically large number of binaries is evolving
in a cluster with an even ten times larger number of single stars. Due to the
Monte Carlo treatment of the binaries we can at every moment analyze their
external and internal parameters in the cluster as in an N-body simulation.
