Function interpolate 

fn interpolate<R: Scalar>(nodes: &[f64], vals: &[R], x: f64) -> R

Linear interpolation of vals over nodes at x, flat outside, in the mode. The coordinate is passive: the segment and the weight are plain arithmetic and only the node vols carry derivatives.