inf&nbsp;+&nbsp;int_{-inf}^th1
</span><br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;eps2&nbsp;=&nbsp;1&nbsp;+&nbsp;erfc2(th2,m2,sqrt(c2))&nbsp;-&nbsp;erfc2(th1,m2,sqrt(c2));&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;
<br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;<span class=comment>%&nbsp;integral&nbsp;eps1=&nbsp;int_th1^th2
</span><br>
&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;eps1&nbsp;=&nbsp;erfc2(th1,m1,sqrt(c1))&nbsp;-&nbsp;erfc2(th2,m1,sqrt(c1));&