,k+1)=ss(k+1)*v;
0052     P(:,k+1)=R(:,k+1)+ww(k+1)*P(:,k+1);
0053   <span class="keyword">end</span>
0054   <span class="comment">%shifted(kmin+1:kmax+1)=iter;</span>
0055   ok=kmin;
0056   kmin=max(find(abs(ss)&lt;maxss));
0057   <span class="keyword">if</span> kmin&gt;ok, waitbar(kmin/kmax,wb); <span class="keyword">end</span>
0058   <span class="keyword">if</span> isempty(kmin), kmin=0; <span class="keyword">end</span>
0059 <span class="keyword">end</span>
0060 close(wb);
0061 <span class="comment">%sav