1&nbsp;+&nbsp;erfc2(th2,m1,sqrt(c1))&nbsp;-&nbsp;erfc2(th1,m1,sqrt(c1));&nbsp;&nbsp;
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=comment>%&nbsp;integral&nbsp;eps2=&nbsp;int_th1^th2
</span><br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;eps2=erfc2(th1,m2,sqrt(c2))-erfc2(th2,m2,sqrt(c2));
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;inter1=[th1,th2];
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span cl