I'm trying to calculate the effect of J3 on my orbit's parameter (a, e, I, $\Omega$, ...)
I start with the equation of the gravitational potential of the Earth :
$$U = \frac{\mu}{r}\left(1 - \sum^\infty_{l=1} \left(\frac{R}{r}\right)^lJ_lP_l(sin \phi)\right)$$
With
$P_l$: the l-th Legendre's polynomials
(r, $\lambda$, $\phi$) is my spherical referential
So for the J3 perturbation, I have (I keep only the term where $J_3$ appears) :
$$U_3 = \frac{\mu}{r} \left(\frac{R}{r}\right)^3J_3P_3(sin \phi)$$
To determine the effect of this, I have to use Lagrange equations :
And now, I have a problem: How can I rewrite $U_3$ in function of the parameter of my orbit: a, e, I, $\Omega$, $\omega$ and M. Does anyone know how to transform my equation in term os osculating elements? Thanks !
