:ά֮Mandelbrot

ʽ:
[k,x,y]=FractalMandelbrot(f,x0,x1,y0,y1,m,n,Sum,R,Old)
[k,x,y]=FractalMandelbrot(f,x0,x1,y0,y1,m,n,Sum,R)
[k,x,y]=FractalMandelbrot(f,x0,x1,y0,y1,m,n,Sum)
[k,x,y]=FractalMandelbrot(f,x0,x1,y0,y1,m,n)
[k,x,y]=FractalMandelbrot(f,x0,x1,y0,y1,m)
[k,x,y]=FractalMandelbrot(f,x0,x1,y0)
[k,x,y]=FractalMandelbrot(f,x0,x1)
[k,x,y]=FractalMandelbrot(f,x0)
[k,x,y]=FractalMandelbrot(f)

f:ַĺʽ,ıΪx
x0:ʼֵ,ĬΪ-1.5
x1:ֵ,ĬΪ1.5
y0:ʼֵ,Ĭx0һ
y1:ֵ,Ĭx1һ
m:,Ĭ400
n:,Ĭmһ
Sum:,ĬΪ100
R:Mandelbrot뾶,ĬΪ512
Old:ʾԭ,˱Ϊһ1*2ľ.ֱ𴢴ԭĺ,ĬϺԭΪ0

k:m*nľ,ʾյ
x:m*nľ,ʾոƽֵ
y:m*nľ,ʾոƽֵ

ԭ:
Old=[x_0,y_0]ÿһx,yµ
    z.real = x
    z.imag = y
    k = 1
    ʱ( (z.real-x_0)^2+(z.imag-y_0)^2)^0.5<R  And  k<Sum  )µ
    {
       z = f(z)
       k = k+1
    }
    x = z.real
    y = z.imag

ע:{FractalMandelbrot2<\FractalMandelbrot2>}ͬǱصĵһǾ

PS:ʱΪ˻ͼ,Խ z=k-log(log(sqrt(x^2+y^2)))/Log(2) Ϊһֶ

ο:{Mandelbrot<http://linas.org/art-gallery/escape/escape.html>}
:

f="x^2+0.3";
[a,x0,y0]=FractalMandelbrot(f,-1.5,1.5,-1.5,1.5,400,400,100,512);
x={MeshGrid2<\MeshGrid2>}(1,400,400,400);
y=x';
{Cloud<\Cloud>}(x,y,a);

//ִ5ɵõµķͼ

{<http://img2.ph.126.net/OKB6CGrIApIwaIRBKuTRwQ==/6630191454932527043.png>}

//ѵյ뾶Ϊ,ִ2õ½
z=a-log(log(sqrt(x0^2+y0^2)))/Log(2);
cloud(x,y,z)
{<http://img0.ph.126.net/Iqhk6sxyMk2FZqMYigiQmg==/6630499318188335778.png>}