from numpy import zeros,dot
from numpy.matlib import rand

N = 5

A = rand(N,N)

B = rand(N,N)

C = zeros([N,N],float)

for i in range(N):

    for j in range(N):

        for k in range(N):
            
            C[i,j] += A[i,k]*B[k,j]
        
print(A)
print(B)
print(C) 
print(dot(A,B))