precision loss with bessel functions in GSL bindings

[John W. Eaton]

I just checked in the following change.

  http://hg.savannah.gnu.org/hgweb/octave/rev/32b15d5c3147

Jaroslav, this should probably be applied to the 3.2.x branch as well.

jwe