simulated annealing

Phase_walk now works, uses simulated annealing at each test point to
find the binodal curve.