Neuman Edward
Algorithm 46 - Determination of an interpolating quadric spline function