Changed gc solver and added arguments to analyze.

I switched the gc.py find_coexistence function to use fsolve. The analyze
function accepts two new arguments: x0, the initial guess for fsolve, and
two_pass, a Boolean flag that tells analyze whether to try to refine the
coexistence point using a second pass.