t<-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<-inci.matCS1D(Xp,Yp,t,c)
IM
dom.num.greedy(IM)
# \donttest{
dom.num.exact(IM) #might take a long time depending on nx
# }
Idom.num.up.bnd(IM,5)
Arcs<-arcsCS1D(Xp,Yp,t,c)
Arcs
summary(Arcs)
plot(Arcs)
inci.matCS1D(Xp,Yp+10,t,c)
t<-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)
inci.matCS1D(Xp,Yp,t,c)
Run the code above in your browser using DataLab