Should be solving the variational equations correctly; testing of d^2Hdq^2 done in orbit_lin_test
lyapunov should be solving linearization correctly; dealing with overflow issues; lyap exp not correct