golden_ratio: Golden Ratio Search

Description

Golden Ratio search for a univariate function minimum in a bounded interval.

Usage

golden_ratio(f, a, b, ..., maxiter = 100, tol = .Machine$double.eps^0.5)

Arguments

Function or its name as a string.

a, b

endpoints of the interval.

maxiter

maximum number of iterations.

tol

absolute tolerance; default sqrt(eps).

...

Additional arguments to be passed to f.

Value

Return a list with components xmin, fmin, the function value at the minimum, niter, the number of iterations done, and the estimated precision estim.prec

Details

`Golden ratio' search for a univariate function minimum in a bounded interval.

Examples

Run this code

f <- function(x) x * cos(0.1*exp(x)) * sin(0.1*pi*exp(x))
golden_ratio(f, 0, 4, tol=10^-10)  # $xmin    = 3.24848329206212
optimize(f, c(0,4), tol=10^-10)    # $minimum = 3.24848328971188