# 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.2)
IndASdomUBtri(dat,1,Tr)
for (k in 1:n)
print(c(k,IndASdomUBtri(dat,k,Tr,M)))
IndASdomUBtri(dat,k=4,Tr,M)
P<-c(.4,.2)
IndASdomUBtri(P,1,Tr,M)
IndASdomUBtri(rbind(dat,dat),k=2,Tr,M)
dat.fr<-data.frame(a=dat)
IndASdomUBtri(dat.fr,1,Tr,M)
# }
Run the code above in your browser using DataLab