This function solves linear least squares problem with nonnegative variables using the active set algorithm described in: http://www.mathworks.com/matlabcentral/fileexchange/10908-active-set-algorithm