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