Added qmult for quaternion multiplication. Tested for vector operation. Need examples in docstring.
Fix qmult for N (more than 1) inputs.
Fix input_check_Nx1 and input_check_Nx3, flatten output to 1-D array when able