Co-catalyst has several roles in propylene polymerization such as site activation, chain transfer, site transformation and Poisson deactivation (or Poisson scavenger). These functionalities strongly impact the most important of the final product properties indices namely the number and weight average. Molecular weight, polydispersity index, melt flow index and even the yield of Ziegler-Natta catalyst. Therefore, the Co-catalyst role as an independent process variable should be carefully investigated in the polymerization. In this study, investigating the co-catalyst effects on the final properties indices and the yield of the catalyst as the existing gap was defined as the aims of this study. We attempt to reply to the problems in a proper way in the absence of hydrogen and at the optimal temperature reaction by proposing a newly developed mathematical model. This model was programmed in MATLAB & Simulink software, then was validated by using experimental data. Through the model, the optimal amount of the used co-catalyst could be determined conveniently.