t = qmt.timeVec(T=10, Ts=0.01)
axis = np.column_stack([np.cos(t), np.zeros_like(t), np.sin(t)])
quat = qmt.quatFromAngleAxis(np.sin(t), axis, debug=True, plot=True)