Simulated orbits of the two-body system projected onto the xy plane for the mass ratio m₁/m₂ = 3/2. The mesh has a size of 128³, and there are 100 time steps, each with four iterations for the MOND force. Red (m₁) and cyan (m₂) are the orbits obtained from numerical integration of Eqs. (19). Purple and blue are the orbits calculated with the PM code. The difference between the PM code and the ODEs is a slow dephasing of the orbits.