IDL-Funktion zur Berechnung von Julia-Mengen

FUNCTION julia, N, c=c
maxitter=1000
maxZ=10000
if not keyword_set(c) then c = complex(-0.194, 0.6557)
temp = intarr(N, N)
for re=0, N-1 do $
    for im=0, N-1 do begin
      z=dcomplex((re-N/2)*2.*RErange/N, (im-N/2)*2*IMrange/N)
      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
return, temp
end