Quelle  matrixmult.cob

       identification division.
       program-id. matrixmult.
      *  aus Wirth, systematisches Programmieren p. 85 
       author. "JD".
       date-written. 11.9.2005.
       date-compiled.
       data division.
       working-storage section.
       78 n   pic 9(4) value 10.
       78 m   pic 9(4) value 10.
       78 p   pic 9(4) value 10.
       77 s   pic 9(4)V9(4).
       77 i   pic 9(4).
       77 j   pic 9(4).
       77 k   pic 9(4).
       linkage section.
       01 fa occurs m.
         05 a pic 9(4)V9(4) occurs p.
       01 fb occurs p.
         05 b pic 9(4)V9(4) occurs n.
       01 fc occurs m.
         05 c pic 9(4)V9(4) occurs n.
       procedure division using fa fb fc.
      * multiply matrices c = a * b
         perform varying i from 1 by 1 until i >= n
           perform varying j from 1 by 1 until j >= n
             move 0 to s
             perform varying k from 1 by 1 until k >= p
               compute s = s + a(i,k) * b(k,j)
             end-perform
             move s to c(i,j)
           end-perform
         end-perform.
      * now c * a * b
        end-program matrixmult.