_ = metbeta.plot(
            arange(1, n_iter + 1),
            medianbd,