use builtin circumcenters from Delaunay triangulation
remove superfluous custom circumcenter algorithm
made Z argument for voronoi function optional
documentation and better integration with axes
pylab demo script