Vector or cross product of two pure quaternions
Y = vector_product(A, B)
vector_product(A, B) returns the vector (or cross) product of the two arguments, elementwise on the two arrays, which must be pure quaternions.
>> vector_product(qi, qj) == qk ans = 1