Function golden_section 

pub fn golden_section(
    f: impl Fn(f64) -> f64,
    a: f64,
    b: f64,
    tol: f64,
) -> (f64, f64)

Golden-section search for a unimodal minimum on [a, b].

Returns (x_min, f_min).