pub fn lagrange_interp(xs: &[f64], ys: &[f64], x: f64) -> f64
Lagrange polynomial interpolation at x.