pub fn interp(x: &[f64], xp: &[f64], fp: &[f64]) -> Vec<f64>