  In this article, a brief description of Discrete Mechanics and Variational
Integrators which preserve the symplectic structure of the flow will be
provided and a Newton-Raphson algorithm that can be used to solve implicit
equations on the SO(3) manifold will be developed. These techniques will be
used to simulate the rotational dynamics of a rigid body evolving on the Lie
Group SO(3).
