);
0336       <span class="keyword">else</span>
0337         <span class="comment">% expand by divs=cut(4*mu2,[2,4])</span>
0338         divs = max(2.,min(4.,4.*mu2));
0339         alp = alp2*divs;
0340       <span class="keyword">end</span>
0341       
0342       bSafeguard = 0;
0343     <span class="keyword">else</span>
0344       <span class="comment">% interpolation phase</span>
0345       <span class="keyword">if</span> (alp2 == 0.)
0346         <span class="comment">% only long steps, mu&lt;1/2, contr