Abstract: The key facility for producing ethyl acetate is a complex distillation column with a reboiler performing as an esterification reactor. An algorithm, employing Block-Hegners algorithm together with a relaxation method has been taken for the simulation of this reactor-three phase distillation column system. More stable convergence has thus been obtained. Furthermore, by using the sparse matrix approach, less storage requirements are needed, allowing the problem to be solved in a microcomputer. Twelve case studies are conducted using this algorithm. Calculated results are compared with data measured in an operating plant.