q2dr function

Get derivative of 3D rotation matrix from quaternion