added condiiton to make sure we only use the real-parameter Bessel funciton in the correct range
Added Bessel function for real parameter from SciPy in order to compute exceptional eigenvalues. I decided to do this because debugging my own routine bewcame tedious...
added some debugging printouts
fixed bug