function [c,ceq] = exactSol(x,Solution,i,j)
for k = 1:length(j)
    c(k) = x'*(Solution(i).Z.Q-Solution(j(k)).Z.Q)*x + ...
        (Solution(i).Z.c-Solution(j(k)).Z.c)'*x + ...
        Solution(i).Z.d-Solution(j(k)).Z.d - 1e-3;
end
ceq = [];
end