A Boundary Element Method (BEM) solver based on the solution of boundary integral equations of potential and electric field has been developed to simulate 3D electrostatic configuration in gaseous detectors. Use of analytical solution of the integral equations governing electric potential and field for estimating influence coefficients of the BEM solver has empowered it to provide extremely precise estimates of the potential and field for a given geometry. The nearly exact BEM (neBEM) solver has been implemented in order to simulate physical and weighting potential and field configurations in several gaseous detectors like MultiWire Proportional Counter and Time Projection Chamber. The efficacy of the solver for simulating 3D electrostatic configuration in composite systems containing both conductors and layered dielectrics has been demonstrated for some of the MicroPattern Gas Detectors and Resistive Plate Chamber. It should be noted that the method treats the dielectric interfaces to be in a steady state with polarization charges only. The reasons why the neBEM can be a preferred tool for electrostatic simulation of gaseous detectors to other 2D or 3D numerical solvers are discussed on the basis of present results.