Derivatives of eccentricity and argument of periapsis with respect to time
[d_ecc_dt,d_pom_dt] = CL_op_frozenOrbitDer(sma,ecc,inc,pom, [er,mu,zonals])
Computes the derivatives of eccentricity and argument of periapsis with respect to time that correspond to the gravitational effects due to J2 and J3.
Semi-major axis [m] (1x1 or 1xN)
Eccentricity (1x1 or 1xN)
Inclination [rad] (1x1 or 1xN)
Argument of pariapsis [rad] (1x1 or 1xN)
(optional) Equatorial radius [m] (default is %CL_eqRad)
(optional) Gravitational constant [m^3/s^2] (default value is %CL_mu)
(optional) Vector of zonal coefficients J1 to Jn, troncated to J3. Default is %CL_j1jn(1:3)). (1xNz)
Derivative of eccentricity [s^-1] (1xN)
Derivative of argument of periapsis [rad/s] (1xN)
CNES - DCT/SB
1) "Frozen orbits in the J2+J3 problem", Krystyna Kiedron and Richard Cook, AAS 91-426.