  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.
