On 10/10/2013 10:33 PM, i2000s wrote: > Hi there, > > I found there is a bug in Bessel functions. For example, the Hankel function besselh >should be defined as a combination of besselj and 1i*bessely, but in Matlab, >this relation does not hold. I suspect, close to the origin, the besselh >and bessely functions may not be correct. Does anyone have a tempary >solution on this? I am using besselh and bessely to do integrals, which always fail for my studies. Thanks. > > Testing Commands: > >>> besselj(1,0.0001)+bessely(1,0.0001)*1i-besselh(1,1,0.0001) > > ans = > > -3.8982e-13
I do not see a bug really, it seems to be just a floating point noise issue? Unless Both Mathematica and Maple and Matlab happen to have the same bug. The result is close to zero.