gth(Mod.z)-1;
0054 <span class="comment">% Model estimation</span>
0055 minKonf=min(abs(N.k)); 
0056 <span class="comment">% Sigma_q=Mittelwert der geringsten Eindringtiefe o. Mittel</span>
0057 rq=median(N.r(find(abs(N.k)==minKonf)));
0058 <span class="comment">%rq=mean(N.r);</span>
0059 <span class="keyword">if</span> rq&gt;20, rq=round(rq/10)*10; <span class="keyword">end</span>
0060 Mod.Bg=ones(size(Mod.z))*rq;
0061 mode=1;
0062 <span class="keyword">switch</span> mode,
0063 <span class="keyword">case</