step or short step, overwrite alp2</span>
0237       alp2 = alp;   
0238       f2 = f;
0239       mu2 = mu;
0240       omega2 = 2*abs(1.-mu2)/alp2;
0241     <span class="keyword">end</span>
0242     
0243     <span class="keyword">if</span> (f2 == f0 &amp;&amp; f1 &gt;= f0 &amp;&amp; k &gt;= 2)
0244       <span class="comment">% check whether f=f0 steps should be considered long</span>
0245       <span class="comment">% ksh=7 indicates f=f0 for all large alp</span>
0246       <span class="comment">% f1=f0 &