1
2      DOUBLE PRECISION FUNCTION DBESY1(X)
3C***BEGIN PROLOGUE  DBESY1
4C***DATE WRITTEN   770701   (YYMMDD)
5C***REVISION DATE  820801   (YYMMDD)
6C***CATEGORY NO.  C10A1
7C***KEYWORDS  BESSEL FUNCTION,DOUBLE PRECISION,ORDER ONE,SECOND KIND,
8C             SPECIAL FUNCTION
9C***AUTHOR  FULLERTON, W., (LANL)
10C***PURPOSE  Computes the d.p. Bessel function of the second kind of
11C            order one.
12C***DESCRIPTION
13C
14C DBESY1(X) calculates the double precision Bessel function of the
15C second kind of order for double precision argument X.
16C
17C Series for BY1        on the interval  0.          to  1.60000E+01
18C                                        with weighted error   8.65E-33
19C                                         log weighted error  32.06
20C                               significant figures required  32.17
21C                                    decimal places required  32.71
22C***REFERENCES  (NONE)
23C***ROUTINES CALLED  D1MACH,D9B1MP,DBESJ1,DCSEVL,INITDS,XERROR
24C***END PROLOGUE  DBESY1
25
26
27