>>>>> "David" == David Joyner <wdjoyner at gmail.com> writes:
David> Hi:
David> The manual
David> http://maxima.sourceforge.net/docs/manual/en/maxima_16.html#SEC52
David> says i0(z) has been deprecated
David> and to use bessel_i(0,z) instead. However,
David> i0(x) still works (fortunately) since bessel_i(0,x)
David> does not:
David> (%i2) i0(1);
David> (%o2) 1.266065877752009
David> (%i3) bessel_i(0,1);
David> (%o3) bessel_i(0, 1)
David> (%i5) bessel_i(0,1.0);
David> (%o5) bessel_i(0, 1.0)
Hmm. And setting numer doesn't help in this case. And bessel_y and
bessel_k suffer from the same problem.
Fortunately, I think the fix is simple. For some reason the code
checks if the order is strictly positive. If not, then it's not
numerically evaluated. Changing the test to non-negative order fixes
this issue.
I'll check in the change shortly.
Ray