in code/figure_1/penalties.py [0:0]
def sample(n=100000, e=1): x = np.random.randn(n, 1) * e y = x + np.random.randn(n, 1) * e z = y + np.random.randn(n, 1) return np.hstack((x, z)), y