pts = np.random.rand(10, 3)
dt = Delaunay(pts)
tets = dt.points[dt.simplices]
vol = np.sum(tetrahedron_volume(tets[:, (0)], tets[:, (1)], tets[:, (2)],
    tets[:, (3)]))
Volume of convex hull with QHull from SciPy