Scipy: how to curve_fit an implicity scalar function?

Question

in a nutshell: is it possible to concatenate scipy.optimize.curve_fit with scipy.optimize.bisect (or fsolve, or whatever) for implicit sca