r<-2
c<-.4
a<-0; b<-10;
nx<-10; ny<-4
set.seed(1)
Xp<-runif(nx,a,b)
Yp<-runif(ny,a,b)
IM<-IncMatPE1D(Xp,Yp,r,c)
IM
dom.greedy(IM)
IndUBdom(IM,6)
dom.exact(IM)
Arcs<-ArcsPEMI(Xp,Yp,r,c)
Arcs
summary(Arcs)
plot(Arcs)
IncMatPE1D(Xp,Yp+10,r,c)
r<-2
c<-.4
a<-0; b<-10;
#nx is number of X points (target) and ny is number of Y points (nontarget)
nx<-20; ny<-4; #try also nx<-40; ny<-10 or nx<-1000; ny<-10;
Xp<-runif(nx,a,b)
Yp<-runif(ny,a,b)
IncMatPE1D(Xp,Yp,r,c)
Run the code above in your browser using DataLab