A modified Pagano method is developed for the three-dimensional (3D) free vibration analysis of simply-supported, multilayered functionally graded material (FGM) circular hollow cylinders with a constant rotational speed with respect to the meridional direction of the cylinders. The material properties of each FGM layer constituting the cylinders are regarded as heterogeneous through the thickness coordinate, and then specified to obey a power-law distribution of the volume fractions of the constituents, and the effects of centrifugal and Coriolis accelerations, as well as the initial hoop stress due to rotation, are considered. The Pagano method, which was developed for the static and dynamic analyses of multilayered composite plates, is modified in that a displacement-based formulation is replaced by a mixed formulation, the complex-valued solutions of the system equations are transferred to the real-valued solutions, a successive approximation method is adopted to extend its application to FGM cylinders, and a propagator matrix method is developed to reduce the time needed for its implementation. These modifications make the Pagano method feasible for multilayered FGM cylinders, and the computation in the implementation is independent of the total number of the layers, thus becoming less time-consuming than usual.