RErange=float(r1-r0)
IMrange=float(i1-i0)
maxitter=1000
maxZ=2000
Nre = N
Nim = Nre*IMrange/RErange
temp = intarr(Nre,Nim)
re_step = RErange/N
im_step = IMrange/N
c = complex(r0, i0)
for re=0, Nre-1 do begin
c = complex(float(c)+re_step, i0)
for im=0, Nim-1 do begin
z=complex(0,0)
c = c + complex(0, im_step)
i=0L
repeat begin
i = i +1
z = z^2 + c
endrep until (abs(z) gt maxZ)
or (i gt maxitter)
temp(re, im)=i
endfor
endfor
return, temp
end