bsp;&nbsp;&nbsp;&nbsp;
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=comment>%&nbsp;finds&nbsp;out&nbsp;label&nbsp;for&nbsp;the&nbsp;interval&nbsp;[th1,th2].
</span><br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;class=classify((th2+th1)/2,p1,p2,m1,m2,c1,c2);
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=keyword>if</span>&nbsp;class==2
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=comment>%&nbsp;integral&nbsp;eps2&nbsp;=&nbsp;int_th2^