This post goes over the functions in the CR3BP folder of my open source Matlab Astrodynamics Library as well as why Jacobi values drift in ODE45