CL_op_frozenOrbitDer — Eccentricity and perigee drift over time for an orbit close to a frozen orbit
[d_ecc_dt,d_pom_dt] = CL_op_frozenOrbitDer(sma,ecc,inc,pom,[er,[mu,[zonals]]])
You can optionally give the zonal coefficients J1/J2/J3 to be used : default is %CL_j1jn(1:3).
semi major axis [m] (1x1 : same semi major axis for all) or (1xM)
eccentricity (1x1 : same eccentricity for all) or (1xM)
inclination [rad] (1x1 : same inclination for all) or (1xM)
argument of perigee (1x1 : same argument of perigee for all) or (1xM)
(optional) equatorial radius [m] (default is %CL_eqRad)
(optional) geocentric gravitational constant [m^3/s^2] (default value is %CL_mu)
(optional) vector of zonals coefficients J1 to Jn (troncated to J3) to be used (default is %CL_j1jn(1:3)) (1 x N)
(optional depending on model) eccentricity derivative with respect to time [s^-1] (1xM)
(optional depending on model) periapsis argument derivative with respect to time [rad/s] (1xM)