identification division.
       program-id. newton.
      * Autor JD
      * Dies ist das Newtonverfahren
      * gesucht x mit f(x)=0
      * es konvergiert nicht immer!!
       author. "JD".
       date-written. 27.12.2006.
       date-compiled.
       data division.
       working-storage section.
         77 x         pic 9(4)V9(8).
         77 f         pic 9(4)V9(8).
         77 d         pic 9(4)V9(8).
         78 eps       value 1e-8.
         77 xn        pic 9(4)V9(8).
         77 counter   pic 9(8).
       procedure division.
         display "start"
         move eps to f
         move 0 to counter
         accept xn
         perform with test after
         until f*f<eps or counter >1E5
           move xn to x
           perform func
           perform funcdiff
           compute xn = xn -f / d
           add 1 to counter
         end-perform.
         display "done, Zero(x**4+3x**2+5)="xn.
         stop run.
        func.
          compute f=x*x*x*x*x+3*x*x+5
        funcdiff.
          compute d=4*x*x*x+6*x
        end-program newton.