C-----------------------------------------------------------------------
C IMSL Name:  BSIS/DBSIS (Single/Double precision version)
C
C Purpose:    Evaluate a sequence of Modified Bessel functions of the
C             first kind with real order and real positive arguments.
C
C Usage:      CALL BSIS (XNU, X, N, BSI)
C
C Example 1:
C                                  Declare variables
      INTEGER    N
      PARAMETER  (N=10)
C
      INTEGER    K, NOUT
      REAL       BSI(N), X, XNU
      EXTERNAL   BSIS, UMACH
C                                  Compute
      XNU = 0.0
      X   = 10.0
      CALL BSIS (XNU, X, N, BSI)
C                                  Print the results
      CALL UMACH (2, NOUT)
      DO 10  K=1, N
         WRITE (NOUT,99999) XNU+K-1, X, BSI(K)
   10 CONTINUE
99999 FORMAT (' I sub ', F6.3, ' (', F6.3, ') = ', F10.3)
      END
C  I sub  0.000 (10.000) =   2815.717
C  I sub  1.000 (10.000) =   2670.988
C  I sub  2.000 (10.000) =   2281.519
C  I sub  3.000 (10.000) =   1758.381
C  I sub  4.000 (10.000) =   1226.491
C  I sub  5.000 (10.000) =    777.188
C  I sub  6.000 (10.000) =    449.302
C  I sub  7.000 (10.000) =    238.026
C  I sub  8.000 (10.000) =    116.066
C  I sub  9.000 (10.000) =     52.319