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