experience and innovation

-------------------------------------------------------------------------------
   Simulate the linear and angular states of a rigid body using "AC.m". 
   Four different geometries are available to choose from. 
   A rotor may be included or omitted.

   The motion is simulated using two different cases. 
     case 0   --    The origin of the body frame coincides with the CG
     case 1   --    The origin of the body frame is translated to a non-CG point

   The inertias for case 0 and case 1 are different, according to the parallel
   axis theorem. 

   The initial angular velocity is different for each case IF the inertia matrix 
   has non-zero off-diagonal terms.

   No aerodynamic forces or torques are considered.

   VERIFICATION
   ============

   * The trajectories of both the CG and a non-CG point on the body are computed 
     for each case, and the errors are plotted. 

   * The total angular momentum (body + rotors) is plotted. It should be constant.

-------------------------------------------------------------------------------
   Form:
   RigidBodyDemo( geometry, rotor )
-------------------------------------------------------------------------------

   ------
   Inputs
   ------
   geometry       (1)     Number 1,2,3,4 -- indicating the type of geometry to use
                             1 - Hemisphere                (inertia is pure diagonal)
                             2 - Cylinder                  (inertia is pure diagonal)
                             3 - Rectangular Tetrahedron   (inertia has off-diagonal terms)
                             4 - Half-Cone                 (inertia has off-diagonal terms)

   rotor          (1)     Flag 0,1 -- indicating whether to include a rotor or not

   -------
   Outputs
   -------
   none

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