function [X,lambda,how] = qpemulation(H,f,A,B,vlb,vub,X,neqcstr,verbosity,varargin)

% vlb is assumed to be zero
if neqcstr > 0 
   error('Needs QP algorithm')
end

n = size(H,1);
k = cond(H);
newH = [0.5*H+k*eye(size(H)) k*eye(size(H));k*eye(size(H)) k*eye(size(H))];
newc = [-0.5*f + k*vub; k*vub];
alpha = fnnls(newH, newc);
X = alpha(1:n);
