Initial commit for Roche, oribital simulator in JavaScript

It performs poorly in both terms of accuracy and speed.
Add a formulation document of the simulation
Add midpoint method or second-order Runge-Kutta method

Now the orbits look more stable.