|
identification division.
program-id. mandelbrot.
* Autor Jens-D. Doll
* Dies ist die Mandelbrot-Folge
* Zn+1 = Zn**2 + C
* in der komplexen Ebenee
* Es wird jeweils die Konvergenz
* der Folge nach 100 Schritten untersucht
author. "JD".
date-written. 27.8.2006.
date-compiled.
data division.
working-storage section.
77 Conv pic 9(1).
77 Cre pic 9(4)V9(4).
77 Cim pic 9(4)V9(4).
77 Zre pic 9(4)V9(4).
77 Zim pic 9(4)V9(4).
77 n pic 9(4).
77 nmax pic 9(4) value 100.
procedure division.
display "start"
perform varying Cre from -1 by 0.01 until Cre >= 1
perform varying Cim from -1 by 0.01 until Cim >= 1
perform varying n from 1 by 1 until n > nmax
compute Zre = (Zre*Zre - Zim*Zim) + Cre
compute Zim = Zim + Cim
end-perform
if Cre*Cre+Cim*Cim>1 then move 0 to Conv
else move 1 to Conv end-if
call plot using Cre Cim Conv
end-perform
end-perform
display "done".
end-program mandelbrot.
¤ Dauer der Verarbeitung: 0.15 Sekunden
(vorverarbeitet)
¤
|
Laden
Fehler beim Verzeichnis:
in der Quellcodebibliothek suchen
Die farbliche Syntaxdarstellung ist noch experimentell.
|
