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