besselk - Modified K sub ALPHA Bessel functions of the second kind.
besselk(alpha,x) computes K sub ALPHA modified Bessel functions of the second 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 besselk(alpha(j),x(i)).
If ice is equal to 2 exponentialy scaled Bessel functions is computed
besselk(0.5:3,1:4) besselk(0.5:3,1:4,2)