ジュリア集合はマンデルブロのzとCを交換したもの

• こちらからの続き
• Cをこていしたときに発散しないzの分布

library(fields) # for tim.colors
library(caTools) # for write.gif
munit<-20 # 解像度指定 (30は重い、20くらいが適当)
m = munit^2 # grid size
#Z = complex( real=rep(seq(-1.8,0.6, length.out=m), each=m ), imag=rep(seq(-1.2,1.2, length.out=m), m ) )
Z = complex( real=rep(seq(-5,5, length.out=m), each=m ), imag=rep(seq(-5,5, length.out=m), m ) )
Z = matrix(Z,m,m)

C =complex(r=0.5,i=0.5) # z(n+1)=z(n)^2+Cのz
K<-20 # 漸化式の「漸化」回数(gifアニメの枚数)
X = array(0, c(m,m,K))
for (k in 1:K) {
#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))
}

image(X[,,K], col=tim.colors(256)) # show final image in R


write.gif(X,"Mandelbrot.gif",col=tim.colors(256), delay=100)