nag_bessel_i_alpha (s18ejc) returns a sequence of values for the modified Bessel functions or for real , non-negative and .
nag_bessel_i_alpha (s18ejc) evaluates a sequence of values for the modified Bessel function of the first kind
, where
is real and nonzero and
is the order with
. The
-member sequence is generated for orders
when
. Note that
is replaced by
when
. For positive orders the function may also be called with
, since
when
. For negative orders the formula
is used to generate the required sequence.
All constants in the underlying functions are specified to approximately 18 digits of precision. If denotes the number of digits of precision in the floating-point arithmetic being used, then clearly the maximum number of correct digits in the results obtained is limited by . Because of errors in argument reduction when computing elementary functions inside the underlying functions, the actual number of correct digits is limited, in general, by , where represents the number of digits lost due to the argument reduction. Thus the larger the values of and , the less the precision in the result.
Not applicable.
None.