# 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