Although the head is more closely represented as an ellipsoid than a sphere, calculation in ellipsoidal coordinates is difficult. This paper presents a four shell ellipsoidal model, employing multipole expansion in ellipsoidal coordinates, for EEG, MEG, and evoked potential applications. Computational detail and insight into efficient calculation of the Lamé functions of the first and second kind are provided to demonstrate feasibilty. The Lamé function of the second kind, derived from the Lamé function of the first kind, can be computed at higher degrees by means of partial fraction expansion.