LunarPolarOrbit:

Path: LunarMissions/LunarEphemeris

% Computes a polar lunar orbit
--------------------------------------------------------------------------
   Form:
     x = LunarPolarOrbit( a, e, jD )
--------------------------------------------------------------------------

   -------
   Inputs
   -------
   a  	(1,1)	Semi-major axis
   e  	(1,1)	Eccentricity
   jD	(1,1)	Julian date

   -------
   Outputs
   -------
   x  	(6,1) [r;v] (km,km/s) ECI

--------------------------------------------------------------------------