Skip to content

fractals

3 messages · Erin Hodgess, Martin Maechler, robin hankin

Original
Erin Hodgess
# 
Dear R People:

Does anyone have any code for Fractals, chaos,
or anything like that, please?

This is strictly for demo purposes...decorative only.

This is R version 1.5.1 for Windows.


Thank you in advance!

Have a great weekend!

Sincerely,
Erin Hodgess
Associate Professor
Department of Computer and Mathematical Sciences
University of Houston - Downtown
1 Main Street 
Houston, TX 77002
mailto: hodgess at uhddx01.dt.uh.edu
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
1 day later
Martin Maechler
# 
Erin> Dear R People:
    Erin> Does anyone have any code for Fractals, chaos,
    Erin> or anything like that, please?

    Erin> This is strictly for demo purposes...decorative only.

On a vacation evening, I've produced quite nice Mandelbrot
sets (+ coloring schemes) -- all in R -- just to see if it was
fast enough without C nor special bitmap displays.  

I had been pleasantly surprised how fast and beatiful the
results were. Feedback (to me privately) is welcome.

Martin Maechler <maechler at stat.math.ethz.ch>	http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum  LEO C16	Leonhardstr. 27
ETH (Federal Inst. Technology)	8092 Zurich	SWITZERLAND
phone: x-41-1-632-3408		fax: ...-1228			<><

-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: Mandelbrot.R
Url: https://stat.ethz.ch/pipermail/r-help/attachments/20020812/6e0bbe7f/Mandelbrot.pl
robin hankin
# 
Hi everyone.  

Try:

R>  image(sylvester.matrix(9,pattern=cbind(c(1,0),c(1,1))))

where 

R> sylvester.matrix <- function (n, pattern = rbind(c(1, 1), c(1, -1))) 
{
    n <- as.integer(n)
    stopifnot(n >= 0)
    if (n == 0) {
        return(as.matrix(1))
    }
    else {
        return(kronecker(pattern, sylvester.matrix(n - 1, pattern)))
    }
}