I\'m having trouble with the matrix exponential calculation using scipy.linalg.expm.
The code is like this:
linA.expm(sp.kron(number(N), np.identity(2)