       % n = dimension of x (#cols of A)
if ~exist('m','var'), m=10*n; end               % m = dimension of y (#rows of A)
if ~exist('k0','var'), k0=n; end                % k0 = true sparsity (generative)
if ~exist('k','var'), k=max([2 fix(n/3)]); end  % k = desired sparsity (recovered)
if ~exist('noise','var'), noise = 0.1; end      % noise = std dev of additive noise
if ~exist('nCases','var'), nCases = 20; end    % nCases = no. of MC trials
if ~exist('newA','var'), newA = 0; end          % newA = mi