QToConeClock:

--------------------------------------------------------------------------
   Computes cone and clock angles from an inertial to body quaternion.
   The cone and clock system is defined with one axis along the sun (x)
   vector, one along the cross product of the orbit normal and sun vector (y)
   and the third vector completes the orthogonal set (z). The clock angle
   is measured from the y axis.

   Since version 7.
--------------------------------------------------------------------------
   Form:
   [cone, clock] = QToConeClock( qIToB, r, v, s )
--------------------------------------------------------------------------

   ------
   Inputs
   ------
   qIToB         (4,n) Quaternions
   r             (3,:) Orbit vector
   v             (3,:) Velocity vector
   s             (3,:) Sun unit vector (from sc to sun)

   -------
   Outputs
   -------
   cone          (1,:) Cone angle (rad)
   clock         (1,:) Clock angle (rad)

--------------------------------------------------------------------------
	Reference: Dachwald, B., "Low Thrust Trajectory Optimization and 
              Interplanetary Mission Analysis Using 
              Evolutionary Neurocontrol," PhD Thesis
              Universitat Bundeswehr Munchen, FRG, April 2004, p. 8.
--------------------------------------------------------------------------

Children:

Common: Quaternion/Q2Mat
Math: Linear/Cross
Math: Linear/Unit