# R code: Game of life on random 2d voronoi tesselation library(deldir) set.seed(1) n = 1000 pts = data.frame(x = runif(n), y = runif(n)) v = deldir(pts$x, pts$y) tiles = tile.list(v) edges = v$dirsgs[, c("ind1", "ind2")] neighbors = split(c(edges$ind2, edges$ind1), c(edges$ind1, edges$ind2)) state = rbinom(n, 1, 0.05) par(mar=rep(0,4)) while(1) { nb = sapply(1:n, function(i) sum(state[neighbors[[as.character(i)]]], na.rm = TRUE)) state = as.numeric(state + nb == 2) plot(0, 0, type = "n", xlim = c(0,1), ylim = c(0,1)) for (i in 1:n) { polygon(tiles[[i]]$x, tiles[[i]]$y, col = grey(1-state[i])) } Sys.sleep(.5) }