Find vector x that minimizes the function f(x), constrained within bounds, using the Broyden–Fletcher–Goldfarb–Shanno Bounded (BFGS-B) algorithm.
The missing gradient is evaluated numerically (forward difference).
For more options and diagnostics consider to use
BfgsBMinimizer directly.