A discrete-layer model is presented and applied to the free vibration of layered anisotropic spheres with coupling among the elastic, electric, and magnetic fields. Through-thickness approximations in the radial direction are pre-integrated and combined with independent approximations in the azimuthal and circumferential directions to provide estimates of the natural frequencies for a variety of representative geometries. Results are in excellent agreement with existing analytical studies and additional results are presented for higher-order spheroidal modes. Predictions of the level of influence of magnetoelectric coupling are also given. (C) 2016 Acoustical Society of America.