function?vorarea
%UNtitleD?Summary?of?this?function?goes?here
%多邊形面積
%???Detailed?explanation?goes?here
clear;clc;
m=5;

x=[4228325869];
y=[8423433769];
%??[vc]=voronoin（[x;y]‘）???%?不設置邊界時的情況，v是所有的voronoi邊的頂點坐標，c{i}就是第i個點周圍的voronoi邊的頂點下標
pbound=Polyhedron（[0?0;0?100;100?100;100?0]）;
VV=mpt_voronoi（[x;?y]‘bound‘pbound）;
for?iPoly?=?1:m
????%?make?points?in?order
????k?=?convhull（VV.Set（iPoly）.V（:1）VV.Set（iPoly）.V（:2））;
????%?make?clockwise
????[pXpY]?=?poly2cw（VV.Set（iPoly）.V（k1）VV.Set（iPoly）.V（k2））;
???