Implement BOBYQA (Bound Optimization BY Quadratic Approximations)

[SimonBlanke]

Implement the BOBYQA algorithm, a derivative-free optimization method for problems with bound constraints. It extends NEWUOA by accommodating box constraints on the variables, while still employing a quadratic approximation strategy within a trust region.

For information about the mathematical background see #73.