% quadprogtest
H = [   1003.1  4.3     6.3     5.9;
        4.3     2.2     2.1     3.9;
        6.3     2.1     3.5     4.8;
        5.9     3.9     4.8     10  ];
c = zeros(4,1);

A = [ -0.17 -0.11 -0.10 -0.18 ];
b = -0.10;

Aeq = [ 1 1 1 1 ];
beq = 1;

lb = zeros(4,1);
ub = inf(4,1);
x0 = [1; 0; 0; 1];

[x, f, exitflag, output, lambda] = quadprog(H, c, A, b, Aeq, beq, lb, ub, x0);

% x
% f
% lambda.lower
% lambda.upper
lambda.eqlin
% lambda.ineqlin
