In case someone else is interested, I'm pretty much done with brentq_expanding. https://github.com/statsmodels/statsmodels/pull/723#diff-2 It doesn't work as well without specifying anything as I hoped. But specifying soft bounds, works pretty well in my test cases. Josef