DOUBLE PRECISION FUNCTION PDA_DBESJ1 (X, STATUS)
  
  
  ***BEGIN PROLOGUE  PDA_DBESJ1
  ***PURPOSE  Compute the Bessel function of the first kind of order one.
  ***LIBRARY   SLATEC (FNLIB)
  ***CATEGORY  C10A1
  ***TYPE      DOUBLE PRECISION (BESJ1-S, PDA_DBESJ1-D)
  ***KEYWORDS  BESSEL FUNCTION, FIRST KIND, FNLIB, PDA_ORDER ONE,
               SPECIAL FUNCTIONS
  ***AUTHOR  Fullerton, W., (LANL)
  ***DESCRIPTION
  
   PDA_DBESJ1(X) calculates the double precision Bessel function of the
   first kind of order one for double precision argument X.
  
   Series for BJ1        on the interval  0.          to  1.60000E+01
                                          with weighted error   1.16E-33
                                           log weighted error  32.93
                                 significant figures required  32.36
                                      decimal places required  33.57
  
      STATUS   Returned error status.
               The status must be zero on entry. This
               routine does not check the status on entry.
  
  ***REFERENCES  (NONE)
  ***ROUTINES CALLED  PDA_D1MACH, PDA_D9B1MP, PDA_DCSEVL, PDA_INITDS, PDA_XERMSG
  ***REVISION HISTORY  (YYMMDD)
     780601  DATE WRITTEN
     890531  Changed all specific intrinsics to generic.  (WRB)
     890531  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)
     910401  Corrected error in code which caused values to have the
             wrong sign for arguments less than 4.0.  (WRB)
     950404  Implement status.  (HME)
  ***END PROLOGUE  PDA_DBESJ1