A method based on Jacobi Integral variational equation for computing Earth-Moon trajectories in the four-body problem. (December 2019)