code: def unit_poly_verts(theta):
      
      x0, y0, r = [0.5] * 3
      verts = [(r*np.cos(t) + x0, r*np.sin(t) + y0) for t in theta]
      return verts
docs: Return vertices of polygon for subplot axes.
  
      This polygon is circumscribed by a unit circle centered at (0.5, 0.5)