if (FALSE) {
tau<-1.5
c<-.4
a<-0; b<-10; int<-c(a,b);
#nx is number of X points (target) and ny is number of Y points (nontarget)
nx<-20; ny<-5; #try also nx<-40; ny<-10 or nx<-1000; ny<-10;
set.seed(1)
xf<-(int[2]-int[1])*.1
Xp<-runif(nx,a-xf,b+xf)
Yp<-runif(ny,a,b)
NumArcsCS1D(Xp,Yp,tau,c)
NumArcsCS1D(Xp,Yp,tau,c=.3)
NumArcsCS1D(Xp,Yp,t=2,c)
}
Run the code above in your browser using DataLab