VarEqP:

Path: Orbit/OrbitMechanics

% Computes the right hand side of the equinoctial variational equations.
 Based on the equinoctial formulation with the acceleration
 in polar coordinates. These equations are singular
 when i = ? 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.
--------------------------------------------------------------------------

Children:

Orbit: OrbitCoord/K2L
Orbit: OrbitCoord/ML2K

Back to the Orbit Module page