nag_bessel_k1_scaled (s18cdc) evaluates an approximation to , where is a modified Bessel function of the second kind. The scaling factor removes most of the variation in .
The function uses the same Chebyshev expansions as nag_bessel_k1 (s18adc), which returns the unscaled value of .
4 References
Abramowitz M and Stegun I A (1972) Handbook of Mathematical Functions (3rd Edition) Dover Publications
5 Arguments
1:
x – doubleInput
On entry: the argument of the function.
Constraint:
. If x is too close to zero, there is a danger of overflow, and a failure will occur.
2:
fail – NagError *Input/Output
The NAG error argument (see Section 3.6 in the Essential Introduction).
6 Error Indicators and Warnings
NE_REAL_ARG_LE
On entry, x must not be less than or equal to 0.0: .
is undefined and the function returns zero.
NE_REAL_ARG_TOO_SMALL
On entry, x must be greater than : .
The function returns the value of the function at the smallest permitted value of the argument.
7 Accuracy
Relative errors in the argument are attenuated when propagated into the function value. When the accuracy of the argument is essentially limited by the machine precision, the accuracy of the function value will be similarly limited by at most a small multiple of the machine precision.
8 Further Comments
None.
9 Example
The following program reads values of the argument from a file, evaluates the function at each value of and prints the results.