CL_rot_quat2eul — Quaternion to Euler '123' angles
[phi,theta,psi] = CL_rot_quat2eul(q)
NOTE: One must be aware of singularities in the Euler angle parametrization when the pitch (theta) approaches +-%pi or +-%pi/2. These cases must be handled specifically: results may not be accurate!