Function rgsl::polynomials::divided_difference_representation::poly_dd_hermite_init

pub fn poly_dd_hermite_init(
    dd: &mut [f64],
    za: &mut [f64],
    xa: &[f64],
    ya: &[f64],
    dya: &[f64]
) -> Value

This function computes a divided-difference representation of the interpolating Hermite polynomial for the points (x, y) stored in the arrays xa and ya of length size. Hermite interpolation constructs polynomials which also match first derivatives dy/dx which are provided in the array dya also of length size. The first derivatives can be incorported into the usual divided-difference algorithm by forming a new dataset z = {x_0,x_0,x_1,x_1,…}, which is stored in the array za of length 2size on output. On output the divided-differences of the Hermite representation are stored in the array dd, also of length 2size. Using the notation above, dd[k] = [z_0,z_1,…,z_k]. The resulting Hermite polynomial can be evaluated by calling gsl_poly_dd_eval and using za for the input argument xa.