[Boost-users] Boost_1_54_0: boost::math::sph_bessel(n, x) crash at x=0 n>0

Dear developpers, I'am running Boost 1_54_0 on Scientific Linux slc-6.1, gcc-4.4.5, 64 bits machine. I face a problem with computation of spherical bessel function, although the problem seems well known here is a test program { int L; try { for(L=100;L<2010;L++) { boost::math::sph_bessel(L, 0.); } } catch (std::exception& sex) { cout << "Error j_"<0 jL(0) = 0 but 1) I'm using a program which compute integral of product of bessel function and sometimes it needs to call the function at 0. So I would expect the BOOST library to be protected against that. Am I wrong? Best regards Jean-Eric