input value
evaluated Bessel function
var v = y1( 0.0 );
// returns -Infinity
v = y1( 1.0 );
// returns ~-0.781
v = y1( -1.0 );
// returns NaN
v = y1( Infinity );
// returns 0.0
v = y1( -Infinity );
// returns NaN
v = y1( NaN );
// returns NaN
Computes the Bessel function of the second kind of order one.
Notes
x
is very poor. Full accuracy is achieved at1.0e-308
but trends progressively to zero at5e-324
. This suggests that underflow (or overflow, perhaps due to a reciprocal) is effectively cutting off digits of precision until the computation loses all accuracy at5e-324
.