This is a 2d Delaunay mesh generator that uses the advancing front technique to tesselate the domain with triangles. For anisotropic metric, the Delaunay criterion is modified to the elements' Steiner ellipses not containing any other points.
Plotting: The PLPlot library is used for plotting: http://plplot.sourceforge.net