  We present a state-of-the-art N-body code which includes a detailed treatment
of stellar and binary evolution as well as the cluster dynamics. This code is
ideal for investigating all aspects relating to the evolution of star clusters
and their stellar populations. It is applicable to open and globular clusters
of any age. We use the N-body code to model the blue straggler population of
the old open cluster M67. Preliminary calculations with our binary population
synthesis code show that binary evolution alone cannot explain the observed
numbers or properties of the blue stragglers. On the other hand, our N-body
model of M67 generates the required number of blue stragglers and provides
formation paths for all the various types found in M67. This demonstrates the
effectiveness of the cluster environment in modifying the nature of the stars
it contains and highlights the importance of combining dynamics with stellar
evolution. We also perform a series of N = 10000 simulations in order to
quantify the rate of escape of stars from a cluster subject to the Galactic
tidal field.
