function v = quadratic(G,c,vec);
%
% Evaluate a quadratic 
%
  v = c'*vec(:) + vec(:)'*G*vec(:) / 2;
%
% return
%
 
