A<-c(0,0); B<-c(1,0); C<-c(1/2,sqrt(3)/2);
n<-10 #try also n<-20
set.seed(1)
dat<-runifTe(n)$gen.points
M<-c(.6,.2) #try also M<-c(1,1,1)
NumArcsPETe(dat,r=1.25,M)
NumArcsPETe(dat,r=1.5,M)
NumArcsPETe(dat,r=2,M)
NumArcsPETe(rbind(dat,c(0,1)),r=2,M)
NumArcsPETe(c(.2,.3),r=2,M)
NumArcsPETe(dat,r=1.5,M);
dat.fr<-data.frame(a=dat)
NumArcsPETe(dat.fr,r=1.5,M);
Run the code above in your browser using DataLab