New implementation for Bessel and hankel functions
Bessel and Hankel function implementation have been changed for improved speed. No changes to test or documentation.
No code duplicates are created. Now the coefficients are computed through functions rather than passing functions during the computations. This seems to significantly improve the speed.