The inverse magnetoencephalography and electroencephalography problems for spherical models have been extensively discussed in the literature. Using the spherical multiple-shell model, we derive novel vector-valued and singularity-free integral equations for both problems based on the quasi-static Maxwell's equations. These equations are solved via a Fourier series expansion. We call this procedure the Edmonds approach, since an orthonormal system based on the Edmonds-vector-spherical harmonics is used for the Fourier series. Employing the associated singular-value decomposition, we provide a complete answer to the non-uniqueness question of these two problems: only the harmonic part of the solenoidal component of the neuronal current is visible to the simultaneous use of magnetoencephalography and electroencephalography. The remaining components of the current are invisible to both techniques. We state Picard's condition for the existence of a solution and derive an explicit formula for the best-approximate solution of the neuronal current. In comparison to previous approaches, the Edmonds approach requires the fewest a-priori assumptions on the neuronal current. Finally, we show that the results obtained by means of the Edmonds approach are consistent with results derived earlier via the Helmholtz decomposition.