URL:
<https://savannah.gnu.org/bugs/?55412>
Summary: Bessel function returns incorrect value
Project: GNU Scientific Library
Submitted by: vladk
Submitted on: Mon 07 Jan 2019 05:59:39 PM UTC
Category: Documentation
Severity: 3 - Normal
Operating System:
Status: None
Assigned to: None
Open/Closed: Open
Release: any
Discussion Lock: Any
_______________________________________________________
Details:
If you run
# include <stdio.h>
# include <gsl/gsl_sf_bessel.h>
int main (void)
{
double x = 5.0;
double y = gsl_sf_bessel_j0 (x);
printf ("J0(%g) = %.18e\n", x, y);
return(0);
}
you get
JO(5) = -1.917848549326277019e-01
That is correct answer according to the definition sin(x)/x.
However, the documentation says that the answer is
J0(5) = -1.775967713143382642e-01
_______________________________________________________
Reply to this item at:
<https://savannah.gnu.org/bugs/?55412>
_______________________________________________
Message sent via Savannah
https://savannah.gnu.org/
