home *** CD-ROM | disk | FTP | other *** search
/ OS/2 Professional / OS2PRO194.ISO / os2 / prgramer / adaptor / examples / dalib / replic / rep2.f < prev    next >
Encoding:
Text File  |  1993-04-28  |  574 b   |  34 lines
  1.       program replicate
  2.  
  3.       double precision a(10,10), ra(5,5), x
  4.       integer i, j, n
  5.  
  6. cmf$  layout ra(:serial,:serial)
  7.  
  8.       forall (i=1:10,j=1:10) a(i,j) = i + j
  9.  
  10.       ra = 0
  11.       ra = a(4:8,2:6)
  12.  
  13.       n = 0
  14.       do i = 4, 8
  15.         do j = 2, 6
  16.            n = n + i + j
  17.            print *, 'a(',i,',',j,') = ', a(i,j)
  18.         end do
  19.       end do
  20.  
  21.       do i = 1, 5
  22.         do j = 1, 5
  23.            print *, 'ra(',i,',',j,') = ', ra(i,j)
  24.         end do
  25.       end do
  26.  
  27.       x = sum (ra)
  28.  
  29.       print *, x
  30.  
  31.       print *, 'Result is ', x, ' should be ', n
  32.  
  33.       end
  34.