# Wiener process approximated with a random walk
n <- 1000
s <- 0.1
plot((0:n)/n, cumsum(c(0, 2 * s * rbinom(n,1,0.5) - s)),
  type="l", main="Wiener Process", xlab="time", ylab="distance")