Neuman Edward
Algorithms 70-71 - Determination of a quadratic spline function with given values of the integrals in subintervals