Upon comparison, all three methods are shown to produce consistent predictions of equilibrium copolymerization when given the same input parameters across a vast input space of copolymerization reactivity, comonomer feed composition, and reaction rates. Correspondingly, all three methods can accurately fit model equilibrium copolymerization data. While all three methods are consistent, we make a recommendation to use numerical integration of a copolymer equation as it is computationally most efficient and its results are not affected by additional parameters such as number of equations solved, targeted chain length, or number of chains simulated. Finally, the utility of the methods developed herein is demonstrated through the fitting of seven copolymerization reactions, including two which have previously been fit under invalid assumption of irreversibility for one or both comonomers. For these cases, use of an accurate method provides statistically significantly different measurements of rate coefficients that agree with canonical knowledge of comonomer reactivities. Altogether, this work advances the understanding of equilibrium copolymerization and establishes the accuracy and precision of easyto-use methods to measure accurate reactivity information from equilibrium copolymerization data.