Function spin_matrix 

pub fn spin_matrix(omega: &Vec3) -> Mat3

Spin matrix for angular velocity omega: the time-derivative of a rotation matrix R satisfies Ṙ = omega_hat * R where omega_hat = skew(omega).