# NOT RUN {
A<-c(0,0,0); B<-c(1,0,0); C<-c(1/2,sqrt(3)/2,0); D<-c(1/2,sqrt(3)/6,sqrt(6)/3)
tetra<-rbind(A,B,C,D)
n<-10 #try also n<-20
dat<-runif.tetra(n,tetra)$g
M<-"CM" #try also M<-"CC"
r<-1.25
PEdom.tetra(dat,tetra,r,M)
PEdom.tetra(rbind(dat,c(5,5,5)),tetra,r,M)
P1<-c(.5,.5,.5)
PEdom.tetra(P1,tetra,r,M)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab