DOUBLE PRECISION FUNCTION DBESJ0(X)
C***BEGIN PROLOGUE  DBESJ0
C***DATE WRITTEN   770701   (YYMMDD)
C***REVISION DATE  820801   (YYMMDD)
C***CATEGORY NO.  C10A1
C***KEYWORDS  BESSEL FUNCTION,DOUBLE PRECISION,FIRST KIND,ORDER ZERO,
C             SPECIAL FUNCTION
C***AUTHOR  FULLERTON, W., (LANL)
C***PURPOSE  Computes the d.p. Bessel function of the first kind of
C            order zero.
C***DESCRIPTION
C
C DBESJ0(X) calculates the double precision Bessel function of
C the first kind of order zero for double precision argument X.
C
C Series for BJ0        on the interval  0.          to  1.60000E+01
C                                        with weighted error   4.39E-32
C                                         log weighted error  31.36
C                               significant figures required  31.21
C                                    decimal places required  32.00
C***REFERENCES  (NONE)
C***ROUTINES CALLED  D1MACH,D9B0MP,DCSEVL,INITDS
C***END PROLOGUE  DBESJ0