fn roots<T:std::ops::Add<Output = T> +
std::ops::Sub<Output = T> +
std::ops::Div<Output = T> +
std::ops::Mul<Output = T> +
std::ops::Neg<Output = T> +
num_traits::FromPrimitive +
PartialOrd +
std::fmt::Debug +
Copy + Clone +
num_traits::Zero>(f:Polynomial<T>, left:T, right:T) -> Vec<T>
{
    let rr = base_roots(f.clone(),left,right);
    rr.into_iter().map(|x| { root_by_newton_raphson_near(x, f.clone())}).collect()
}