Function muller 

pub fn muller(
    f: impl Fn(f64) -> f64,
    x0: f64,
    x1: f64,
    x2: f64,
    tol: f64,
    max_iter: usize,
) -> Option<f64>

Muller’s method — quadratic interpolation, can find complex roots (returns real part here).