Wrote a new function R_2vect() to generate a rotation matrix from 2 vectors.
This rotation matrix is not unique because of axially symmetry allowing the rotation axis to be any
vector within the symmetry plane, so the cross product with the axis-angle convention is used to
define this matrix.
