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