Matlab code for a three body gravitational simulation. This code uses ODE45 to solve the 18 differential equations in a three body problem. Most initial conditions result in chaotic results so a few solutions or orbits have been provided.
Chaotic Result
Most initial conditions will result in chaotic paths. Solutions are hard to predict and are very complex.
Ying-Yang Solution
A few unique and beautiful solution paths exist. Note that these are real, possible orbits that might exist somewhere in the universe!
