Abstract To improve the application of Point Interpolation Method (PIM) in Element Free Galerkin Method (EFG) is the aim of this study. The trait of EFG is using overlap of influence domain between different nodes to construct discretization nodes¡¦ connection. EFG just uses nodal data, but not element. For constructing shape function, EFG has two types of methods, Fitting and Interpolation. Fitting uses Moving Least Square Method (MLS). MLS-EFG has stable effect on numerical analysis; however, users who use it need to choose more numerical parameters and do more computation. Besides, users can not apply boundary conditions directly when using MLS-EFG. Interpolation method applies nodal coordinates to proceed computation, and it called PIM. Boundary conditions could be used directly and less computation is needs while using PIM. However, the coefficient of interpolation function of sample is singular. This study tries to construct Coordination Point Interpolation Method. It owns advantages of both methods that mentioned above, and extra numerical parameters are not needed. It applies the notion of influence domain of MLS-EFG, then search correlative efficient nodes which are contained in near field of sample. The correlative efficient nodes make up matrix that con cause inverse matrix. In addition, via numerical simulations, it shows that CPIM has excellent convergence and accurate solution, and is better that MLS-EFG.