R = [H*x < k]
Q = projection(R,x(1:2))
plot(R);hold on;
plot(Q,'k');