AttitudeProfileToTorque:

--------------------------------------------------------------------------
   Compute an equivalent torque sequence from an attitude profile.
   The torque is calculated using an equivalent single-axis rotation for each
   step in the profile.
   This function has a built-in demo for LVLH tracking in an elliptical orbit.

   Since version 7.
--------------------------------------------------------------------------
   Form:
   [torque,omega] = AttitudeProfileToTorque( q, t, inertia, w0 )
--------------------------------------------------------------------------

   -------
   Inputs
   -------
   q        (4,n)   Quaternion profile
   t        (1,n)   Time vector
   inertia  (3,3)   Vehicle inertia
   w0       (3,1)   Initial rotational velocity (default is zeros)

   -------
   Outputs
   -------
   torque   (3,n-1)   Torque sequence
   omega    (3,n)   Angular velocity

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