We introduce the function $Z(x; \xi, \nu) := \int_{-\infty}^x \varphi(t-\xi)\cdot \Phi(\nu t)\ \text{d}t$, where $\varphi$ and $\Phi$ are the pdf and cdf of $N(0,1)$, respectively. We derive two recurrence formulas for the effective computation of its values. We show that with an algorithm for this function, we can efficiently compute the second-order terms of Bonferroni-type inequalities yielding the upper and lower bounds for the distribution of a max-type binary segmentation statistic in the case of small samples (where asymptotic results do not work), and in general for max-type random variables of a certain type. We show three applications of the method --- (a)~calculation of critical values of the segmentation statistic, (b)~evaluation of its efficiency and (c)~evaluation of an estimator of a point of change in the mean of time series.
Bonferroni inequality, segmentation statistic, Z-function
62E17, 05A20