[Maxima] bessel_i quadpack bug 5.26.0
woollett at charter.net
Wed Feb 29 15:38:31 CST 2012
On Feb. 29, 2012, Raymond Toy wrote:
>This could probably be fixed by changing how bessel_i(1,%i*x) is computed.
>bessel_i(1,3*%i) has a small real part >even though we know it must be
>For the other bessel functions at %i*x the result is complex, except
>bessel_j. For bessel_j, it looks like we are >careful and return a purely
>Or maybe tell maxima that bessel_i(n,%i*x) = bessel_j(n,x). Maxima doesn't
>seem to know that property. Maybe the >simplifier should honor %iargs to
>do this transformation.
the relation Mathematica has on page
bessel_i(nu,%i*z) = (%i*z)^nu * bessel_j(nu,z) / z^nu
but I haven't checked this numerically.
More information about the Maxima