Quelle  mandelbrot.cob   Sprache: Cobol

       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.