N = 2;
H = rand (N);
Q = zeros (N, 1);
A = ones (1, N);
B = 1;
LB = zeros (N, 1);
# This works
UB = ones (N, 1) / N;
qp ([], H, Q, A, B, LB, UB);
# This triggers an error
UB = ones (N, 1) / N - 0.001;
qp ([], H, Q, A, B, LB, UB);

