Introduction - If you have any usage issues, please Google them yourself
gradient_QU.m
For example, invoking the partial sparse-simplex method with the same setting as the last
example (of the greedy sparse-simplex method) also results with the correct solution:
>>[X,fun_val]=partial_sparse_simplex(@(x)f_QU(A,c,x),@(x)gradient_QU(A,c,x),...
@(x,S)g_QU(A,c,x,S),p,500,[1 1 1 ones(7,1)])