A computer program has been developed which is capable of calculating phase boundaries in a ternary system with multiple phases when binary thermodynamic data are known. The important features of the program are as follows: (1) The binary data are expressed in the form of Redlich - Kister equation or Legendre polynomial. (2) The ternary free energy of mixing is calculated using three different model equations which are equivalent to Kohler, Toop and the shortest distance composition path equations. (3) The phase boundary is determined by minimizing the ternary free energy of mixing using simplex and hill-climbing methods. The calculated phase diagrams of the Co-Ni-Fe system at 873˚K and 1073˚K and of the Cr-Ni-Fe system at 1373˚K are discussed in comparison with the measured ones.