Rigid Bidy Rotation
Quaternion
rotation matrix R
zeta-axis
matrix H
constraint stabilization
external torque
dynamic equation of rotation using quaternion
solve the dynamic equation of rotation using quaternion