def get_random_data(n=50, seed=0): np.random.seed(seed) data = np.random.choice(10000, (n, 1), replace=False) D = pdist(data, 'euclidean') Z = linkage(D, 'ward') return data, Z, D