Many widely-used polymers are made via free-radical polymerization. Mathematical models of polymerization reactors have many applications such as reactor design, operation, and intensification. The method of moments has been utilized extensively for many decades to derive rate equations needed to predict polymer bulk properties. In this article, for a comprehensive list consisting of more than 40 different reactions that are most likely to occur in high-temperature free-radical homopolymerization, moment rate equations are derived methodically. Three types of radicals—secondary radicals, tertiary radicals formed through backbiting reactions, and tertiary radicals produced by intermolecular chain transfer to polymer reactions—are accounted for. The former tertiary radicals generate short-chain branches, while the latter ones produce long-chain branches. In addition, two types of dead polymer chains, saturated and unsaturated, are considered. Using a step-by-step approach based on the method of moments, this article guides the reader to determine the contributions of each reaction to the production or consumption of each species as well as to the zeroth, first and second moments of chain-length distributions of live and dead polymer chains, in order to derive the overall rate equation for each species, and to derive the rate equations for the leading moments of different chain-length distributions. The closure problems that arise are addressed by assuming chain-length distribution models. As a case study, β-scission and backbiting rate coefficients of methyl acrylate are estimated using the model, and the model is then applied to batch spontaneous thermal polymerization to predict polymer average molecular weights and monomer conversion. These predictions are compared with experimental measurements.