For the effective analysis of an engineering problem, meshless methods which require only positioning finite points without the element meshing recently have been proposed and being studied extensively. Meshless methods have difficulty in imposing essential boundary conditions directly, because non-interpolate shape functions originated from an approximation process are used. So some techniques, which are Lagrange multiplier method, modified variational principles and coupling with finite elements and so on, were introduced in order to impose essential boundary conditions. In spite of these methods, imposition of essential boundary conditions have still many problems like as non-positive definiteness, inaccuracy and negation of meshless characteristics. In this paper, we propose a new method which modifies shape function. Through numerical tests, convergence, accuracy and validity of this method are compared with the standard EFGM which uses Lagrange multiplier method or modified variational principles. According to this study, the proposed method shows the comparable accuracy and efficiency.