After the early work by and Thiele and Geddes 1933 many methods have been described in the literature to determine all process variables of a conventional distillation column. They include the so-called θ-method , the sum rate method and the relaxation method . Also Newton type methods have been used to solve the corresponding MESH equations Naphtali and Sandholm 1971; Ishii and Otto