function [res_near,res_far] = BR_run_study_ind(test_KB_ind)
%

ALPHA = 0.001;
TOL = 0.001;
MAX_ITER = 40000;

num_KB = length(test_KB_ind);

for k = 1:num_KB
    k
    fname = ['KB',num2str(k),'F'];
    KB = test_KB_ind(k).KB;
    ap = test_KB_ind(k).ap;
    BR_KB2F(KB,fname);
    % near initial point
    x0 = ap' + min(1,max(0,0.1*(0.5-rand(length(ap),1))));
    [x,e,xt,et] = BR_KB_Jacobian(KB,fname,ALPHA,TOL,MAX_ITER,x0);
    res_near(k).xt = xt;
    res_near(k).et = et;
    [c,ct] = BR_res_convergence(xt);
    res_near(k).c = c;
    res_near(k).ct = ct;
    % far initial point
    x0 = ap' + min(1,max(0,1-ap));
    [x,e,xt,et] = BR_KB_Jacobian(KB,fname,ALPHA,TOL,MAX_ITER,x0);
    res_far(k).xt = xt;
    res_far(k).et = et;
    [c,ct] = BR_res_convergence(xt);
    res_far(k).c = c;
    res_far(k).ct = ct;
end
