# NOT RUN {
P <- runifpoint(15)
X <- dirichlet(P)
plot(X, col="red")
B <- bdist.tiles(X)
# identify tiles that do not touch the boundary
plot(X[B > 0], add=TRUE, col="green", lwd=3)
# }
Run the code above in your browser using DataCamp Workspace