Quelle  merkur.cob

       identification division.
       program-id. Merkur.
      * aus Prof. Ulrich Gähdes Vorlesung
      * Begriffsdefinition und Hypothesenbildung
      * am 19. Januar 2011 gegebenen Beispiel
      * über die Merkuranomalie
       author. Jens.
       date-written. 19.1.2011.
       date-compiled.
       data division.
       working-storage section.
         78 g value 9.81.
         78 c value 9.81.
         78 a value 9.81.
      *  -------------------
         77 i pic 9(2).
         77 f1 pic 9(4)V99.
         77 f2 pic 9(4)V99.
         77 diff pic 9(4)V99.
       linkage section.
         77 r pic 9(4)V99.
         77 Masse pic 9(4)V99.
         77 m pic 9(4)V99.
         77 f pic 9(4)V99.
       procedure division using r m returning f.
          perform Newton.
          compute diff=f;
          perform Clairant.
          compute diff=diff-f.
          display "diff Newton-Clairant="diff 
      *  -------------------
          stop run. 
      *  -------------------
      *  für alle r, Masse gilt
      *  -------------------
       Clairant.
         compute f1=function pow(r,-2)
         compute f2=function pow(f1,2)
         compute f=g*Masse*f1+a*f2
      *  -------------------
      *  für alle r gilt
      *  -------------------
       Hall.
         move 1 to f1
         perform varying i from 1 until i>m multiply f1 by r end-perform
         compute f=c*1/f1
      *  -------------------
      *  für alle r,m,Masse gilt
      *  -------------------
       Newton.
         compute f1=function pow(r,-2)
         compute f=g*Masse*m*f1
        end-program Merkur.