working on initial conditions. Spherical orbit obtained to good accuracy. Added a function to compute H derivatives numerically with flux set to 0.