We investigate a sequence of polynomials orthogonal with respect to the weight function, which appears in the error statistics in high speed optical fiber communication systems. Computation of these integrals by classical approach shows unexpectedly dubious accuracy. That is why we had to construct the new quadratures, to get more reliable results in our applications. Numerical procedure for quadrature construction is considered and its stability is discussed in detail.