Although, the finite element method numerically is efficient it exhibits difficulties whenever the remeshing of the analysis domain must be performed. For such problems, utilizing meshless computation methods is very promising. But, their large computational cost, which arises from following a time-consuming process to find shape functions, is one of the most important factors limiting the use of these methods. In this paper, we introduce a direct approach, based on properties required for any shape function, to prepare the shape function and propose a new meshless method. The proposed method does not need any predetermined basis or weighting functions and can be performed faster and more efficiently. Another advantage of the introduced method is its capability to apply desirable changes to the shape function. Application of the proposed approach for electrostatic field computation and verification of the obtained results using theoretically known solution is also presented. Copyright © 2010 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]