N-Body Simulation and Visualization of Stellar Clusters

We have constructed a gravitational N-body simulation to study the evolution of stellar clusters. Our integration is based on the Hermite Individual Timestep Scheme (HITS) which accurately follows close interactions and collisions, and conserves energy to within acceptable limits. We discuss preliminary results of simulations that include a stellar mass distribution based on a three power law initial mass function. We also present a visualization technique that can be used as a tool in analyzing simulations and discuss plans for future studies.