Quelle  fraktal.cob

       identification division.
       program-id. fraktal.
      * Autor Jens-D. Doll
      * T(z) = (a*z+b*z)/(c*z+d*z)
      * unfertig, siehe Bild der Wissenschaft 5-2008 o.ä.
       author. "JD".
       date-written. 4.7.2006.
       date-compiled.
       data division.
       working-storage section.
         77 Eps       pic 9(4)V9(8).
         77 i         pic 9(4).
         77 maxi      pic 9(4).
         77 a         pic 9(4).
         77 b         pic 9(4).
         77 c         pic 9(4).
         77 d         pic 9(4).
         77 z         pic 9(4)V9(4).
         77 zi        pic 9(4)V9(4).
         77 u         pic 9(4)V9(4).
         77 ui        pic 9(4)V9(4).
       procedure division.
         display "start, enter a,b,c,d"
         accept a b c d
         move 0 to z zi i
         perform until d<eps or i>maxi
           compute z=a*z
           compute zi=b*z
           compute z=z+b*z
           compute zi=zi+b*z
           compute u=a*z
           compute ui=b*z
           compute u=u+b*z
           compute ui=ui+b*z
           compute z=(z+zi) / (u+ui)
           add 1 to i   
         end-perform
         if i<maxi then display "done, sequence converges"
         else           display "done, diversion!" .
        end-program fraktal.