rotate
expression function
Returns a 4×4 rotation matrix from an axis and angle.