We employ the extended boundary condition method (EBCM) and construct a solution for the problem of electromagnetic (EM) scattering by anisotropic core-shell bodies of revolution (BoRs). In particular, two different core-shell configurations are examined: the gyroelectric-isotropic and the isotropic-gyroelectric setup. To construct the solution, we employ two groups of integral representations (IRs)—one group for each configuration solved—in conjunction with the discrete eigenfunction (DE) expansion of the fields in terms of spherical vector wave functions (SVWFs) for the gyroelectric regions. We demonstrate the validity and the computational performance of the method by comparisons with the HFSS commercial software for various core-shell setups such as spheroidal, cylindrical, and combined spherical-cylindrical BoRs. We also employ ADDA, a particular version of the discrete dipole approximation (DDA) method, to trace the boundaries of validity of the EBCM. Finally, we present an application of the method to the study of magnetically-tunable spheroidal THz antennas. The method can be used in a variety of potential EM applications including microwaves, functional photonics structures, as well as nanoantenna engineering.