# credit to the Stack Overflow user in the source link for this example

import pandas as pd
import numpy as np

shape = (50, 4460)
data = np.random.normal(size=shape)
data[:, 1000] += data[:, 2000]

df = pd.DataFrame(data)
c = df.corr() # the correlation matrix
s = c.unstack()
so = s.sort_values(kind = "quicksort") # series containing all possible  couples