import numpy as np
from scipy.spatial import Delaunay

#do triangulation
points = np.random.rand(4, 3)
tri = Delaunay(points)

#indices
tri.simplices

#points
tri.points