Next stepsΒΆ
The algorithmic possibilities for piecewise linear and cubic spline approximation are explored in a different way in Van Loan [VL00]. In that source, a binary search is used to find the interval for evaluating the piecewise polynomial interpolant.
Further details regarding the derivation of the cubic spline equations, with an emphasis on minimizing memory usage, may be found in a number of sources, e.g., Burden and Faires [BF01], Cheney and Kincaid [CK12], and Atkinson and Han [AH04]. Comprehensive theoretical results can be found in de Boor [deBoor78].
On a historical note, Carl de Boor was elected to several National Academies of Science (USA, Poland and Germany, e.g.) and was awarded the (USA) National Medal of Science in 2003 for his work on splines. Splines continue to be important for computer-aided design and computer graphics, among other applications.
An excellent and pragmatic introduction to finite difference methods is by Fornberg [For98]. Numerical integration is a large topic unto itself; one longer introduction to it is by Davis and Rabinowitz [DR14].