The numerical accuracy of the boundary element (BE) method used to solve the volume conduction problem of nested compartments, each having a homogeneous conductivity, is studied. The following techniques for improving this accuracy are discussed: the handling of the auto solid angle element omega ii, the overall refinement of the level of discreteness, the use of a locally refined discrete grid, the isolated problem approach, and an adaptive refined computation of the discrete surface integrals involved in the BE method. The effects of these techniques on the numerical accuracy of the computed electrical potentials are illustrated by taking a volume conductor consisting of four concentric spheres representing the head since for this model an analytical (exact) solution is available. The techniques are of importance for numerically computed electroencephalograms (EEG's) since the numerically computed surface EEG's are severely affected by the relatively low conductivity of the compartment representing the skull.