Besj0 Besj0_sub Besj1 Besj1_sub Besjn Besjn_sub Besjnd Besjnd_sub Besjnu Besjnu_sub Besjnud Besjnud_sub Besy0 Besy0_sub Besy1 Besy1_sub Besyn Besyn_sub Besynd Besynd_sub Besynu Besynu_sub Besynud Besynud_sub Cbesh Cbesj Cbesy Sbesjn Sbesjn_sub Sbesjnu Sbesjnu_sub Sbesyn Sbesyn_sub Sbesynu Sbesynu_sub

Cbesy() Sub Cbesy ( Z As  Complex, Nu As  Double, N As  Long, Y() As  Complex, Info As  Long, Optional Kode As  Long = 1  ) Sequence of the Bessel functions of the second kind Yν(z) (fractional order) (complex argument) Purpose This routine computes the sequence of the Bessel functions of the second kind Yν(z) for complex 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 the Bessel functions of the second kind. Kode = 1: Yν(z) (ν = Nu to Nu+N-1) Kode = 2: exp(abs(im(z)))*Jν(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. (default = 1) = 1: No scaling. = 2: Returns exponentially scaled values. Reference SLATEC