N-MODY: a code for collisionless N-body simulations in modified   Newtonian dynamics
  We describe the numerical code N-MODY, a parallel particle-mesh code for collisionless N-body simulations in modified Newtonian dynamics (MOND). N-MODY is based on a numerical potential solver in spherical coordinates that solves the non-linear MOND field equation, and is ideally suited to simulate isolated stellar systems. N-MODY can be used also to compute the MOND potential of arbitrary static density distributions. A few applications of N-MODY indicate that some astrophysically relevant dynamical processes are profoundly different in MOND and in Newtonian gravity with dark matter. 
Pasquale Londrillo (1), Carlo Nipoti (2) ((1) INAF-Bologna   Astronomical Observatory, (2) Bologna University)