Anew boundary element method (BEM) is presented for 2-D dc resistivity modeling with a point source. When compared with previously published techniques, the new method has two main features: (1)The normal derivative of potential has been eliminated from the integral equation. The formulation of the present method is simpler and requires less memory and time than the previous published methods. (2) Multiple subsurface inhomogeneous bodies can be modeled.For a simple testing model, the maximum relative error of reciprocity test is 0.24% and the average relative error is 0.05%. For the same model, the maximum relative difference between the BEM solution and the finite-element method solution is 1.14% and the average relative difference is 0.73%. For a field geoelectrical profile, the responses of the constructed model agree with the observed data quite well.