x = sdpvar(3,1);
plot(H*x < k);hold on
plot(H*x < k,x(1:2),'k');