[WIP] Fixes the Hankel overflow

[joeydumont]

This PR aims to fix the Hankel functions overflow near abs(z) = 0. However, this first attempt does not work, and even introduces errors when trying to evaluate the Hankel function as J(z)±iY(z).

I should verify that this issue also happens in the original FORTRAN code.