The element free Galerkin (EFG) meshless method is very effective for solving electrostatic problems. However, it leads to shape functions which do not satisfy the Kronecher delta property. Therefore, it is necessary to use additional techniques to impose the essential boundary conditions (EBCs) and to treat the material discontinuities. In this work, the EFG method is used with interpolating shape functions. Therefore, the Kronecher delta property is satisfied, the EBCs are enforced directly in the discrete system, and the material discontinuities treatment is improved. The method is evaluated in the analysis of a cylindrical coaxial capacitor filled with two concentric and different dielectric materials where obtained results confirmed its efficiency. © 2015 Wiley Periodicals, Inc. Microwave Opt Technol Lett 57:1390–1395, 2015