pub fn rotation_2d(radians: f64) -> Array2<f64>
generates a matrix describing 2d rotation around origin