This is a generic 1-D rational spline package.
Adopted from
Spaeth, H., 1983, Spline-Algorithmen zur Konstruktion glatter
Kurven und Flaechen, Muenchen, R. Oldenbourg
Integration algorithm by Stift (1990)
The procedure Create
establishes the spline coefficients for the
tabulated function Y(X)
The function Point
returns the interpolated value Y(XX)
The function Integrate
returns the integral over the tabulated function
between the X(1) and X(N)
The parameter PP
determines the stiffness of the interpolation.
Its default value is 1.0
This generic package is completely self-contained,
the interface to the outside world is found in the
specification part.
a) With the GNAT compiler
gcc -c -gnatcv Rational_Spline.ads
gcc -c -gnatcv Rational_Spline.adb
b) Usually this package will be WITHed and thus
automatically compiled with the calling unit