mul1.errorbar(
    range(1, n_nonzero_coefs + 1),
    rmse_uniform2.mean(1),
    yerr=rmse_uniform2.std(1),