This is a great submission! A minor change was required for my application. Reason: The triangulation tends to be slightly over-refined on complex boundaries. For some applications this is a bonus, but for hydrodymanic modelling it is not really necessary. I therefore replaced the line in 'quadtree.m' (about line 64) which reads:
L = 2.0*dist2poly(wm,edgexy,2.0*len);
with the following:
This helps prevent over refinement and a nice mesh. Excellent work Darren!