I'm trying to take into account the Earth geopotential effect on the orbit. The only term I could consider, is J2.
3.Why does the shown equation give 1 km error in 1 day in comparison with GMAT (configuration is below)? The coefficients for Eq radius and J2 are the same. The integrator is 9th order Runge-Kutta with tol=10e-13.