C-----------------------------------------------------------------------
C IMSL Name:  BSY1/DBSY1 (Single/Double precision version)
C
C Purpose:    Evaluate the Bessel function of the second kind of
C             order one.
C
C Usage:      BSY1(X)
C
C Example 1:
C                                  Declare variables
      INTEGER    NOUT
      REAL       BSY1, VALUE, X
      EXTERNAL   BSY1, UMACH
C                                  Compute
      X     = 3.0
      VALUE = BSY1(X)
C                                  Print the results
      CALL UMACH (2, NOUT)
      WRITE (NOUT,99999) X, VALUE
99999 FORMAT (' BSY1(', F6.3, ') = ', F6.3)
      END
C  BSY1( 3.000) =  0.325