How to find the area of voronoi shapes

clc
clear all
close all  
n= 60;
  x=720*rand(1,n);
  y=360*rand(1,n);
  voronoi(x,y)
  [vx,vy] =voronoi(x,y);
  plot(x,y,'r+',vx,vy,'b-'); 
  axis equal;
With this example of a voronoi diagram, is it possible to find the areas of each respective voronoi shape?  If so, how?

[Sean]

figure(2)
for j=1:length(vx(1,:))
   line([vx(1,j) vx(2,j)],[vy(1,j) vy(2,j)])
  end
  axis equal

Outlines the voronoi diagram.