I made a page explaining a bit the Mandelbrot set:
http://iki.fi/warp/Mandelbrot/

  Perhaps I should put a link to it in the povVFAQ?

-- 
char*i="b[7FK@`3NB6>B:b3O6>:B:b3O6><`3:;8:6f733:>::b?7B>:>^B>C73;S1";
main(_,c,m){for(m=32;c=*i++-49;c&m?puts(""):m)for(_=(
c/4)&7;putchar(m),_--?m:(_=(1<<(c&3))-1,(m^=3)&3););}    /*- Warp -*/

Post a reply to this message

"Warp" <war### [at] tagpovrayorg> wrote in message
news:3a6db8fe@news.povray.org...
>   I made a page explaining a bit the Mandelbrot set:
> http://iki.fi/warp/Mandelbrot/
>
>   Perhaps I should put a link to it in the povVFAQ?
>

Nice - now how about the idiots version?

Post a reply to this message

Warp wrote:
> 
>   I made a page explaining a bit the Mandelbrot set:
> http://iki.fi/warp/Mandelbrot/
> 
>   Perhaps I should put a link to it in the povVFAQ?
> 

Looks good, i did not yet have time to read it, but the background color
of the code is disgusting :-)

Christoph

-- 
Christoph Hormann <chr### [at] gmxde>
IsoWood include, radiosity tutorial, TransSkin and other 
things on: http://www.schunter.etc.tu-bs.de/~chris/

Post a reply to this message

Warp wrote:

>   I made a page explaining a bit the Mandelbrot set:
> http://iki.fi/warp/Mandelbrot/

Too late, I gave my extra credit presentation for calculus (on fractals,
specifically julia sets) yesterday. ;)

--
David Fontaine  <dav### [at] faricynet>  ICQ 55354965
My raytracing gallery:  http://davidf.faricy.net/

Post a reply to this message

"Tom Melly" <tom### [at] tomandlucouk> wrote in message
news:3a6dc02b$1@news.povray.org...

Actually, I took some more time last night to read it properly, and it's
more of an idiot's version than I first thought. At first glance, I ran
screaming for the hills when I saw the opening equations.

What I can't find (although I'm not saying it isn't there) is an explanation
of what the imaginary number is (the square-root of -1).

BTW Here is my ascii line-art mandelbrot code for POV - the trouble is that
it doesn't work and I don't understand why not. Can anyone see what's wrong
with it?

#declare Start = -2.00;
#declare End = 2.00;
#declare Iterations = 5;

#declare Imaginary = -Start;
#while(Imaginary >= -End)
  #declare Real = Start;
  #while(Real <= End)
    #declare Inside = true;
    #declare TReal = Real;
    #declare TImaginary = Imaginary;
    #declare Count = 1;
    #while(Count <= Iterations)
      #declare TTReal = TReal;
      #declare TReal = (TReal * TReal) - (TImaginary * TImaginary);
      #declare TImaginary = 2 * TTReal * TImaginary;
      #if((TReal * TReal) + (TImaginary * TImaginary) > 4)
        #declare Inside = false;
        #declare Count = Iterations + 1;
      #else
        #declare Count = Count + 1;
      #end
    #end
    #if(Inside)
      #debug "X"
    #else
      #debug " "
    #end
    #declare Real = Real + ((End - Start)/100);
  #end
  #debug "\n"
  #declare Imaginary = Imaginary - ((End - Start)/50);
#end

Post a reply to this message

Tom Melly <tom### [at] tomandlucouk> wrote:
:       #declare TReal = (TReal * TReal) - (TImaginary * TImaginary);
:       #declare TImaginary = 2 * TTReal * TImaginary;

  You are making Z=Z^2 here, not Z=Z^2+c.

-- 
char*i="b[7FK@`3NB6>B:b3O6>:B:b3O6><`3:;8:6f733:>::b?7B>:>^B>C73;S1";
main(_,c,m){for(m=32;c=*i++-49;c&m?puts(""):m)for(_=(
c/4)&7;putchar(m),_--?m:(_=(1<<(c&3))-1,(m^=3)&3););}    /*- Warp -*/

Post a reply to this message

Christoph Hormann <chr### [at] gmxde> wrote:
: Looks good, i did not yet have time to read it, but the background color
: of the code is disgusting :-)

  Thanks :)

  (Any better suggestion?)

-- 
char*i="b[7FK@`3NB6>B:b3O6>:B:b3O6><`3:;8:6f733:>::b?7B>:>^B>C73;S1";
main(_,c,m){for(m=32;c=*i++-49;c&m?puts(""):m)for(_=(
c/4)&7;putchar(m),_--?m:(_=(1<<(c&3))-1,(m^=3)&3););}    /*- Warp -*/

Post a reply to this message

In article <3a6db8fe@news.povray.org>, war### [at] tagpovrayorg says...
>   I made a page explaining a bit the Mandelbrot set:
> http://iki.fi/warp/Mandelbrot/
> 
>   Perhaps I should put a link to it in the povVFAQ?
> 
I was wondering whether the inner part of what you call the Mandelbrot 
set /is/ part of the set. Isn't the Mandelbrot set just an infinitely 
large set of points forming an infinitely large (continuous?) line; the 
outline as it were? So anything outside or inside this set of points (or 
line of infinite length) would not be part of it?

-- 
Regards,  Sander

Post a reply to this message

"Sander Stols" <san### [at] stolscom> wrote in message
news:MPG### [at] NEWSPOVRAYORG...

> I was wondering whether the inner part of what you call the Mandelbrot
> set /is/ part of the set. Isn't the Mandelbrot set just an infinitely
> large set of points forming an infinitely large (continuous?) line; the
> outline as it were? So anything outside or inside this set of points (or
> line of infinite length) would not be part of it?
>
> --
> Regards,  Sander

IMHO, no. Each point is either inside or outside the set. The boundary has
no special significance in that sense. The only point about the boundary is
that it represents the area of least certainty - you can never be sure if,
with one more iteration, a point might suddenly turn out to be outside the
set.

Post a reply to this message

"Warp" <war### [at] tagpovrayorg> wrote in message
news:3a6ec2ea@news.povray.org...
>
>   You are making Z=Z^2 here, not Z=Z^2+c.
>

Doh! Many thanks - fixed it now.

The final code is:
#declare XStart = -1.5;

#declare XEnd = 5;

#declare YStart = -3/2;

#declare YEnd = 5;

#declare Iterations = 50;

#declare Imaginary = -YStart;

#while(Imaginary >= -YEnd)

    #declare Real = XStart;

    #while(Real <= XEnd)

        #declare Inside = true;

        #declare TReal = Real;

        #declare TImaginary = Imaginary;

        #declare Counter = 1;

        #while(Counter <= Iterations)

            #declare T2Real = TReal;

            #declare TReal = ((TReal * TReal) - (TImaginary * TImaginary));

            #declare TImaginary = (2 * T2Real * TImaginary);

            #if(((TReal * TReal) + (TImaginary * TImaginary)) > 4)

                #declare Inside = false;

                #declare Counter = Iterations;

            #end

            #declare TReal = TReal + Real;

            #declare TImaginary = TImaginary + Imaginary;

            #declare Counter = Counter + 1;

        #end

        #if(Inside)

            #debug "X"

        #else

            #debug " "

        #end

        #declare Real = Real + ((XEnd - XStart)/130);

    #end

    #debug "\n"

    #declare Imaginary = Imaginary - ((YEnd - YStart)/80);

#end

Post a reply to this message