# d9b1mp.f

          SUBROUTINE D9B1MP (X, AMPL, THETA)
    C***BEGIN PROLOGUE  D9B1MP
    ....
    ....
    Warning:   this routine is not intended to be user-callable.
    ....
    ....
    C***SUBSIDIARY
    C***PURPOSE  Evaluate the modulus and phase for the J1 and Y1 Bessel
    C            functions.
    C***LIBRARY   SLATEC (FNLIB)
    C***CATEGORY  C10A1
    C***TYPE      DOUBLE PRECISION (D9B1MP-D)
    C***KEYWORDS  BESSEL FUNCTION, FNLIB, MODULUS, PHASE, SPECIAL FUNCTIONS
    C***AUTHOR  Fullerton, W., (LANL)
    C***DESCRIPTION
    C
    C Evaluate the modulus and phase for the Bessel J1 and Y1 functions.
    C
    C Series for BM1        on the interval  1.56250E-02 to  6.25000E-02
    C                                        with weighted error   4.91E-32
    C                                         log weighted error  31.31
    C                               significant figures required  30.04
    C                                    decimal places required  32.09
    C
    C Series for BT12       on the interval  1.56250E-02 to  6.25000E-02
    C                                        with weighted error   3.33E-32
    C                                         log weighted error  31.48
    C                               significant figures required  31.05
    C                                    decimal places required  32.27
    C
    C Series for BM12       on the interval  0.          to  1.56250E-02
    C                                        with weighted error   5.01E-32
    C                                         log weighted error  31.30
    C                               significant figures required  29.99
    C                                    decimal places required  32.10
    C
    C Series for BTH1       on the interval  0.          to  1.56250E-02
    C                                        with weighted error   2.82E-32
    C                                         log weighted error  31.55
    C                               significant figures required  31.12
    C                                    decimal places required  32.37
    C
    C***SEE ALSO  DBESJ1, DBESY1
    C***REFERENCES  (NONE)
    C***ROUTINES CALLED  D1MACH, DCSEVL, INITDS, XERMSG
    C***REVISION HISTORY  (YYMMDD)
    C   770701  DATE WRITTEN
    C   890531  Changed all specific intrinsics to generic.  (WRB)
    C   890531  REVISION DATE from Version 3.2
    C   891214  Prologue converted to Version 4.0 format.  (BAB)
    C   900315  CALLs to XERROR changed to CALLs to XERMSG.  (THJ)
    C   900720  Routine changed from user-callable to subsidiary.  (WRB)
    C   920618  Removed space from variable name and code restructured to
    C           use IF-THEN-ELSE.  (RWC, WRB)
    C***END PROLOGUE  D9B1MP