,0);
rand(<span class="string">'state'</span>,0);

A0 = randn(m,n);
[U,S,V] = svd(A0);
S= diag(fliplr(logspace(-0.7,1,n)));
A0 = U(:,1:n)*S*V';
A1 = randn(m,n);  A1 = A1/norm(A1);
A2 = randn(m,n);  A2 = A2/norm(A2);

Aperb0 = [A1;A2];
p = 2;

b = U(:,1:n)*randn(n,1) + .1*randn(m,1);

<span class="comment">% we consider LS problems || (A0 + u1*A1 + u2*A2) x - b||</span>
<span class="comment">% where  ||u|| leq rho</span>

<span class="comment">% Nominal Solution</span>
xnom = A0\b;

<span class="comment">% T