From: Tuszynski, Jaroslaw W. <JAROSLAW.W.TUSZYNSKI_at_saic.com>

Date: Wed 05 Oct 2005 - 05:04:57 EST

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)) # show final image in R write.gif(X, "Mandelbrot.gif", col=tim.colors(256), delay=100) # drop "Mandelbrot.gif" file from current directory on any web brouser to see the animation

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Oct 05 05:06:34 2005

Date: Wed 05 Oct 2005 - 05:04:57 EST

Hi,

I was playing with Mandelbrot sets and come up with the following code, I thought I would share:

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)) # show final image in R write.gif(X, "Mandelbrot.gif", col=tim.colors(256), delay=100) # drop "Mandelbrot.gif" file from current directory on any web brouser to see the animation

Jarek

====================================================\==== Jarek Tuszynski, PhD. o / \ Science Applications International Corporation <\__,| (703) 676-4192 "> \ Jaroslaw.W.Tuszynski@saic.com ` \ ______________________________________________R-help@stat.math.ethz.ch mailing list

https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html Received on Wed Oct 05 05:06:34 2005

*
This archive was generated by hypermail 2.1.8
: Fri 03 Mar 2006 - 03:40:35 EST
*