# \donttest{
a<-0; b<-10
c<-.4
r<-2
#nx is number of X points (target) and ny is number of Y points (nontarget)
nx<-15; ny<-4; #try also nx<-40; ny<-10 or nx<-1000; ny<-10;
set.seed(1)
Xp<-runif(nx,a,b)
Yp<-runif(ny,a,b)
PEdom.num1D(Xp,Yp,r,c)
PEdom.num1D(Xp,Yp,r,c=.25)
PEdom.num1D(Xp,Yp,r=1.25,c)
# }
Run the code above in your browser using DataLab