Probing Yukawian gravitational potential by numerical simulations. I.
  Changing N-body codes
  In the weak field limit general relativity reduces, as is well known, to the
Newtonian gravitation. Alternative theories of gravity, however, do not
necessarily reduce to Newtonian gravitation; some of them, for example, reduce
to Yukawa-like potentials instead of the Newtonian potential. Since the
Newtonian gravitation is largely used to model with success the structures of
the universe, such as for example galaxies and clusters of galaxies, a way to
probe and constrain alternative theories, in the weak field limit, is to apply
them to model the structures of the universe. In the present study, we consider
how to probe Yukawa-like potentials using N-body numerical simulations.
