We verified the effect of the Taylor expansion order on the accuracy of the Taylor-series expansion method of moment (TEMOM) model in both the free molecular and continuum-slip regimes. The ordinary differential equations for moments with fourth-order Taylor expansion were first derived for fractal-like agglomerates, which were further compared to the existing TEMOM model with third-order Taylor expansion. We confirmed that the TEMOM model with a fourth-order Taylor expansion is less accurate than that with a third-order Taylor expansion. Moreover, the scope of application of the TEMOM model with a fourth-order Taylor expansion is limited. The existing TEMOM model with a third-order Taylor expansion was verified as the most reliable model for solving population balance equations for agglomerates undergoing Brownian coagulation.