Re: [R] Animation of Mandelbrot Set

From: Romain Francois <francoisromain_at_free.fr>
Date: Wed 05 Oct 2005 - 05:43:59 EST

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@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
>
>

-- 
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    ~~~~~~
~~~~~~~~ ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~

______________________________________________
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:59:55 2005

This archive was generated by hypermail 2.1.8 : Sun 23 Oct 2005 - 18:18:10 EST