Numerical integration is one of the key techniques for numerical methods based on the nodes calculation. In order to meet the calculation requirements of the true meshless methods
an adaptive numerical integration algorithm for different integration domains associated with the position of the nodes is proposed. The method can be used to determine the integration region and its Gaussian integral points automatically. Two examples
i.e. a slope body which is up against vacancy and a square plate with a hole
are analyzed. The results show that the program based on the algorithms proposed have good boundary re-shape character and high accuracy. The program can be used as a general program module for the numerical methods and applied to practice. It is easily to interpret the program to other computer languages since its basic algorithms.