numerics: root and fsolve in scipy.optimize