Methods based on Diophantine approximations lead to a simple and easy construction of rational vectors the multiple of which, reduced modulo 1, form finite sequences of points with certain properties of equipartition over the unit square. Given a function of bounded variation over this square, it is suggested that computing the average of its values at the points of such a sequence can be a practical method of numerical integration. Precise bounds for the error ore obtained. In the general case, these are of the order of the product of the reciprocal and of the logarithm of the number of points; in the case of a function satisfying stated conditions of regularity and periodicity, they are of the order of the logarithm of the number of points divided by an appropriately high power of this number. A slight sharpening of some well-known results on equipartiton is obtained incidentally.