  The past few years have seen dramatic improvements in the scope and realism
of star cluster simulations. Accurate treatments of stellar evolution, coupled
with robust descriptions of all phases of binary evolution, have been
incorporated self-consistently into several dynamical codes, allowing for the
first time detailed study of the interplay between stellar dynamics and stellar
physics. The coupling between evolution, dynamics, and the observational
appearance of the cluster is particularly strong in young systems and those
containing large numbers of primordial binary systems, and important inroads
have been made in these areas, particularly in N-body simulations. I discuss
some technical aspects of the current generation of N-body integrators, and
describe some recent results obtained using these codes.
