In this simulation of three gravitating Newtonian particles, a blue "planet" orbits a red "star" until they meet the "dark star." The physics is simple, yet the numerical behavior is extreme. Indeed, Mathcad's numerical integrators have difficulty modeling the motions correctly.
Solution is sensitive to the choosing numerical methods for solving
Same initial condition but the method is not RADAU (see the first video) but Adams/BDF: