def compute_vc(self, points):
"""compute the Voronoi cell"""
#compute S and V
S = ConvexHull(points).area
V = ConvexHull(points).volume
#voronoi cell
eta = S ** 3 / (36 * np.pi * V ** 2)
return eta