Path: Orbit/OrbitMechanics
% Computes the right hand side of the Gauss variational equations.
 Based on the Gauss formulation with the acceleration in polar
 coordinates. These equations are singular when i = 0.
--------------------------------------------------------------------------
   Form:
   dElDt = VarGausP( el, aD, tol, mu )
--------------------------------------------------------------------------
   ------
   Inputs
   ------
   el                    Elements vector [a,i,W,w,e,M]
   aD                    Disturbance vector [r,theta,z] (default is zeros)
   tol                   Tolerance for Kepler's equation solver (0.0001)
   mu                    Gravitational parameter (3.98600436e5;)
   -------
   Outputs
   -------
   dElDt                 Derivative to the elements
--------------------------------------------------------------------------
   References: Battin, R.H., An Introduction to the Methods and Mathematics
               of Astrodynamics, AIAA, New York, 1987, pp. 484-490.
--------------------------------------------------------------------------
SC: BasicOrbit/E2Nu SC: BasicOrbit/M2EEl
Back to the Orbit Module page