The issue appears in the latest stable scipy (0.12.0) (with numpy 1.6.1 and Python 2.7.3) on a 64 Bit Linux. When using scipy.optimize.root with the default method 'hybr' (or equivalently ...