Three Dimensional Numerical Code for the Expanding Flat Universe

The current distribution of galaxies may contain clues to the condition of the universe when the galaxies condensed and to nature of the subsequent expansion of the universe. The development of this large scale structure can be studied by employing N-body computer simulations. The present paper describes the code developed for this purpose. The computer code calculates the motion of collisionless matter acting under the force of gravity in an expanding flat universe. The test run of the code shows the error less than 0.5% in 100 iterations.