mul2.errorbar(
    range(1, n_nonzero_coefs + 1),
    rmse_gaussian3.mean(1),
    yerr=rmse_gaussian3.std(1),