fixed a 1/z issue in computing bessel functions