Computationally localizing electrical current sources of the electroen
cephalographic signal requires a volume conductor model which relates
theoretical scalp potentials to the dipolar source located within the
modeled brain, The commonly used multishell spherical model provides t
his source-potential relationship using a sum of infinite series whose
computation is difficult, This paper provides a closed-form approxima
tion to this sum based on an optimal fitting to the weights of the Leg
endre polynomials, The second-order (third-order) approximation algori
thm, implemented by a provided C-routine; requires only 100 (140) floa
ting point operations to compute a single scalp potential in response
to an arbitrary current dipole located within a four-shell spherical v
olume conductor model. This cost of computation represents only 6.3% (
8.9%) of that required by the direct method, The relative mean square
error, measured by using 20 000 random dipoles distributed within the
modeled brain, is only 0.29% (0.066%).