Lloyd's Algorithm
Runs Lloyd's algorithm on the particles at xy positions (Px,Py) within the boundary polygon crs for numIterations iterations.
showPlot = true will display the results graphically.
Lloyd's algorithm starts with an initial distribution of samples or points and consists of repeatedly executing one relaxation step:
1. The Voronoi diagram of all the points is computed.
2. Each cell of the Voronoi diagram is integrated and the centroid is computed.
3. Each point is then moved to the centroid of its Voronoi cell.
Run with no input to see example.
Uses http://www.mathworks.com/matlabcentral/fileexchange/34428-voronoilimit,
which requires the Polybool function of the mapping toolbox to run.