An Analysis of Numerical Integration Techniques for the Modeling of Open Galactic Clusters.

While testing the integration technique, it was found that one of the constants of integration, the energy, was not conserved. This non-conservation of energy produces an enhancement of the evolutionary trends. Therefore, this integration scheme, is not applicable to the problem of integrating individually tailored evolutionary models. Since these individualized models based on the observational data could not be produced, an analysis was done of the errors still inherent in a numerical integration of the N-body problem using Wielen's polynomial technique. It was found that the energy error drives the evolution to form a denser nucleus and a more extended halo than would occur in a real cluster.