import numpy as np
from scipy.spatial import Delaunay
import mathplotlib.pyplot as plt
points = np.array([
[2, 4],
[3, 4],
[3, 0],
[2, 2],
[4, 1],
])
simplices = Delaunay(points).simplices
plt.triplot(points[:, 0], points [:, 1], simplices)
plt.scatter(points[:, 0], points [:, 1], color = 'r')
plt.show()