Ticket #2733: bessel_jn.patch

boost/math/special_functions/detail/bessel_jn.hpp

@@ -32 +32 @@
     {
         return bessel_j0(x);
     }
-    if (n == 1)
+    if ((n == 1) || (n == -1))
     {
-        return bessel_j1(x);
+        factor = n;  // J_{-1}(z) = -1*J_1(z)
+        value = bessel_j1(x);
+        value *=factor;
+        return value;
     }
     if (n < 0)
     {