fn slopes<R: Scalar>(ks: &[f64], ys: &[R]) -> Vec<R>
The divided differences d_i = (y_{i+1} − y_i)/h_i over the knot sequence.
d_i = (y_{i+1} − y_i)/h_i