import numpy as np matA = [[1, 2], [3, 4]] matB = [[5, 6], [7, 8]] arrKP = np.kron(matA, matB) print (arrKP)