.^ Array power
(Quaternion overloading of standard MATLAB® function)
Y = X .^ Y
X .^ Y raises the quaternion X to the quaternion power Y, elementwise. It can also handle cases where one or both parameters is a scalar (a single quaternion).
The general case is handled using logarithms, but certain special cases of the right parameter are computed using more direct formulae for greater speed and accuracy.
>> randv .^ 2 ans = -1 + 0 * I + 0 * J + 0 * K