Community Tip - Did you get called away in the middle of writing a post? Don't worry you can find your unfinished post later in the Drafts section of your profile page. X
Description:
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.
Odesolve: method Radau
Same initial condition but the method is Adams/BDF
Same initial condition but the method is Fixed