Besi0 Besi1 Besin Besind Besinu Besk0 Besk1 Beskn Besknd Besknu Besknud Cbesi Cbesk Sbesin Sbesinu Sbeskn Sbesknu

Cbesi() Sub Cbesi ( 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 first kind Iν(z) (fractional order) (complex argument) Purpose This routine computes the sequence of the modified Bessel functions of the first kind Iν(z) for complex argument z and real nonnegative orders ν, ν+1, ... A scaling option is available to help avoid overflow. Parameters [in] Z Argument z. [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 the modified Bessel functions of the first kind. Kode = 1: Iν(z) (ν = Nu to Nu+N-1) Kode = 2: exp(-abs(re(z)))*Iν(z) (ν = Nu to Nu+N-1) [out] Info = 0: Successful exit. = -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 re(Z) too large on Kode = 1. = 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