Choice of specific interpolation points in quadrature formula for integrals containing Bessel functions

A method for calculation of integrals containing Bessel functions is suggested which based on linear interpolation applied to the function entering the integral along with the Bessel function. A specific choice of interpolation points allows excluding the contribution into integral of all interpolation points besides the first one. The study of errors leads to determination of bounds of the considered integral.