Arbitrary Square Bounded Voronoi Diagram

Compute individual Voronoi cell area of 2D point sets bounded in an arbitrary square

This function compute the individual Voronoi cell area of point sets bounded in an arbitrary square.

Inputs:
x : M x 1 array of x-coordinates
y : M x 1 array of y-coordinates
toggleplot : 1 to turn on figures, 0 to turn off figures;
[x1,x2,y1,y2] : 1 x 4 array of square edges

Outputs:
CellArea : M x 1 array of Voronoi cell area bounded in an arbitrary square

Usage:
CellArea=SquareBV(x,y,toggleplot,[x1,x2,y1,y2])
CellArea=SquareBV(x,y,toggleplot)