buffon.needle = function(n = 100, sec = 0.1, l = 1.0, d = 1.0) { x <- y <- xd <- yd <- phi <- ctr <- numeric(0) sum <- 0; PI <- rep(NA, n) layout(matrix(c(1, 2), 2)) for (j in 1:n) { # Draw parallel lines xlim <- c(-0.5 * l, 1.5 * l); ylim <- c(0, 2 * d) plot(1, xlim = xlim, ylim = ylim, type = 'n', xlab = '', ylab = '', axes = FALSE) abline(h = c(0.5 * d, 1.5 * d), lwd = 2, col = 'blue') # Simula