Si bien no hace mucho publicaba una entrada sobre el triángulo de Sierpinsky, mi tocayo Carlos Ortega (y ahora gentil colaborador) nos ha proporcionado un enlace en este blog a un pedazo de código que bien vale la pena replicar aquí para el solaz (y tal vez, incluso, provecho) de los lectores de estas páginas. Es:
library(fields) # for tim.colors
library(caTools) # for write.gif
m = 400 # grid size
C <- complex(
real=rep(seq(-1.8,0.6, length.out=m), each=m ),
imag=rep(seq(-1.2,1.2, length.out=m), m ) )
C <- matrix(C,m,m)
Z <- 0
X <- array(0, c(m,m,20))
for (k in 1:20) {
Z <- Z^2+C
X[,,k] <- exp(-abs(Z))
}
image(X[,,k], col=tim.colors(256))
write.gif(X, "Mandelbrot.gif", col=tim.colors(256), delay=100)(extraído de aquí).
El resultado, infinitamente mejor que la tele:
