Bessel functions only support BigFloat arguments for integer orders

Hi! I'm using the besselj and bessely ( )functions on a function to solve their zeros, as is shown in the source code of boostlibrary (`https://github.com/Ceboc/BesselZeros`). As i want to resolve roots with large values in certain x's interval, x_min < x_roots < x_max, where with .

The problem that i found is that is not consistent how bessely works.

That it's, if i use a i must give an bigfloat argument, but if i want to use a non entire value for , doesn't works.

I was using a wrap of the equivalent boost functions on CxxWrap, but i have problems with the parallelization, so i wrote my own functions to find the zeros of bessels.