PDA_DBESJ1 - Bessel function of first kind of order one.

Origin

SLATEC / CAMSUN

Implementation Status:

The routine and its subsidiaries will now return an error status as supplied by PDA_XERMSG.

      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