function pts = BR_MP_gen_sols(step,thresh)
%

pts = [];

for x = 0:step:1
    for y = 0:step:1
        Se = norm(MPg([0.7;x;y]));
        [p1,p2,p3,p4] = BR_cond2ft(y,0.7,x);
        p_min = min([p1,p2,p3,p4]);
        p_max = max([p1,p2,p3,p4]);
        if ~isempty(pts)&length(pts(:,1))>=119
            tch = 0;
        end
        if Se<thresh&p_min>=0&p_max<=1
            pts = [pts; x,y,Se];
        end
    end
end
