帮我实现一个Delaunay算法

Delaunay算法是一种用于生成凸包的算法，主要用于计算离散点集的三角剖分。你需要输入离散点集，然后使用Delaunay算法生成凸包。Delaunay三角剖分是一个生成的三角形网格，它满足以下两个条件：1.所有的外接圆内不包含离散点集中的任何点。2.将每个三角形的外接圆放在这个三角形的所有点集中都是唯一的。

下面是一个Delaunay算法的代码示例：

def delaunay(points):
"""
Compute the Delaunay triangulation of a set of points.
"""
tri = Delaunay(points)

# Plot triangulation
plt.triplot(points[:,0], points[:,1], tri.simplices.copy())
plt.plot(points[:,0], points[:,1], 'o')
plt.show()