docs: def rotation_matrix_2D(th):
      '''
      rotation_matrix_2D(th) yields a 2D numpy rotation matrix th radians counter-clockwise about the
      origin.
      '''
      s = np.sin(th)
      c = np.cos(th)
      return np.array([[c, -s], [s, c]])