Quaternion on Mathematica

Import

<< Quaternions`

Multiplication

q1 = Quaternion[1, 2, 3, 4];
q2 = Quaternion[5, 6, 7, 8];
q1 ** q2