Abstract: This paper modifies a mathematical model for distillation with two liquid phases so as to suit for thermally coupleddistillation with two liquid phases and advances an algorithm based upon the modified Naphtali-Sandholm method.Thestructures of Jacobi matrix and its submatrices of all types are given.The multiple phase flash method combined with acomparison of free energies for determination of the liquid phase split in stage is used.Seperation of butanol-water-propanol ternary system is employed to justify the method mentioned above.Satisfactory results are obtained.