SUBROUTINE PDA_DPOLVL (NDER, XX, YFIT, YP, N, X, C, WORK, IERR)
  
  
  ***BEGIN PROLOGUE  PDA_DPOLVL
  ***PURPOSE  Calculate the value of a polynomial and its first NDER
              derivatives where the polynomial was produced by a previous
              call to PDA_DPLINT.
  ***LIBRARY   SLATEC
  ***CATEGORY  E3
  ***TYPE      DOUBLE PRECISION (POLYVL-S, PDA_DPOLVL-D)
  ***KEYWORDS  POLYNOMIAL EVALUATION
  ***AUTHOR  Huddleston, R. E., (SNLL)
  ***DESCRIPTION
  
       Abstract -
          Subroutine PDA_DPOLVL calculates the value of the polynomial and
       its first NDER derivatives where the polynomial was produced by
       a previous call to PDA_DPLINT.
          The variable N and the arrays X and C must not be altered
       between the call to PDA_DPLINT and the call to PDA_DPOLVL.
  
       ******  Dimensioning Information *******
  
       YP   must be dimensioned by at least NDER
       X    must be dimensioned by at least N (see the abstract )
       C    must be dimensioned by at least N (see the abstract )
       WORK must be dimensioned by at least 2*N if NDER is .GT. 0.
  
       *** Note ***
         If NDER=0, neither YP nor WORK need to be dimensioned variables.
         If NDER=1, YP does not need to be a dimensioned variable.
  
  
       *****  Input parameters
         ***  All TYPE REAL variables are DOUBLE PRECISION ***
  
       NDER - the number of derivatives to be evaluated
  
       XX   - the argument at which the polynomial and its derivatives
              are to be evaluated.
  
       N    - *****
              *       N, X, and C must not be altered between the call
       X    - *       to PDA_DPLINT and the call to PDA_DPOLVL.
       C    - *****
  
  
       *****  Output Parameters
         ***  All TYPE REAL variables are DOUBLE PRECISION ***
  
       YFIT - the value of the polynomial at XX
  
       YP   - the derivatives of the polynomial at XX.  The derivative of
              order J at XX is stored in  YP(J) , J = 1,...,NDER.
  
       IERR - Output error flag with the following possible values.
            = 1  indicates normal execution
  
       ***** Storage Parameters
  
       WORK  = this is an array to provide internal working storage for
               PDA_DPOLVL.  It must be dimensioned by at least 2*N if NDER is
               .GT. 0.  If NDER=0, WORK does not need to be a dimensioned
               variable.
  
  ***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  (NONE)
  ***REVISION HISTORY  (YYMMDD)
     740601  DATE WRITTEN
     890531  Changed all specific intrinsics to generic.  (WRB)
     891006  Cosmetic changes to prologue.  (WRB)
     891006  REVISION DATE from Version 3.2
     891214  Prologue converted to Version 4.0 format.  (BAB)
     920501  Reformatted the REFERENCES section.  (WRB)
  ***END PROLOGUE  PDA_DPOLVL