bugGNU Scientific Library - Bugs: bug #55412, Bessel function returns incorrect...


You are not allowed to post comments on this tracker with your current authentication level.

bug #55412: Bessel function returns incorrect value

Submitted by:  Vladimir Khodygo <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

Mon 07 Jan 2019 05:59:39 PM UTC, original submission:  

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);

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

Vladimir Khodygo <vladk>


No files currently attached


Depends on the following items: None found

Items that depend on this one: None found


Carbon-Copy List
  • -email is unavailable- added by vladk (Submitted the item)



    No changes have been made to this item

    Back to the top

    Powered by Savane 3.5