LocallyOptimalTraj:

--------------------------------------------------------------------------
   Computes locally optimal solar sail control angles.
   Constraints are applied to keep the cone angle within 35 degrees but
   they may be commented out.

   The four types of maneuvers are:

   'semi-major axis'
   'eccentricity'
   'inclination'
   'ascending node'

   Since version 7.
--------------------------------------------------------------------------
   Form:
   [alpha, delta, alphaHat] = LocallyOptimalTraj( type, r, v, mu, dirFlag, ascendingNode )
--------------------------------------------------------------------------

   ------
   Inputs
   ------
   type           (1,:)  Type (see above)
   r              (3,:)  Position vector (km)
   v              (3,:)  Velocity vector (km/s)
   mu             (1,1)  Gravitational parameter (km^3/s^2) 
                         (Default is the sun)
   dirFlag        (1,1)  Directional flag, +/- 1
   ascendingNode  (1,1)

   -------
   Outputs
   -------
   alpha          (1,:)  In-plane angle (rad)
   delta          (1,:)  Out-of-plane angle (rad)
   alphaHat       (1,:)  Angle of required force direction

--------------------------------------------------------------------------
	  Reference: Colin R. McInnes, "Solar Sailing: Technology, Dynamics and Mission
              Applications", Springer Praxis, London, 1999, pp. 136-148. 
--------------------------------------------------------------------------

Children:

Common: Database/Constant
Common: Graphics/Plot2D
OrbitMiniToolbox: Support/M2Nu
SC: BasicOrbit/RV2El
SC: BasicOrbit/RVFromKepler