Function File: theta = arg (q)

Compute the argument or phase of quaternion q in radians. theta is defined as atan2 (sqrt (q.x.^2 + q.y.^2 + q.z.^2), q.w). The argument theta lies in the range (0, pi).