function err = KBMPis(x)
%

F(1) = -0.7+x(1);
F(2) = -0.7+(1-x(1))+x(2)-(1-x(1))*x(2);
err = norm(F);
