fugacity-coefficient functions are now working. It appears that pure-component coefficients, multiplied by pressure and mole fraction to give ideal-solution fugacity, will offer the best method for solving phase equilibrium.
branches/jacob/models/johnpye/fprops/test/init_mix
branches/jacob/models/johnpye/fprops/test/init_mix.c