This paper presents a new method of numerical integration combining properties of classical numerical methods and stochastics ones. It discusses an estimator of the numerical value of a one-dimensional integral which is exact for polynomials of a certain degree and which is unbiased under repeated use.