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,