function v = BR_f_MP(a,b);
%

v = (-0.7+a).^2 + (-0.7 + (1-a) + b - (1-a).*b).^2;
