X = [ 0.5    0
      0      0.5
     -0.5   -0.5
     -0.2   -0.1
     -0.1    0.1
      0.1   -0.1
      0.1    0.1 ];
dt = DelaunayTri(X);
[V,R] = voronoiDiagram(dt);
for i=1:length(R)
if all(R{i}~=1)   
patch(V(R{i},1),V(R{i},2),i); % use color i.
end
end