scipy.optimize.fsolve () internally uses minpack hybrj () which can terminate abnormally even when well within tolerance. This appears to be a hybrj () bug when the tolerance is very tight and the ...