Last-modified: 2019-01-27 (日) 01:45:41 (26d)

---

---

†

†

i * i = -1

†

†

†

†

†

globals [ xmin ymin xmax ymax  itr colscale ]

†

to minit 
	set itr 100
	set xmin (cenx - (rsize * 0.5)) / scale
  	set xmax (cenx + (rsize * 0.5)) / scale
  	set ymin (ceny - (rsize * 0.5)) / scale
  	set ymax (ceny + (rsize * 0.5)) / scale
end

†

to mcalc

  ask-patches [

    let [:cx (xmin + xcor * (xmax - xmin)) ]
    let [:cy (ymin + ycor * (ymax - ymin)) ]

    let [:zx 0]
    let [:zy 0]
    let [:tx :cx]
    let [:ty :cy]
    let [:cnt 0]

    repeat itr [
      ;; z = z*z + c iteration
      set :tx (:zx * :zx) - (:zy * :zy)
      set :ty (:zx * :zy) + (:zy * :zx)
      set :zx (:zx + :tx + :cx)
      set :zy (:zy + :ty + :cy)

      let [ :a sqrt ((:zx * :zx) + (:zy * :zy))]

      if :a > 4.0 [ setpc  :cnt stop ]
      set :cnt (:cnt + 1)
    ]
    setpc black
 ]
end

†

to mandelbrot
  ca
  minit
  mcalc
end

GUI †

†

†

---

†

2005/06/17 †

patches-own [ xmin ymin xmax ymax  itr colscale ]

globals [ xmin ymin xmax ymax  itr colscale ]

mcalc

---