wA = mix new &amp; old A matrices</span>

figure
xr = [-6.5:0.1:6.5];  <span class="comment">% discretized x range</span>
x = zeros(n,1);
Px = zeros(1,length(xr));

<span class="comment">% base design matrix</span>
A0 = randn(m,n); <span class="comment">% Gaussian A</span>
<span class="comment">%A0 = rand(m,n);  % Uniform A</span>
<span class="comment">%A0 = sign(randn(m,n));  % Binary A</span>

<span class="keyword">if</span> noise==0, noise = 1e-6; <span class="keyword">end</span>  <span class="comment">%