delaunay function

Delaunay triangulation for points