fit_coal_cst function

Maximum likelihood fit of the equilibrium model