docs: def polygon_diameter(points):
      ''' Compute the maximun euclidian distance between any two points
      in a list of points
      '''
      return max(point_dist(p0, p1) for (p0, p1) in combinations(points, 2))