Watch this Author's files
I found a partial solution to my enquiry, voronoin has to be called with the option 'QJ'
This file is great, however there seem to be a problem with colinear points. I got an error message while testing for :
x = [0.25; 0.75; 0.25; 0.75];
y = [0.25; 0.25; 0.75; 0.75];