Community Tip - Need to share some code when posting a question or reply? Make sure to use the "Insert code sample" menu option. Learn more! 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