Abstract: A mathematical model and a subsequent algorithm, simulating reactive distillation processes are presented in this paper. The model takes into accountreactions in series that occur within a distillation column and the algorithm is worked out with block tridiagonal matrix technique. Suggested existences of offtridiagonal submatrices and/or submatrices on the border render the solution of the block tridiagonal matrix equation better general ised. This would be a helpful tool in design computation. It would also be useful for observing the effect of variation in flow configurations between stages on intermediate product yield. The saponii ication of propylene chloro-hydrin is illustrated to verify the employment of the presented method. Satisfactory results are obtained.