# \donttest{
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)
r<-1.5
n<-3 #try also n<-20
Xp<-runif.std.tetra(n)$g #try also Xp[,1]<-Xp[,1]+1
plotPEregs.std.tetra(Xp[1:3,],r)
P1<-c(.1,.1,.1)
plotPEregs.std.tetra(rbind(P1,P1),r)
# }
Run the code above in your browser using DataLab