关于项目
8不完全伽马及相关 功能计算

§8.25计算方法

目录
  1. §8.25(i)序列展开
  2. §8.25(ii)正交
  3. §8.25(iii)渐近展开
  4. §8.25(iv)连续分数
  5. §8.25(v)重复关系

§8.25(i)序列展开

尽管§§8.7,8.19(iv),8.21(vi)收敛于的所有有限值z(z),他们是在以下情况下使用很麻烦|z(z)|由于收敛速度慢和取消。对于大型|z(z)|相应的渐近展开式(通常是发散的)。另请参见卢克(1975第101-102页)泰姆牌手表(1994年b).

§8.25(ii)正交

请参阅Allasia和Besenghi(1987年b)用于数值计算Γ(,z(z))来自(8.6.4)通过梯形法则。

§8.25(iii)渐近展开

迪多纳托和莫里斯(1986)描述了一种计算算法P(P)(,x个)(,x个)对于0,x个0、和+x个0来自§中的一致展开8.12.该算法提供14S精度。文中还给出了一个数值反演程序,用于计算x个(具有10S精度),当P(P)(,x个)根据牛顿定律(§3.8(ii)). 另请参见泰姆牌手表(1987,1994年b).

§8.25(iv)连续分数

计算γ(,z(z))Γ(,z(z))通过连分式描述如下琼斯和王位(1985)高茨基(1979年b, §§4.3, 5)。另请参阅雅各布森等。(1986)泰姆牌手表(1996年b,第280页).

§8.25(v)重复关系

涉及不完全伽马函数的展开通常需要生成序列P(P)(+n个,x个),(+n个,x个),或γ(+n个,x个)用于固定n个=0,1,2,一个有效的程序,部分基于递推关系(8.8.5)和(8.8.6),是中描述的高茨基(1979年b,1999).

用于计算的稳定递归格式E类第页(x个)描述了在里面米勒(1960)对于x个>0和整数第页。对于x个>0而且是真的第页看见阿莫斯(1980亿)奇科利等。(1987,1988)。请参阅奇科利等。(1990年)Stegun和Zucker(1974).