pub fn newton_eval(xs: &[f64], coef: &[f64], x: f64) -> f64
Evaluate the Newton form (Horner over the nested products) at x.
x