Community Tip - Did you get an answer that solved your problem? Please mark it as an Accepted Solution so others with the same problem can find the answer easily. 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