Add Python version of the simulation

It's initially using Euler method.
It's not particularly faster than JavaScript, but it's easier to write
codes with overloaded arithmetic operators for vectors.
Add Second-order Runge-Kutta method to roche-glut.py