load "tetgen" 
int nn=100; 
real[int] xx(nn),yy(nn),zz(nn);

for(int i=0; i<nn; ++i)
{
  xx[i] = randreal1();
  yy[i] = randreal1();
  zz[i] = randreal1();
}
mesh3 Th=tetgconvexhull(xx,yy,zz,region=1,label=3);
plot(Th);
