A mathematical model is developed for the molecular weight distribution (MWD) of free-radical styrene polymerization in a simulated semi-batch reactor system. The generation function technique and moment method are employed to establish the MWD model in the form of Schultz-Zimmdistribution. Both static and dynamic models are described in detail. In order to achieve the closed-loop MWD shaping by output probability density function (PDF) control, the dynamic MWD model is further developed by a linear B-spline approximation. Based on the general form of the B-spline MWD model, iterative learning PDF control strategies have been investigated in order to improve the MWD control performance. Discussions on the simulation studies show the advantages and limitations of the methodology