The amended triple diagonal matrix algorithm and the feasible method of solving fugacity coefficient and activity coefficient were used to calculate the column.