A construção de um modelo matemático da reação de copolimerização em emulsão de acrilato de butila e estireno em reator semi-batelada e isotérmico apresentou bons resultados quando comparado aos experimentos. O modelo contemplou a solução do balanço populacional, tanto para a distribuição de tamanho de partículas quando para a distribuição média de radicais por partícula. O problema descrito foi resolvido através da solução numérica de um conjunto de equações algébrico-diferenciais. A solução do modelo de partição de monômeros entre a fase aquosa e a fase partícula foi resolvida através da solução de sistema de equações não lineares a cada passo de integração, enquanto que as equações diferencias foram resolvidas pelo método de Runge-Kutta. Pôde-se validar experimentalmente a concentração de monômeros na fase aquosa, diâmetro médio de partículas, número de partículas por litro de emulsão, número médio de radicais por partícula e a distribuição de tamanho de partículas ao longo do tempo