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);