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