Module 828 in TOMSDNSPLIN1: Discrete nonlinear spline interpolation. Given a sequence of data values, constructs a discrete approximation to a nonlinear interpolating spline function. Minimizes the discretized strain energy by a Sobolev gradient descent method. (See R. Renka, ACM TOMS 29 (2003) pp. 458-468.).