Function peroxide::numerical::root::bisection

pub fn bisection<F: Fn(AD1) -> AD1>(
    f: F, 
    interval: (f64, f64), 
    times: usize, 
    tol: f64
) -> Result<f64, RootError>

Bisection method to find root

Usage

extern crate peroxide;
use peroxide::fuga::*;

fn main() -> Result<(), RootError> {
    let x = bisection(f, (0f64, 4f64), 100, 1e-15)?;
    assert!((x - 3f64).abs() < 1e-15);
    Ok(())
}

fn f<T: AD>(x: T) -> T {
    x.powi(2) - x * 2f64 - 3f64
}