Abstract: In this paper, we present an analytical T -matrix solution of Mie type for electromagnetic (EM) scattering by a radially multilayered gyromagnetic anisotropic sphere. A sphere located in free space is divided into N regions, each region is characterized by scalar permittivity and gyromagnetic permeability tensor. Based on the vector spherical wave functions (VSWFs) theory, the incident field, scattered field, and the field inside the sphere are expanded as linear combination of VSWFs. The scattering coefficients and unknown expansion coefficients of each layer are calculated by matching the boundary conditions on the discontinuous surfaces. According to the derived formulations, the numerical results of radar cross sections (RCS) given by proposed analytical method is in good agreement with those computed by finite element method (FEM). Several calculation examples for some specific cases are also presented to investigate the EM scattering characteristic of multilayered gyromagnetic sphere.