A method for producing a computationally efficient system that reduces the number of iterations required to generate a conductivity image pattern of a subsurface object, and its attendant conductivity distribution, through a solution to the system of field equations that simultaneously satisfies all of the boundary conditions and conserves internal current flux densities.