Quadratic spline subroutine package (Q102416)