n> (q &gt;= sml*sml &amp;&amp; q &lt;= sml)
0382           <span class="comment">% contract by factor q &lt;= sml</span>
0383           alp = q*alp1;
0384         <span class="keyword">else</span>
0385           <span class="comment">% contract by sml</span>
0386           alp = sml*alp1;
0387         <span class="keyword">end</span>
0388       <span class="keyword">else</span>
0389         <span class="comment">% force interiority condition for interpolation step</span>
0390         alpd = alp2+q*(alp1-alp