# \donttest{
A<-c(0,0); B<-c(1,0); C<-c(1/2,sqrt(3)/2)
Tr<-rbind(A,B,C)
n<-10 #try also n<-20
Xp<-runif.tri(n,Tr)$g
M<-as.numeric(runif.tri(1,Tr)$g) #try also M<-c(1,1,1)
r<-1.4
PEdom.num.tri(Xp,Tr,r,M)
IM<-inci.matPEtri(Xp,Tr,r,M)
dom.num.greedy #try also dom.num.exact(IM)
gr.gam<-dom.num.greedy(IM)
gr.gam
Xp[gr.gam$i,]
PEdom.num.tri(Xp,Tr,r,M=c(.4,.4))
# }
Run the code above in your browser using DataLab