pub fn cubic_hermite_interpolation<F: IntegrateFloat>(
x: &[F],
y: &[Array1<F>],
dy: &[Array1<F>],
x_new: F,
) -> Array1<F>pub fn cubic_hermite_interpolation<F: IntegrateFloat>(
x: &[F],
y: &[Array1<F>],
dy: &[Array1<F>],
x_new: F,
) -> Array1<F>