from matplotlib import pyplot as plt
from ortoolpy import tsp
M = m.coefs_[0]
N = M.shape[1]
dist = {(i, j): np.linalg.norm(M[:, i] - M[:, j])
for i in range(N)
for j in range(N)
if i != j}
_, result = tsp(range(N), dist)
plt.imshow(M[:, result])
plt.show()