!!! Not adapted to run on all machines, creates a numerical solution to the solar system for a given set of entities, each entity created with the 'body' class and added to the system state, a choice of integrators can be selected, the most advisable being the dynamic timestep that ensures the error on each step is never above a user defined value 'tolerance' inside integrators.py,