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