Angular momentum vs. time t for the simulation shown in Fig. 2. These are constants of motion of the ODE system Eq. (19) and numerically constant (with machine precision, not shown) when integrated with Leapfrog. In the PM code, the total angular momentum (blue) is well conserved, but the single-particle angular momenta (orange and red) have large errors when the distance to the origin is smaller than one pixel.