pub fn cubic_hermite_spline(
    node_x: &[f64],
    node_y: &[f64]
) -> CubicHermiteSpline