A local–node numerical integration scheme for meshless methods is presented in this work. The distinguishing characteristic of the introduced scheme is that not support mesh or grid to perform numerical integration is needed, besides the fact that gauss cubature points for each node are generated in a properly fashion and that the extension of the methodology to high dimensions is straightforward. The numerical integration is computed with the Gauss–Hermite cubature formulas, and the partition of unity is employed to introduce the Gaussian weight in a natural way. Selected numerical tests in two-dimensions are used to illustrate the validity of the proposed methodology. Although the obtained results are encouraging, the behavior of the integration error is not still well understood when the dimensionless parameter which control the width of the Gaussian kernel varies.