Function roty 

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

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

Input:

angle - Angle of rotation in radians