Fix Bug 83566 - cyl_bessel_j returns wrong result for x>1000 for high orders

[Ed Smith-Rowland]

Michele,

So I think the patch you had plus tests
testsuite/tr1/5_numerical_facilities/special_functions/09_cyl_bessel_j/check_value.cc
testsuite/tr1/5_numerical_facilities/special_functions/11_cyl_neumann/check_value.cc
plus I'd like the new tests also in
 Â  testsuite/special_functions/08_cyl_bessel_j/check_value.cc
 Â  testsuite/special_functions/10_cyl_neumann/check_value.cc
for C++17.

Put that up with a Changelog and see if we're all good.

Thank you for going through this and other math issues.

Ed