>>> import numpy as np

>>> matrixA = np.matrix([
... [0, 0, 0, 0, 1, 0, 1, 0, 0, 0],
... [0, 0, 0, 0, 0, 0, 1, 0, 1, 0],
... [0, 0, 0, 0, 0, 0, 0, 1, 0, 1],
... [0, 0, 0, 0, 1, 0, 0, 0, 1, 0],
... [1, 0, 0, 1, 0, 0, 0, 0, 0, 0],
... [0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
... [1, 1, 0, 0, 0, 0, 0, 1, 0, 0],
... [0, 0, 1, 0, 0, 0, 1, 0, 0, 0],
... [0, 1, 0, 1, 0, 0, 0, 0, 0, 0],
... [0, 0, 1, 0, 1, 0, 0, 0, 0, 0]])

>>> matrixA
matrix([[0,  0,  0,  0,  1,  0,  1,  0,  0,  0],
        [0,  0,  0,  0,  0,  0,  1,  0,  1,  0],
        [0,  0,  0,  0,  0,  0,  0,  1,  0,  1],
        [0,  0,  0,  0,  1,  0,  0,  0,  1,  0],
        [1,  0,  0,  1,  0,  0,  0,  0,  0,  0],
        [0,  0,  0,  0,  0,  0,  0,  0,  0,  0],
        [1,  1,  0,  0,  0,  0,  0,  1,  0,  0],
        [0,  0,  1,  0,  0,  0,  1,  0,  0,  0],
        [0,  1,  0,  1,  0,  0,  0,  0,  0,  0],
        [0,  0,  1,  0,  1,  0,  0,  0,  0,  0]])

>>> matrixB = matrixA

>>> matrixC = matrixA*matrixB # The Matrix Multiplication bit

>>> matrixC
matrix([[2,  1,  0,  1,  0,  0,  0,  1,  0,  0],
        [1,  2,  0,  1,  0,  0,  0,  1,  0,  0],
        [0,  0,  2,  0,  1,  0,  1,  0,  0,  0],
        [1,  1,  0,  2,  0,  0,  0,  0,  0,  0],
        [0,  0,  0,  0,  2,  0,  1,  0,  1,  0],
        [0,  0,  0,  0,  0,  0,  0,  0,  0,  0],
        [0,  0,  1,  0,  1,  0,  3,  0,  1,  0],
        [1,  1,  0,  0,  0,  0,  0,  2,  0,  1],
        [0,  0,  0,  0,  1,  0,  1,  0,  2,  0],
        [1,  0,  0,  1,  0,  0,  0,  1,  0,  1]])