Mathematical models of reactors for the polymerization of methylmethacrylate (MMA) have been developed and analyzed to elucidate reactor dynamics and to determine conditions for improved operation. The effects of mixing and heat transfer in an MMA polymerization reactor system have been explored by the development of an imperfect mixing model. To model imperfect mixing in polymerization, a reactor configuration using two tanks in parallel was used. Bifurcation diagrams developed using numerical analysis of the model have been drawn with two variable parameters, an exchange ratio, σ, and a volume ratio, κ. We use feed and coolant temperatures as bifurcation parameters. If variable parameters are small, the lower solution branch of the steady state solutions is quite different from that of a simple model that assumes perfect macro-mixing as bifurcation parameters change. If σ increases (κ=0.1, σ=l.0), the shape of a steady state solution curve differs significantly from that of a simple model as the feed temperature decreases.