A general route is shown to calculate the entropy production sigma as function of time t in a closed system during reversible polymerization. We treat the polymer molecules to behave nonideal and apply exemplarily the classical Flory-Huggins theory to get explicit expressions for the activity coefficient. At the beginning of the polymerization the system is in a nonequilibrium state where chemical reactions take place that irreversibly drive the system towards equilibrium with sigma approaching zero in the limit t-->infinity. The time-dependent course of the entropy production is explicitly calculated for two cases where the reaction starts (i) from monomer molecules polymerizing to a defined number average chain length xn,eq and (ii) from monodisperse polymer molecules reacting with each other under the constrain that xn is the same at the beginning and the end of the reaction. In both cases we find that the nature of the activity coefficient has an important effect on the curvature of sigma which may considerably differ from that of an ideal behavior.