# R program to display the estimated tau values from "canes.c".
# Should be run *AFTER* running "canes.c".

taufilein = matrix( scan("tauvals"), ncol=3, byrow=T )

x = taufilein[,1]
y = taufilein[,2]
tau = taufilein[,3]
K = 30

fact = 10 / K  # to make total display be 10 seconds
numpersec = 10

plot(x, y, pch="-")

for (i in 1:(K*numpersec)) {
  Sys.sleep(fact/numpersec)
  points(x[tau<=i/numpersec], y[tau<=i/numpersec], col="red", pch="*")
}

# hist(tau, breaks=100)