マンデルブロー集合の描画 Rendering Mandelbrot set

pythonの練習でマンデルブロー集合を描画するコードを書いた。

Fig:実行結果　The result of execution

from scipy import *

from pylab import *

STOP_VALUE = 4

CX_AREA=[-1,0.5]

CY_AREA=[-1,1]

STEP = 0.01

def judge_conv_mandelbrot_seq(c):

    mandelbrot_seq = 0.+0.j

    for n in range(0,50):

        mandelbrot_seq=mandelbrot_seq**2+c

        if (abs(real(mandelbrot_seq))>STOP_VALUE)or(abs(imag(mandelbrot_seq))>STOP_VALUE):

            break

#    if abs(mandelbrot_seq)    if abs(mandelbrot_seq)>STOP_VALUE:

        return 0

    else:

        return 1

    

mandelbrot_set=[]

for cx in arange(CX_AREA[0],CX_AREA[1],STEP):

    for cy in arange(CY_AREA[0],CY_AREA[1],STEP):

        if judge_conv_mandelbrot_seq(complex(cx,cy)):

             mandelbrot_set.append(complex(cx,cy))

        else:

            pass

plot(real(mandelbrot_set),imag(mandelbrot_set),',')



show()