With the large number of stars in a galaxy, the particles in the simulation model typically do not map one-to-one onto physical objects. Part of the investigation is therefore to establish how many particles are necessary to simulate the specific macroscopic behaviour of the N-body system. The mathematical model adopter in this research is a model of gravitational forces between N particles which are integrated without any simplifying assumptions (such as the restriction of the radius of interaction). However, the contributions of general relativity and the finite speed of gravity is being ignored, because the dynamical timescales used in simulation are long compared to the light crossing times. Also, the space-time curvature induced by the particles and the particle velocities are small when averaged over the simulation time-step. The boundary conditions of the simulations are defined as periodic (or toroidal), so that one edge of the simulation volume matches up with the opposite edge.