function [KB00,KB01,KB10,KB11] = BR_incon
%

THRESH = 0.001;

KB(1).clauses = [1,3];
KB(2).clauses = [-1,2];
KB(3).clauses = [-2,3];
KB(1).prob = 0;
KB(2).prob = 0;
KB(3).prob = 0;

KB00 = [];
KB01 = [];
KB10 = [];
KB11 = [];

for p1 = 0.1:0.1:0.9
    KB(1).prob = p1;
    for p2 = 0.1:0.1:0.9
        KB(2).prob = p2;
        for p3 = 0.1:0.1:0.9
            if p1>0.6&p2>0.6&p3>0.3
                tch = 0;
            end
            KB(3).prob = p3;
            P = BR_Nilsson(KB);
            P_sum = sum(P);
            PS1 = BR_Nilsson_query_clause(KB(1).clauses,P);
            PS2 = BR_Nilsson_query_clause(KB(2).clauses,P);
            PS3 = BR_Nilsson_query_clause(KB(3).clauses,P);
            cond1 = abs(P_sum-1)<THRESH;
            cond2 = abs(PS1-KB(1).prob)<THRESH;
            cond3 = abs(PS2-KB(2).prob)<THRESH;
            cond4 = abs(PS3-KB(3).prob)<THRESH;
            X = BR_linprog(KB);
            if ~(cond1&cond2&cond3&cond4)&isempty(X)
                KB00 = [KB00;p1,p2,p3];
            end
            if ~(cond1&cond2&cond3&cond4)&~isempty(X)
                KB01 = [KB01;p1,p2,p3];
            end
            if (cond1&cond2&cond3&cond4)&isempty(X)
                KB10 = [KB10;p1,p2,p3];
            end
            if (cond1&cond2&cond3&cond4)&~isempty(X)
                KB11 = [KB11;p1,p2,p3];
            end
        end
    end
end
