It seems to me that the hankelh1 function ([image: H_\nu(z)]) has a bug for small argument values, for example:
> hankelh1(1,2e-32) 1.9490859162596872e15 - 3.1830988618379058e31im while > besselj(1,2e-32) + im*bessely(1,2e-32) 1.0e-32 - 3.1830988618379067e31im Iam using Version 0.5.0-dev+3667 under Fedora Linux 23. Cheers, E.
