identification division.
       program-id. pow-fast.
      * aus Dirk Hoffmann,
      * Software-Qualität
      * p339, Springer 2008
       author. "JD".
       date-written. 28.1.2009.
       date-compiled.
       data division.
       working-storage section.
        77 k pic 9(4).
        77 p pic 9(16).
        linkage section.
        77 a   pic 9(4).
        77 n   pic 9(4).
        77 y   pic 9(16).
       procedure division using a n returning y.
      * computes r = i to the power of n
         move n to k
         move a to p
         move 1 to y
         perform until k <= 0
           if function mod(k,2)=0 then
             multiply p by p
             divide 2 into k
           else
             multiply p by y
             subtract 1 from k
           end-if
         end-perform.
        end program pow-fast.