SUBROUTINE PDA_DPLINT (N, X, Y, C, STATUS)
  
  
  ***BEGIN PROLOGUE  PDA_DPLINT
  ***PURPOSE  Produce the polynomial which interpolates a set of discrete
              data points.
  ***LIBRARY   SLATEC
  ***CATEGORY  E1B
  ***TYPE      DOUBLE PRECISION (POLINT-S, PDA_DPLINT-D)
  ***KEYWORDS  POLYNOMIAL INTERPOLATION
  ***AUTHOR  Huddleston, R. E., (SNLL)
  ***DESCRIPTION
  
       Abstract
          Subroutine PDA_DPLINT is designed to produce the polynomial which
       interpolates the data  (X(I),Y(I)), I=1,...,N.  PDA_DPLINT sets up
       information in the array C which can be used by subroutine PDA_DPOLVL
       to evaluate the polynomial and its derivatives and by subroutine
       PDA_DPOLCF to produce the coefficients.
  
       Formal Parameters
       *** All TYPE REAL variables are DOUBLE PRECISION ***
       N  - the number of data points  (N .GE. 1)
       X  - the array of abscissas (all of which must be distinct)
       Y  - the array of ordinates
       C  - an array of information used by subroutines
       STATUS  - Returned error status.
                 The status must be zero on entry. This
                 routine does not check the status on entry.
       *******  Dimensioning Information  *******
       Arrays X,Y, and C must be dimensioned at least N in the calling
       program.
  
  ***REFERENCES  L. F. Shampine, S. M. Davenport and R. E. Huddleston,
                   Curve fitting by polynomials in one variable, Report
                   SLA-74-0270, Sandia Laboratories, June 1974.
  ***ROUTINES CALLED  PDA_XERMSG
  ***REVISION HISTORY  (YYMMDD)
     740601  DATE WRITTEN
     891006  Cosmetic changes to prologue.  (WRB)
     891006  REVISION DATE from Version 3.2
     891214  Prologue converted to Version 4.0 format.  (BAB)
     900315  CALLs to XERROR changed to CALLs to PDA_XERMSG.  (THJ)
     920501  Reformatted the REFERENCES section.  (WRB)
     950403  Implement status.  (HME)
  ***END PROLOGUE  PDA_DPLINT