 sets = subsets1(1:n,kk);
   Ns = length(sets);
   Xmode = zeros(Ns,n);
   Emode = zeros(Ns,1);

   <span class="keyword">for</span> c = 1:nCases   <span class="comment">% ------------------------------------- MC loop</span>

       <span class="comment">% sample a (blend of new &amp; old) design matrix</span>
       A = (1-newA)*A0 + newA*(randn(m,n));  <span class="comment">% Gaussian update</span>
       <span class="comment">%A = (1-newA)*A0 + newA*(rand(m,n));  % Uniform update</span>

       <span cla