PROGRAM MatrMult 
   IMPLICIT NONE

   REAL(8), DIMENSION(1000,1000) :: A, B, C
   INTEGER :: L
   REAL ::   cput2, cput1

   INTEGER :: i, j, k

   !
   L=1000

   DO i = 1, L
      DO k = 1, L
         A(k,i) = REAL(i-k) / REAL(i+k)
         B(k,i) = 0.01 * REAL(k-i) / REAL(k+i)
      END DO
   END DO

   CALL CPU_TIME(cput1)
   
   DO j = 1, L
      DO i = 1, L
         C(i,j) = 0.0
         DO k = 1, L
            C(i,j) = C(i,j) + A(i,k) * B(k,j)
         END DO
      END DO
   END DO

   CALL CPU_TIME(cput2)
    
   WRITE(*,*) "First 10 elements of C:"
   DO i = 1, 10
      WRITE(*,*) C(i,1)
   END DO
   WRITE(*,*) "Cpu time is: ",(cput2-cput1)
!
   STOP
END PROGRAM MatrMult