第一类贝塞尔函数：Jn（x）  besselj(nu,Z) computes the Bessel function of the first kind, , for each element of the array Z. The order nu need not be an integer, but must be real. The argument Z can be complex. The result is real where Z is positive.

第二类贝塞尔函数：Yn（x）  bessely(nu,Z) computes Bessel functions of the second kind, , for each element of the array Z. The order nu need not be an integer, but must be real. The argument Z can be complex. The result is real where Z is positive.

第三类贝塞尔函数：Hn（x）  Bessel function of the third kind (Hankel function)  besselh(nu,K,Z)

球贝塞尔函数：
function F = sphericalbessel?(n,x)
if  x==0,
     if n==0,
         F = 1;
     else
         F = 0;
    end
else
     F=sqrt(pi/(2*x))*bessel?(n+0.5,x);
end

其中bessel?代表第一类贝塞尔函数或第二类贝塞尔函数