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
dat<-runif.tri(n,Tr)$g
M<-as.numeric(runif.tri(1,Tr)$g) #try also M<-c(1,1,1)
r<-1.4
PEdomtri(dat,Tr,r,M)
IM<-IncMatPEtri(dat,Tr,r,M)
dom.greedy(IM)
dom.exact(IM)
gr.gam<-dom.greedy(IM)
gr.gam
dat[gr.gam$i,]
PEdomtri(rbind(dat,c(5,5)),Tr,r,M)
PEdomtri(dat,Tr,r,M=c(.4,.4))
PEdomtri(rbind(dat,c(5,5)),Tr,r,M=c(.4,.4))
P1<-c(.5,.5)
PEdomtri(P1,Tr,r,M)
Run the code above in your browser using DataLab