quaternions available

Callers can supply an optional argument, q which will store the quaternion corresponding to rotation R. From it, the rotation angle and rotation axis can easily be determined.