JacksonPollock.R

JacksonPollock <- function(n=5,r1=5,r2=0.5,drip=100,col="black"){
  old<-par(no.readonly = TRUE)
  par(mai=c(0,0,0,0))
  plot(c(0,1),c(0,1),type="n")
  pos<-locator(n)
  points(pos$x,pos$y,cex=rexp(n,1/r1),pch=16,col=col)
  d <-rpois(1,drip)
  u <-runif(d)
  for(i in 2:n){
  points(u*pos$x[i-1]+(1-u)*pos$x[i],u*pos$y[i-1]+(1-u)*pos$y[i],
         cex=rexp(d,1/r2),pch=16,col=col)
  }
  par(old)
}