🚀 go-pugleaf

In the base article of this thread, whose reference I have alas
misplaced, a request is made for code to fit a cubic spline through
a data set of non-equidistant points.

One other point...  If you need a higher degree of smoothness than
the usual cubic spline, you might consider *Hermite* splines.  The
(excellent, superb, ++good, highly recommended) book

        David Kahaner, Cleve Moler, Stephen Nash
        "Numerical Methods and Software"
        Prentice-Hall   1989   0-13-627258-4

discusses Hermite splines, and suggests using the PCHIP package to
manipulate them.  The book comes with Fortran code on floppy disk;
I suspect PCHIP itself can be had separately from netlib et al, but
I'm not sure just how or where.

- Jonathan Thornburg
  <jonathan@einstein.ph.utexas.edu> or <jonathan@hermes.chpc.utexas.edu>
  University of Texas at Austin / Physics Dept / Center for Relativity
  and (for a few more months) U of British Columbia / {Astronomy,Physics}