r = rot4y(theta) roty produces a 4x4 rotation matrix representing a rotation by theta radians about the y axis. Argument definitions: theta = rotation angle in radians