function res = BR_2Dprobe3(step)
%

THRESH = 0.001;

res = [];

P = [0.0000
    0.0000
    0.2999
    0.0000
    0.2500
    0.2500
    0.2001
    0.0000];

vals = [0:step:1];
num_vals = length(vals);

% [0.7] A, [0.5] ~AvB, [0.5] ~BvC
V = [0 0 0 0 1 1 1 1; 1 1 1 1 0 0 1 1; 1 1 0 1 1 1 0 1; 1 1 1 1 1 1 1 1];
pc = [0.7;0.5;0.5;1];
index = 0;

wb = waitbar(0,'Probe');
for ind1 = 1:num_vals
    ind1
    waitbar(ind1/num_vals);
    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);
                for ind5 = 1:num_vals
                    p(5) = vals(ind5);
                    for ind6 = 1:num_vals
                        p(6) = vals(ind6);
                        for ind7 = 1:num_vals
                            p(7) = vals(ind7);
                            for ind8 = 1:num_vals
                                p(8) = vals(ind8);
                                pp = V*p';
                                if abs(sum(p)-1)<THRESH&norm(pp-pc)<THRESH
                                    res = [res; p,sum(p(5:8)),...
                                        sum(p(2:3))+sum(p(7:8)),...
                                        sum(p(2:2:8))];
                                end
                            end
                        end
                    end
                end
            end
        end
    end
end
close(wb);
