Function File: qn = unit (q)

Normalize quaternion to length 1 (unit quaternion).

          q = w + x*i + y*j + z*k
          unit (q) = q ./ sqrt (w.^2 + x.^2 + y.^2 + z.^2)