besselzero。m是函数。其他代码文档。
这段代码是有效的abs(顺序)小于或等于146222.16674537213和370030.762407380分别为第一和第二类型的贝塞尔函数。这是非常快和灵活。200万0(0:10000订单和200根/订单)在200秒可以找到标准的笔记本电脑。准确性是默认设置为12位有效数字。
这段代码是基于另一个提交。我固定的缺陷,改善了文档,提高代码质量,提高了代码的速度。如果你有问题,请让我知道。如果我有时间,我将修复它们。
它是如何工作的:
前三根任何阶贝塞尔可以用一个简单的近似方程。这些方程生成使用最小二乘的根从n = 0:10000订单。近似用于启动哈雷的迭代方法。第四和更高的根可以用理解近似根是定期间隔对于一个给定的顺序。一旦发现,第二和第三根间距可以用之间的距离近似第二和第三根。然后哈雷的方法可以应用于精确定位根。因为该算法取决于良好的前三个零的猜测,如果猜很远然后哈雷的方法将收敛到零错误的将随后导致其他零错误。因此,限制放在贝塞尔秩序,abs (n), 146222.16674537213和370030.762407380分别为第一和第二种。如果n是上面指定的这些限制,则会引发一个错误。