[R] Animation of Mandelbrot Set

[Romain Francois]

Nice !!

Le 04.10.2005 21:04, Tuszynski, Jaroslaw W. a écrit :

>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 at saic.com                     `   \
>
>______________________________________________
>R-help at 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
>  
>
-- 
visit the R Graph Gallery : http://addictedtor.free.fr/graphiques
~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~
~~~~~~      Romain FRANCOIS - http://addictedtor.free.fr         ~~~~~~
~~~~        Etudiant  ISUP - CS3 - Industrie et Services           ~~~~
~~                http://www.isup.cicrp.jussieu.fr/                  ~~
~~~~           Stagiaire INRIA Futurs - Equipe SELECT              ~~~~
~~~~~~   http://www.inria.fr/recherche/equipes/select.fr.html    ~~~~~~
~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~