docs: def rotation(angle):
      """Rotation about the Z axis (in the XY plane)"""
      return N.array([[N.cos(angle),-N.sin(angle),0],
          [N.sin(angle), N.cos(angle),0],
          [0           , 0           ,1]])