Function find_root 

pub fn find_root<T>(a: T, b: T, eq: RefractionEq<T>, tolerance: T) -> Option<T>
where
    T: RealField + Copy,