tau<-2
c<-.4
a<-0; b<-10; int<-c(a,b);
#n is number of X points
n<-10; #try also n<-20
xf<-(int[2]-int[1])*.1
set.seed(1)
Xp<-runif(n,a-xf,b+xf)
Arcs<-arcsCSint(Xp,int,tau,c)
Arcs
summary(Arcs)
plot(Arcs)
Xp<-runif(n,a+10,b+10)
Arcs=arcsCSint(Xp,int,tau,c)
Arcs
summary(Arcs)
plot(Arcs)
Run the code above in your browser using DataLab