cardtrick.R

secretnum = 1:10
stop = 150

faceval = round(floor(13*runif(stop))+0.999)
faceval = pmin(10,faceval)

data = secretnum
draws = 1:stop
for(i in draws){
  secretnum = secretnum -1
  secretnum = ifelse(secretnum==0, faceval[i],secretnum)
  data=cbind(data, secretnum)
}
draws = c(draws, stop+1)

plot(draws, data[1,], type="l")
for(k in 2:10){
   lines(draws, data[k,], type="l")
}