function res = BR_2Dprobe(step)
%

THRESH = 0.001;

res = [];

vals = [0:step:1];
num_vals = length(vals);
V = [0 0 1 1; 1 1 0 1; 1 1 1 1];
pc = [0.7;0.7;1];
index = 0;

for ind1 = 1:num_vals
    ind1
    p(1) = vals(ind1);
    for ind2 = 1:num_vals
        p(2) = vals(ind2);
        for ind3 = 1:num_vals
            p(3) = vals(ind3);
            for ind4 = 1:num_vals
                p(4) = vals(ind4);
                pp = V*p';
                if abs(sum(p)-1)<THRESH&norm(pp-pc)<THRESH
                    res = [res; p,sum(p(3:4)),p(2)+p(4)];
                end
            end
        end
    end
end
