% test with example 16.2 in Nocedal and Wright
G = [6 2 1 ;2 5 2; 1 2 4]; c = [-8;-3;-3];
A = [1 0 1; 0 1 1]; b = [3;0];
opts.tol=1e-10;
[x,la] = eqp(G,c,A,b,opts)
% answer should be x = [2;-1;1], la=[3;-2]