pub fn brent(f: impl Fn(f64) -> f64, a: f64, b: f64, tol: f64) -> Option<f64>
Brent’s method — superlinear convergence without derivative. Requires f(a)*f(b) <= 0.