global numf numg numH

load hwk9 A G
b = ones(size(A,1),1);
Q = G; d = ones(size(A,2),1);

options = optimset('quadprog');
options.Algorithm = 'active-set';

% tic, [x,....] = quadprog(...); toc,

% now form KKT matrix as K and rhs as r
% y = K\r;
% compare differences

% now do more general qp
load hwk8 Q d A b ineq eq

% solve using various algorithm options for quadprog
% give option settings for the best one