</span>
0399             alpd = fact*alp2;
0400           <span class="keyword">end</span>
0401         <span class="keyword">end</span>
0402         
0403         alp = min(alpd,alp1/fact);
0404         <span class="keyword">if</span> (alp ~= alpd || alp &gt;= alp1 || alp&lt;=alp2)
0405           <span class="comment">% step not interior enough, use geometric mean step</span>
0406           alp = alp1*sqrt(alp2/alp1);
0407         <span class="keyword">end</span>
0408       <span class="keyword">end</span>