if (FALSE) {
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)
Xp<-runif.std.tri(n)$gen.points
M<-as.numeric(runif.std.tri(1)$g) #try also M<-c(.6,.2)
r<-1.5
S<-rbind(Xp[1,],Xp[2,])
IndNPETe.domset(S,Xp,r,M)
S<-rbind(Xp[1,],Xp[2,],Xp[3,],Xp[5,],c(.2,.5))
IndNPETe.domset(S,Xp[3,],r,M)
}
Run the code above in your browser using DataLab