- こちらからの続き
- 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)