function z=qab(A,b,x)
z=0.5*x'*A*x-x'*b;
end