added apply_jacobi_inv with dgetrs; ported and tested kepler_step
added symplectic method, direct and indirect terms; needs to be tested
implemented sparse FW/BW subs; hamiltonian is bounded