Merged in eschnett/rk2-central (pull request #6)
This is an RK2 integrator where the first substep has size dt/2.
This is an RK2 integrator where the first substep has size dt/2.