Bessel function scaling + limited range?

[Dmitri A. Sergatskov]

Unless the Matlab compatibility is a must, you can use
a wrapper to GSL's library from octave-forge:


octave:1> bessel_Kn(24, 0.6)
ans =  4.5588e+34
octave:2> bessel_Kn_scaled(24, 0.6)
ans =  8.3068e+34

octave:5> help bessel_Kn_scaled
 -- Loadable Function: Y = bessel_Kn_scaled (N, X)
 -- Loadable Function: [Y, ERR] = bessel_Kn_scaled (...)
     ERR contains an estimate of the absolute error in the value Y.

     This function is from the GNU Scientific Library, see
     `http://www.gnu.org/software/gsl/' for documentation.

/usr/libexec/octave/packages/gsl-1.0.3/x86_64-redhat-linux-gnu-api-v32/gsl_sf.oct


Sincerely,

Dmitri.
--