N-body realizations of compound galaxies

A prescription for constructing N-body models of galaxies consisting of more than one component is described. Spatial density profiles are realized exactly, but the phase space distribution is approximated using moments of the collisionless Boltzmann equation. While this approach is not fully rigorous, empirical tests suggest that it is adequate for studies of, e.g., interacting galaxies and the forced response of galaxies to imposed perturbations such as bars.