# \donttest{
c<-.4
a<-0; b<-10; int = c(a,b)
Mc<-centerMc(int,c)
rel.vert.mid.int(6,int,c)
n<-20 #try also n<-40
xr<-range(a,b,Mc)
xf<-(int[2]-int[1])*.5
Xp<-runif(n,a,b)
Rv<-vector()
for (i in 1:n)
Rv<-c(Rv,rel.vert.mid.int(Xp[i],int,c)$rv)
Rv
jit<-.1
yjit<-runif(n,-jit,jit)
Xlim<-range(a,b,Xp)
xd<-Xlim[2]-Xlim[1]
plot(cbind(Mc,0),main="vertex region indices for the points", xlab=" ",
ylab=" ", xlim=Xlim+xd*c(-.05,.05),ylim=3*range(yjit),pch=".",cex=3)
abline(h=0)
points(Xp,yjit)
abline(v=c(a,b,Mc),lty = 2,col=c(1,1,2))
text(Xp,yjit,labels=factor(Rv))
text(cbind(c(a,b,Mc),.02),c("rv=1","rv=2","Mc"))
# }
Run the code above in your browser using DataLab