pub fn polyvalx(n: usize, p: &[f64], x: f64) -> f64
Evaluation sum(p[i] * x^i, i, 0, N) via Horner’s method.
sum(p[i] * x^i, i, 0, N)