pub fn polyval(x: f64, p: &[f64], m: i32) -> f64
Evaluation sum(p[i] * x^i, i, 0, N) via Horner’s method. N.B. p is of length N+1.
sum(p[i] * x^i, i, 0, N)