Python numba matrix multiplication

numbamatmul.py

import time
import numba
from numba import jit
import numpy as np

#input matrices
matrix1 = np.random.rand(30,30)
matrix2 = np.random.rand(30,30)
rmatrix = np.zeros(shape=(30,30))

#multiplication function
@jit('void(float64[:,:],float64[:,:],float64[:,:])')
def matmul(matrix1,matrix2,rmatrix):
for i in range(len(matrix1)):
for j in range(len(matrix2[0])):
for k in range(len(matrix2)):
rmatrix[i][j] += matrix1[i][k] * matrix2[k][j]

#Calculate running time                                                                                                             start=time.clock()
matmul(matrix1,matrix2,rmatrix)
end=time.clock()

#print results
print end-start
for r in rmatrix:
print(r)