Path: Orbit/OrbitMechanics
% Computes the right hand side of the equinoctial variational equations (polar).
Based on the equinoctial formulation with the acceleration in polar
coordinates. These equations are singular when i = pi/2 and/or h = 0.
--------------------------------------------------------------------------
Form:
dEqDt = VarEqP( eq, aD, tol, mu )
--------------------------------------------------------------------------
------
Inputs
------
eq Elements vector [a,P1,P2,Q1,Q2,l]
aD Disturbance vector [r,theta,z]
tol Tolerance for the Kepler equation solver
mu Gravitational parameter
-------
Outputs
-------
dEqDt Derivative to the elements
--------------------------------------------------------------------------
References: Battin, R.H., An Introduction to the Methods and Mathematics
of Astrodynamics, AIAA, New York, 1987, pp. 492-493.
--------------------------------------------------------------------------
Orbit: OrbitCoord/K2L Orbit: OrbitCoord/ML2K
Back to the Orbit Module page