ReturnMatrix rotd(const Real theta, const ColumnVector & k1, const ColumnVector & k2);
Note: the column vectors k1 and k2 must have a length of at least 3. Only the first 3 elements are used.
Matrix