🚀 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.

The (excellent, superb, ++good, highly recommended) book

        G. E. Forsythe, M. A. Malcolm, and C. B. Moler
        "Computer Methods for Mathematical Computations"
        Prentice-Hall, Englewood Cliffs, 1977

discusses cubic splines in detail, and contains a Fortran subroutine
which does just what you want.

	echo 'send spline from fmm' | Mail netlib@ornl.gov

will get you a copy.

- 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}