该Matlab源代码实现了Bessel函数，通过实现几个不同的渐近公式的计算，找到贝塞尔函数的初始根，并通过牛顿-拉弗森方法（Newton-Raphson）更新。

贝塞尔函数（Bessel functions），是数学上的一类特殊函数的总称。通常单说的贝塞尔函数指第一类贝塞尔函数（Bessel function of the first kind）。

贝塞尔函数的具体形式随上述方程中任意实数或复数α变化而变化（相应地，α被称为其对应贝塞尔函数的阶数）。实际应用中最常见的情形为α是整数n，对应解称为n 阶贝塞尔函数。

尽管在上述微分方程中，α本身的正负号不改变方程的形式，但实际应用中仍习惯针对α和−α定义两种不同的贝塞尔函数（这样做能带来好处，比如消除了函数在α=0 点的不光滑性）。

贝塞尔函数也被称为圆柱函数或圆柱谐波，因为他们是于拉普拉斯方程在圆柱坐标上的求解过程中被发现的。

牛顿法最初由艾萨克·牛顿在《流数法》（Method of Fluxions，1671年完成，在牛顿死后的1736年公开发表）。约瑟夫·拉弗森也曾于1690年在Analysis Aequationum中提出此方法。牛顿法（Newton's method）又称为牛顿-拉弗森方法（Newton-Raphson method），它是一种在实数域和复数域上近似求解方程的方法。方法使用函数f(x)的泰勒级数的前面几项来寻找方程f(x)=0的根。

The Bessel function cross product arise in problems on an annulus, and the roots of these need to be found to solve certain boundary value problems.

This function uses several different asymptotic formulae to find an initial guess at the roots of the desired Bessel function cross product, then updates via Newton's method. The method works over a large range of radius ratio and order.

The function currently requires the zerobess function.