In this work, a robust mesh free method has been presented for the analysis of two dimensional problems. An efficient natural neighbor algorithm for construction of polygonal support domains has been used in this method that takes into account the nonuniform nodal discretization in the element free Galerkin formulation. The use of natural neighbors for determining the compact support is shown to overcome some of the shortcomings of the conventional distance metric based methods. For nonuniform nodal discretization there is a need for evaluating weights that have anisotropic compact supports. The smoothness and conformance of the weight function to the support domain obtained from natural neighbor algorithm is achieved through an efficient conformal mapping procedure such as Schwarz-Christoffel mapping. Numerical examples demonstrate that the proposed mesh free method gives good estimates of the stress/strain fields. © 2008 Springer-Verlag.