# \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)
n<-10 #try also n<-20
set.seed(1)
Xp<-runif.tetra(n,tetra)$g
M<-"CM" #try also M<-"CC"
r<-1.25
Narcs = num.arcsPEtetra(Xp,tetra,r,M)
Narcs
summary(Narcs)
#plot(Narcs)
# }
Run the code above in your browser using DataLab