def naive_matrix_product(A, B):
    """ Naive square matrix multiplication """
    msize = len(A)
    nsize = len(B[0])
    # determine zero matrix
    C = [[0 for _ in range(nsize)] for _ in range(msize)]
    for i in range(msize):
        for j in range(nsize):
            for k in range(msize):
                C[i][j] += A[i][k] * B[k][j]
    return C