c Evgenii B. Rudnyi, http://Evgenii.Rudnyi.Ru
c 29.05.08 etime has been changed to cpu_time
      integer dim1, dim2, dim3
      parameter (dim1 = 1000, dim2 = 1000, dim3 = 1000)
      double precision A(dim1, dim2), B(dim2, dim3), C(dim1, dim3)
c      real start, finish, tarray(2)
      real start, finish
      call srand(86456)
      do i = 1, dim1
        do j = 1, dim2
          A(i, j) = rand()
        enddo
      enddo
      do i = 1, dim2
        do j = 1, dim3
          B(i, j) = rand()
        enddo
      enddo
c etime is not a part of the standard      
c      call etime(tarray, start)
      call cpu_time(start)
      do i = 1, dim1
        do j = 1, dim3
          C(i, j) = 0.
          do k = 1, dim2
            C(i, j) = C(i, j) + A(i, k)*B(k, j)
          enddo
        enddo
      enddo
c      call etime(tarray, finish)
      call cpu_time(finish)
      print *,'time for C(',dim1,',',dim3,') = A(',dim1,',',dim2,') B(',
     1dim2,',',dim3,') is',finish - start,' s'
      end