matrix_multiply.py

import sys, random
from time import time

n = 4096

A = [[random.random() 
        for row in range(n)] 
            for col in range(n)]

B = [[random.random() 
        for row in range(n)] 
            for col in range(n)]

C = [[0
        for row in range(n)] 
            for col in range(n)]

start = time()
for i in range(n):
    for j in range(n): 
        for k in range(n):
            C[i][j] += A[i][k] * B[k][j]
end = time()

print(round(end - start, 6), "seconds")