calculate velocity of center of mass alongside its position
merge PlanetarySystem.centerOfMass into PlanetarySystem.center
add subtractCenterOfMass function; subtract center of mass on calls to PlanetarySystem.reset
fix bugs in subtractCenterOfMass