s17aef
s17aef
© Numerical Algorithms Group, 2002.
Purpose
S17AEF Bessel function J (x)
0
Synopsis
[fun_value,ifail] = s17aef(x<,ifail>)
Description
This routine evaluates an approximation to the Bessel Function of
the first kind J (x).
0
Note: J (-x)=J (x), so the approximation need only consider x>=0.
0 0
For very large x, it becomes impossible to provide results with
any reasonable accuracy, hence the routine fails. Such arguments
contain insufficient information to determine the phase of
________
/ 2
oscillation of J (x); only the amplitude, / -------, can be
0 \/ (pi)|x|
determined and this is returned on soft failure. The range for
which this occurs is roughly related to the machine precision;
the routine will fail if |x|>~1/machine precision
Parameters
s17aef
Required Input Arguments:
x real
Optional Input Arguments: <Default>
ifail integer -1
Output Arguments:
ifail integer