I am using SciPy’s optimize.curve_fit to get polynomial coefficients for an implicit function, and that was successful. I now have y(x). I need x(y). If I instead do a curve_fit for y, x (as opposed to x, y), the fit is not as good near the origin, and it needs to be as good as possible. I am in no danger of over-fitting.