Execute the code on http://www.statisticsblog.com/2012/12/let-it-snow/ directly

Let It Snow.R

# Code by Matt Asher for statisticsblog.com

# Feel free to modify and redistribute

# How many flakes do you want to fall?

flakes = 100

# Width and height of your space

width = 800

height = 600

# Initial wind

wind = 0

# Setup the background of the plot and margins

par(bg = "black")

par(oma=c(0,0,0,0))

par(mar=c(0,0,0,0))

plot(0, 0, col="black", pch=".", xlim=c(0,width), ylim=c(0,height), axes=F)

for(i in 1:flakes) {

startY = height

startX = runif(1,1,width)

xPos = startX

yPos = startY

for(j in 1:height) {

# Optional drift in wind

wind = wind + rcauchy(1,0,.05)

# Update snowflake position

xPos = xPos + rnorm(1,.1,1.5)

yPos = yPos - runif(1,4,20)

# Are we in the space, if so display it

if(xPos > 0 & xPos <= width & yPos > 0 & yPos <= height) {

points(round(xPos), round(yPos), col="white", pch=8)

# System dely, slows down the flakes

Sys.sleep(0.1)

}

}

}