void spline (coef, x, y, n) double seval(x0, x, y, coef, n) double spldif(x0, x, y, coef, n) double spldif2(x0, x, y, coef, n) double *coef, *x, *y, x0; int n;
seval returns the value of the interpolated function at the ordinate x0.
spldif returns the derivative of the function at the ordinate x0.
This cubic spline fit has been described by Forsythe, Malcolm & Moler in: "Computer Methods for Mathematical Computations", pp. 76-79. (see seval.for)
dark-ages written Josh? sep-90 added spldif2() PJT jul-92 document written, recall last interval PJT