_ = methaus.plot(
            arange(1, n_iter + 1),
            medianhfs,