Let's run your scaling test on curve_fit with a full rank linear fit. In [1]: import numpy as np In [2]: from scipy import optimize as opt In [3]: true_p = np.array([3.0, -4.0, 2.0, -6.]) In [4]: x = np.linspace(0, 1, 9) In [5]: fit = lambda x, a, b,...