For a discussion of the difficulties in making numerical simulations with MOND, see this humorous write-up by Prof. Mihos and this brief but more formal discussion.
Rotation curves of the spiral galaxy NGC 6946 (left)
and the gas rich dwarf NGC 1560 (right). The black line is the mass model
for the stars and gas. The blue line is the MOND fit.