Function rotz 

pub fn rotz<T: Float>(angle: T) -> M33<T>

Compute the rotation around the z axis(in cartesian coordinates)

Input:

angle - Angle of rotation in radians