Removed old constraint interface
Refactored the system so it works with indices in place of real particles
Added RK solver and Mid-Point Solver
added extra keys for change in timestep and solver choice