ax, P(:,k+1)=vv; <span class="keyword">end</span>
0025 R=P;
0026 
0027 <span class="comment">%% inner variables</span>
0028 gg=ones(kmax+1,1);     <span class="comment">% \hat\gamma</span>
0029 ww=zeros(kmax+1,1);     <span class="comment">% \hat\omega</span>
0030 ss=beta*ones(kmax+1,1);  <span class="comment">% \hat\sigma Knackpunkt</span>
0031 kmin=0;
0032 iter=0;
0033 wb=waitbar(0,<span class="string">'Shifted CG Inversion...'</span>);
0034 <span class="keyword">while</span> kmin&lt;kmax,   
0035   iter=