function q = quaternion(vector)
%function q = quaternion(vector)
%inputs:  vector  (3X1 array)
%output    q       1X4 vector (pure quaternion)
  q = [0, vector'];
  
  