besseli - Modified I sub ALPHA Bessel functions of the first kind.
besseli(alpha,x) computes I sub ALPHA modified Bessel functions of the first kind, for real, non-negative order alpha and argument x. alpha and x may be vectors. The output is m-by-n with m = size(x,'*'), n = size(alpha,'*') whose (i,j) entry is besseli(alpha(j),x(i)).
If ice is equal to 2 exponentialy scaled Bessel functions is computed
besseli(0.5:3,1:4) besseli(0.5:3,1:4,2)