[en] Boynton's mathematical simulation procedure for multi-component distillation calculations has the advantage that the Jacobian matrix is calculated analytically. The purpose of the present study is to adapt this procedure to hydrogen isotope separation columns by cryogenic distillation. The Boynton's model is modified so that the model can incorporate decay heat of tritium, nonideality of the hydrogen isotope solutions, multiple feeds and multiple sidestreams. Basic equations are derived and the mathematical simulation procedure is briefly explained. (author)