0010 X=zeros(length(A(1,:)),length(alpha));
0011 R=X;
0012 
0013 <span class="comment">%% outer variables</span>
0014 uu=y;         <span class="comment">%% u=\tilde u</span>
0015 vv=(uu'*A)';     <span class="comment">%% v=\tilde v</span>
0016 beta=norm(vv);
0017 v=vv/beta; u=uu/beta;
0018 p=vv; r=vv; uold=0;
0019 maxss=norm(r)*1e-6;
0020 
0021 kmax=length(alpha)-1;
0022 <span class="comment">% P anders</span>
0023 <span class="comment">%shifted=0:kmax;</span>
0024 <span class="keyword">for</span> k = 0:km