In MATLAB write out a program to integrate the equations of motion of a rigid body. The inertia matrix is given by I = [125 0 0; 0 100 0; 0 0 75] which is a diagonal, where diag operator provides a matrix with given elements placed on its diagonal. Consider three cases where the body rotates 1 rad/sec about each principal axis. Integrate the resulting motion and study the angular rates and the resulting attitude (use any attitude coordinates). For each principal axis case, assume first that a pure spin about the principal axis is performed, and then repeat the simulation where a small 0.1 rad/sec motion is present about another principal axis. Discuss the stability of each motion. The code should produce a total of 6 simulations results when it is ran.

In MATLAB write out a program to integrate the equations of motion of a rigid body. The inertia matrix is given by I = [125 0 0; 0 100 0; 0 0 75] which is a diagonal, where diag operator provides a matrix with given elements placed on its diagonal. Consider three cases where the body rotates 1 rad/sec about each principal axis. Integrate the resulting motion and study the angular rates and the resulting attitude (use any attitude coordinates). For each principal axis case, assume first that a pure spin about the principal axis is performed, and then repeat the simulation where a small 0.1 rad/sec motion is present about another principal axis. Discuss the stability of each motion. The code should produce a total of 6 simulations results when it is ran.