Optimize with Fminbox. Can't find the way how to put tolerance and number of iterations

You can use options, this way:

results = optimize(f, lower, upper, initial_x, Fminbox(GradientDescent()), 
Optim.Options(outer_iterations = 2))

See options here: Optim.jl