SUBROUTINE INT1D1 (FI, f, xi, x, FX, h, m, mi, cmu1, c1, clmdam, cm, a, c, d)

ARGUMENTS

OUT

FI(MI)  vector of interpolated values                          real 
FX(M)   spline derivatives of the function at points X         real

IN

F(M)    vector of values of the function                       real
XI(MI)  location of the values of the interpolated vector      real
X(M)    location of the values of input vector F               real
H(M)    coordinate differences H(i)=X(i+1)-X(i) i=1,(M-1)      real
M       dimension of vectors F, FX, X, H, A, D AND C           integer
MI      dimension of vectors FI and XI                         integer 
CMU1, C1,
CLMDAM, CM constants specifying end conditions(see notes)      real
A(M)    working vector                                         real
C(M)    working vector                                         real
D(M)    working vector                                         real

DESCRIPTION

EXAMPLE

NOTES

CMU1, C1, CLMDAM, CM  specify the end conditions for the spline.
        The recommended method is to specify the end point slopes,    
        which is done by setting CMU1=CLMDAM=0.0, and
        putting the desired slopes in C1 and CM. For other use,
        see the source code. Subroutines fd1 and fdm, and 
        functions d1 and dn are available to calculate end point 
        slopes. fd1 and fdm use the analytically differentiated 
        form of the cubic lagrange interpolating polynomial, and 
        may be used for generally spaced points. They are used as...
                   call fd1 (gd1,f,h)
                   call fdm (gdm,f,h,m)
        the derivatives are returned in gd1 or gdm. (f,h,m
        as before). The functions d1 and dn, (dn is an entry
        point in d1), are specialized to uniform spacing,
        and are used as ...      gd1 = d1(f,h,m)
                                 gdm = dn(f,h,m)

SEE ALSO