asyik.f
SUBROUTINE ASYIK (X, FNU, KODE, FLGIK, RA, ARG, IN, Y)
C***BEGIN PROLOGUE ASYIK
....
....
Warning: this routine is not intended to be user-callable.
....
....
C***SUBSIDIARY
C***PURPOSE Subsidiary to BESI and BESK
C***LIBRARY SLATEC
C***TYPE SINGLE PRECISION (ASYIK-S, DASYIK-D)
C***AUTHOR Amos, D. E., (SNLA)
C***DESCRIPTION
C
C ASYIK computes Bessel functions I and K
C for arguments X.GT.0.0 and orders FNU.GE.35
C on FLGIK = 1 and FLGIK = -1 respectively.
C
C INPUT
C
C X - argument, X.GT.0.0E0
C FNU - order of first Bessel function
C KODE - a parameter to indicate the scaling option
C KODE=1 returns Y(I)= I/SUB(FNU+I-1)/(X), I=1,IN
C or Y(I)= K/SUB(FNU+I-1)/(X), I=1,IN
C on FLGIK = 1.0E0 or FLGIK = -1.0E0
C KODE=2 returns Y(I)=EXP(-X)*I/SUB(FNU+I-1)/(X), I=1,IN
C or Y(I)=EXP( X)*K/SUB(FNU+I-1)/(X), I=1,IN
C on FLGIK = 1.0E0 or FLGIK = -1.0E0
C FLGIK - selection parameter for I or K function
C FLGIK = 1.0E0 gives the I function
C FLGIK = -1.0E0 gives the K function
C RA - SQRT(1.+Z*Z), Z=X/FNU
C ARG - argument of the leading exponential
C IN - number of functions desired, IN=1 or 2
C
C OUTPUT
C
C Y - a vector whose first in components contain the sequence
C
C Abstract
C ASYIK implements the uniform asymptotic expansion of
C the I and K Bessel functions for FNU.GE.35 and real
C X.GT.0.0E0. The forms are identical except for a change
C in sign of some of the terms. This change in sign is
C accomplished by means of the flag FLGIK = 1 or -1.
C
C***SEE ALSO BESI, BESK
C***ROUTINES CALLED R1MACH
C***REVISION HISTORY (YYMMDD)
C 750101 DATE WRITTEN
C 890531 Changed all specific intrinsics to generic. (WRB)
C 891214 Prologue converted to Version 4.0 format. (BAB)
C 900328 Added TYPE section. (WRB)
C 910408 Updated the AUTHOR section. (WRB)
C***END PROLOGUE ASYIK