This paper presents a mathematical model to describe the evolution of the molecular weight distribution (MWD) in vinyl chloride (VCM) free‐radical suspension polymerizations performed with a bifunctional initiator, 1,3‐di(2‐neodecanoylperoxyisopropyl) (DIPND). The model yields, as a function of time, the mass balances for the distinct phases, the monomer conversion, the number‐ and mass‐average molecular weights and the complete MWD of both the growing and dead polymer chains. In order to describe the MWD, the model uses probability generating functions (pgf) to transform the mass balance equations into a reduced and finite set of model equations. As shown throughout many examples, the MWD's of the final polymer resin is little sensitive to the presence of the linear symmetrical bifunctional initiator.