I can successfully add and multiply quaternions, but I can't figure out how to raise to the power of a quaternion? Advice?
In my case, this works:
 
Needs["Quaternions`"]
Needs["Quaternions`"]
NumberForm[
 N[(Quaternion[0, 0, 0, 1] + 
       Quaternion[0, 0, 0, -1/(n + t)])^(n + t), 5] 
But this does not work:
 
Needs["Quaternions`"]
NumberForm[
 N[(Quaternion[0, 0, 0, 1] + 
       Quaternion[0, 0, 0, -1/(n + t)])^(Quaternion[0, 0, 0, 1]**(n + t)), 5] 
Suggestions?