# NOT RUN {
A<-c(1,1); B<-c(2,0); C<-c(1.5,2);
Tr<-rbind(A,B,C);
n<-10
set.seed(1)
dat<-runif.tri(n,Tr)$gen.points
M<-as.numeric(runif.tri(1,Tr)$g) #try also M<-c(1.6,1.0)
t<-.5
IndCSdomUBtri(dat,1,t,Tr,M)
for (k in 1:n)
print(c(k,IndCSdomUBtri(dat,k,t,Tr,M)))
IndCSdomUBtri(dat,k=4,t,Tr,M)
P<-c(.4,.2)
IndCSdomUBtri(P,1,t,Tr,M)
IndCSdomUBtri(rbind(dat,dat),3,t,Tr,M)
dat.fr<-data.frame(a=dat)
IndCSdomUBtri(dat.fr,1,t,Tr,M)
# }
# NOT RUN {
# }
Run the code above in your browser using DataLab