Function find_root_bisection

pub fn find_root_bisection<T: MyFloat>(
    a: &T,
    b: &T,
    f: impl Fn(&T) -> T,
    epsilon: Fpt,
) -> Option<T>