Two methods of discretization, collocation method and Galerkin method, have been dominant in existing meshless methods. Although Galerkin method possesses several advantages, one of the major difficulties in the implementation of Galerkin based meshless method is how to evaluate integrals. In the Galerkin method, derivatives in domain integrals are lowered by using the divergence theorem to establish the weak form. The inaccuracy in integration will result in significant error in the solution. However, the ...