Simple SciPy procedure for calculating newtonian dynamics of an n-body system - View it on GitHub