c Evgenii B. Rudnyi, http://Evgenii.Rudnyi.Ru
#ifndef dim1
#define dim1 1000
#endif
#ifndef dim2
#define dim2 dim1
#endif
#ifndef dim3
#define dim3 dim1
#endif
      double precision A(dim1, dim2), B(dim2, dim3), C(dim1, dim3)
      real start, finish, tarray(2)
      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      
      call etime(tarray, start)
      do j = 1, dim3
        do i = 1, dim1
          C(i, j) = 0
        enddo
        do k = 1, dim2
          do i = 1, dim1
            C(i, j) = C(i, j) + A(i, k)*B(k, j)
          enddo
        enddo
      enddo
      call etime(tarray, finish)
      print *,'time for C(',dim1,',',dim3,') = A(',dim1,',',dim2,') B(',
     1dim2,',',dim3,') is',finish - start,' s'
      end