The file lbfgsb.py in Lib/optimize/ has a bug that prevents it working with an argument list. Here's a patch: 179c179 < f[0], g = func_and_grad(x, *args) --- > f[0], g = func_and_grad(x) Cheers, Ed