Elsevier Science Home
Computer Physics Communications Program Library
Full text online from Science Direct
Programs in Physics & Physical Chemistry
CPC Home

[Licence| Download | New Version Template] acyv_v1_0.gz(9 Kbytes)
Manuscript Title: A cubic spline interpolation of unequally spaced data points.
Authors: J. Anderson, R.W.B. Ardill, K.J.M. Moriarty, R.C. Beckwith
Program title: CUSPLN
Catalogue identifier: ACYV_v1_0
Distribution format: gz
Journal reference: Comput. Phys. Commun. 16(1979)199
Programming language: Fortran.
Computer: CDC 6600.
Operating system: CDC SCOPE.
RAM: 18K words
Word size: 60
Peripherals: graph plotter.
Keywords: General purpose, Interpolation, Spline, Cubic spline, Natural cubic spline, Curve length, Differentiation.
Classification: 4.10.

Subprograms used:
Cat Id Title Reference
AAUN_v1_0 APLOT CPC 9(1975)85
AAUN_v1_0 0001 POLAR PLOT AND IMPROVEMENTS CPC 15(1978)437

Nature of problem:
Quite often an experiment results in a series of discrete data points inthe interval (x0,xN). The theorist then wishes to interpolate this data at some point x epsilon(x0,xN) or at a large number of points covering the domain (x0,xN) so that an interpolating curve may be drawn. The theorist may also wish to interpolate the discrete data points in order to i)calculate the length of the curve through the data points ii)calculate the area under the curve through the data points or iii)calculate the derivative of the curve through the data points over the domain (x0,xN). The most efficient method of accomplishing this is with the cubic spline interpolation function.

Solution method:
The cubic spline interpolation function matrix equations are solved to give the second derivatives at the N+1 discrete data points. From these the cubic spline interpolation function equations in each interval x esiplon [xi,xi+1] (i=0,1,...,N-1) are then constructed. The interpolation, or the first derivative of the interpolation, can then easily be carried out at the point x epsilon [xO,xN].