Besi0 Besi0_sub Besi1 Besi1_sub Besin Besin_sub Besind Besind_sub Besinu Besinu_sub Besinud Besinud_sub Besk0 Besk0_sub Besk1 Besk1_sub Beskn Beskn_sub Besknd Besknd_sub Besknu Besknu_sub Besknud Besknud_sub Cbesi Cbesk Sbesin Sbesin_sub Sbesinu Sbesinu_sub Sbeskn Sbeskn_sub Sbesknu Sbesknu_sub

Cbesk() Sub Cbesk ( Z As  Complex, Nu As  Double, N As  Long, Y() As  Complex, Info As  Long, Optional Kode As  Long = 1  ) Sequence of the modified Bessel functions of the second kind Kν(z) (fractional order) (complex argument) Purpose This routine computes the sequence of the modified Bessel functions of the second kind Kν(z) for complex nonzero argument z and real nonnegative orders ν, ν+1, ... A scaling option is available to help avoid overflow. Parameters [in] Z Argument z. (Z <> 0) [in] Nu Order of first member of the sequence. (Nu >= 0) [in] N Number of members in the sequence. (N members of orders Nu to Nu+N-1 to be computed) (N >= 1) [out] Y() Array Y(LY - 1) (LY >= N) Sequence of Bessel functions of the second kind. Kode = 1: Kν(z) (ν = Nu to Nu+N-1) Kode = 2: exp(z)*Kν(z) (ν = Nu to Nu+N-1) [out] Info = 0: Successful exit. = -1: The argument Z had an illegal value. (Z = 0) = -2: The argument Nu had an illegal value. (Nu < 0) = -3: The argument N had an illegal value. (N < 1) = -4: The argument Y() is invalid. = -6: The argument Kode had an illegal value. (Kode <> 1 and Kode <> 2) = nz (0 < nz <= N): The last nz components of Y() set to zero due to underflow. = N+1: Precision warning. (Half precision or less because |Z| or Nu+N-1 is large) = N+2: Precision error. (No precision because |Z| or Nu+N-1 is too large) = N+3: Overflow |Z| too small and/or Nu+N-1 too large. = N+4: Algorithmic error. (Termination condition not met) [in] Kode (Optional) A parameter to indicate the scaling option. = 1: No scaling. = 2: Returns exponentially scaled values. Reference SLATEC