# NOT RUN {
A<-c(0,0); B<-c(1,0); C<-c(1/2,sqrt(3)/2);
Te<-rbind(A,B,C);
n<-10
set.seed(1)
dat<-runifTe(n)$gen.points
M<-as.numeric(runifTe(1)$g) #try also M<-c(.6,.2)
t<-.5
IndCSdomUBTe(dat,1,t,M)
for (k in 1:n)
print(c(k,IndCSdomUBTe(dat,k,t,M)$IndUBdom))
print(c(k,IndCSdomUBTe(dat,k,t,M)$domUB))
IndCSdomUBTe(dat,k=4,t,M)
P<-c(.4,.2)
IndCSdomUBTe(P,1,t,M)
IndCSdomUBTe(rbind(dat,dat),3,t,M)
dat.fr<-data.frame(a=dat)
IndCSdomUBTe(dat.fr,1,t,M)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab