function [c,ct] = BR_res_convergence(xt)
%

len_xt = length(xt(:,1));
ct = zeros(len_xt,1);
x_star = xt(end,:);

for e = 1:len_xt
    x = xt(e,:);
    ct(e) = abs(norm(x-x_star));
end
c = ct(2:end)./ct(1:end-1);
