r = rot4z(theta) rotz produces a 4x4 rotation matrix representing a rotation by theta radians about the z axis. Argument definitions: theta = rotation angle in radians