in code/figure_1/penalties.py [0:0]
def ls(x, y, reg=.1): return Ridge(alpha=reg, fit_intercept=False).fit(x, y).coef_