Identification Division.
       Program-Id. Ising.
      * from Stauffer/Stanley
      * From Newton to Mandelbrot
      * Springer 1996
       Data Division.
       Working-Storage Section.
         01 is        pic S9(4)V9(4) occurs 1680.
         77 vx          pic S9(4)V9(4).
         77 vy          pic S9(4)V9(4).
         77 dt          pic S9(4)V9(4).
       Procedure Division.
        move 40 to L
        move 0.2 to p
        compute L1=L+1
        compute Lp=L*L+L
        compute Lm=Lp+L
        perform varying i from 1 until i>Lm
          move -1 to is(i)
          if function rnd(i)<p then move 1 to is(i)
        end-perform
        perform varying it from 1 until it>100
          move 0 to m
          perform varying i from L1 until i>Lp
            if is(i-1)+is(i+1)+is(i-L)+is(i+L)=0 then
              move -is(i) to is(i)
            add is(i) to m
          end-perform
          display it, m
        end-perform
       End-Program Ising.