关于项目
6指数、对数、正弦和 余弦积分计算

§6.18计算方法

目录
  1. §6.18(i)主要功能
  2. §6.18(ii)辅助功能
  3. §6.18(iii)零点
  4. §6.18(iv)其他参考文献

§6.18(i)主要功能

对于较小或中等值x个|z(z)|幂级数的扩张6.6)或球面贝塞尔函数系列6.10(ii))可以使用。对于大型x个|z(z)|这些系列遭受缓慢的收敛或取消(或两者兼而有之)。然而,这个问题对于球贝塞尔函数系列来说不太严重,因为它们收敛速度更快,并且(除了(6.10.6))在以下情况下未取消z(z)=x个(>0).

对于大型x个|z(z)|,逆阶乘级数的展开式6.10(i))或渐近展开(§6.12)是可用。渐近展开式的可达到精度可以是以指数级的改善显著增加。此外,其他范围的酸碱度z(z)可以通过使用的连续公式来覆盖§6.4.

积分表示的求积是另一种有效的方法。对于例如,高斯-拉盖尔公式(§3.5(v))可以应用于(6.2.2); 看见托德(1954)曾和李(1998)。对于高斯-勒让德公式的应用(§3.5(v))参见图珀和马克(1968).

最后,连分数(6.9.1)可以在以下情况下使用|z(z)|远离原点。z(z)靠近然而,实际轴为负。

§6.18(ii)辅助功能

幂级数、渐近展开式和求积也可用于计算功能(f)(z(z))(z(z))此外,阿克顿(1974)制定了如下复发程序。对于n个=0,1,2,,定义

6.18.1 A类n个 =0e(电子)z(z)1+2(21+2)n个d日,
Bn个 =0e(电子)z(z)1+2(21+2)n个d日,
C类n个 =0e(电子)z(z)(21+2)n个d日.

然后(f)(z(z))=B0,(z(z))=A类0、和

6.18.2 A类n个1 =A类n个+z(z)2n个C类n个,
Bn个1 =2n个Bn个+z(z)A类n个12n个1,
C类n个1 =C类n个+Bn个1,
n个=1,2,,.

A类0,B0、和C类0可以用米勒算法计算3.6(iii)),从初始值开始(A类N个,BN个,C类N个)=(1,0,0)比如,在哪里N个是任意大整数,并通过C类0=1/z(z).

§6.18(iii)零点

的零Ci公司(x个)(x个)可以通过以下方法进行高精度计算牛顿定律(§3.8(ii)),使用渐近展开(6.13.2)作为初始近似值。

§6.18(iv)其他参考文献

对于所处理函数的计算方法的全面调查在本章中,请参见范德兰和特姆(1984,第四章).